Unclassified 


uitu  n  Li  COPY 


AD-A225  737  -jq 


REPORT  DOCUMENTATION  PAGE 

— T-*- 


2b.  DECLASSIFICATION  /  DO 


4.  PERFORMING  ORGANIZA 

FR-11 Jul90 


•  I* T I 


jORT  NUMBER(S) 


1b.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  Pal 
distribution  unH«UajraS°  - 

5.  MONITORING  ORGANIZATION IREPORT  NUMBERfS), 


AfOSR-TR-  9  0  0 


6a.  NAME  Of  PERFORMING  ORGANIZATION 

Prometheus  Inc . 


6b  OFFICE  SY/V180L 
(If  applicable) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Air  Force  Office  of 
Scientific  Research 


6c.  ADDRESS  (Gty,  State,  and  ZIP  Code) 

103  Mansfield  St. 
Sharon,  MA  02067 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 

Building  410 

Bolling  AFB  DC  20332-6448 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


Ji 


8b.  OFFICE  SYMBOL 
(If  applicable) 


JIL 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F4  96 2 0-88 -C- 0028 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 


10.  SOURCE  OF  FUNDING  NUMBERS 


\  \ 

\  MM  t  •  ’  ,  .  .  It 

M  "1 _ si _ 

PROGRAM 
ELEMENT  NO. 

•  i  i  j .  -  i 

PROJECT 

NO. 

TASK 

NO 

i  1 

WORK  UNIT 
ACCESSION  NO 

11.  TITLE  (Include  Security  Classification) 

New  Approaches  to  Beamforming, 

Null  Steering,  and 

Filtering  (Unclass 

ified ) 

i2.  personal  author(s)  Barrett,  Benedetto,  Boyd,  Byrnes,  Giroux,  Khraishi,  Miller, 
Newman.  Ostheimer. _ Rov.  Shapiro _ 

13a.  T$^r^fEPORT 

’*™*ffl®»ra90033l 

14.  1^6  OFjRE£ORT  ^r.  Month,  Day) 

15  PAGE  COUNT 

207 

16.  SUPPLEMENTARY  NOTATION 


1  17.  COSATI  CODES  | 

FIELD 

GROUP 

SUB-GROUP 

18.  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

beamforming,  null  steering,  filtering,  shading 
coefficients,  antennas,  convex  programming,  arrays : 
peak  factor,  wavelets,  spectrum  estimation  ,  [  ' 


iLjiecessarv  and  identify  by  block  number) 


19.  ABSTRACT  _ _  _  _  _  _ 

This  final  report  consists  of  ten  sections.  The  first  is  a  report  describing  our  abilities  in  the  co-site  radio  frequency  interference 
(RFI)  problem.  The  second  is  the  paper  Heisenberg  Wavelets  and  the  Uncertainty  Principle,  which  will  be  expanded  upon  and 
submitted  to  the  appropriate  technical  journal.  The  third  is  a  revised  version  of  the  paper  A  Multidimensional  Wiener- Wintner 
Theorem  and  Spectrum  Estimation,  which  will  appear  shortly  in  the  Transactions  of  the  American  Mathematical  Society.  The  fourth 
is  the  expository  paper  Problems  on  Polynomials  with  Restricted  Coefficients  Arising  from  Questions  in  Antenna  Theory,  to 
appear  shortly  in  the  Proceedings  of  the  1989  NATO  Advanced  Study  Institute  on  Fourier  Analysis  and  its  Applications.  T  lie  fifth  is  the 
paper  An  Ideal  Omnidirectional  Transmitting  Array,  and  Optimal  Peak  Factor  Array,  for  Less  Than  Half-wavelength  Spacing, 
to  be  submitted  shortly  to  the  IEEE  Transactions  on  Antennas  and  Propagation.  The  sixth  is  the  revised  paper  Properties  on  ilie  Unit 
Circle  of  Polynomials  with  Unimodular  Coefficients,  which  just  appeared  in  the  Proceedings  of  the  American  Mathematical  Society. 
The  seventh  gives  an  overview  of  awd,  which  is  our  new  method  of  choosing  shading  coefficients,  based  on  a  convex  programming 
approach.  The  eighth  is  the  latest  revision  of  the  paper  A  Specification  Language  Approach  to  Solving  Convex  Antenna  Weight 
Design  Problems,  which  will  be  submitted  shortly  to  the  appropriate  technical  journal.  The  ninth  is  version  3.0  of  the  awd  User  s 
Manual,  incorporating  several  improvements  which  have  been  added  recently.^cvcral  specific  awd  examples  which  were  done  for  the 
Naval  Underwater  Systems  Center  constitute  the  tenth  and  final  section. 

20.  DISTRIBUTION /AVAILABILITY  OF  A8STRACT 


0 UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  □  DTIC  USERS 


21 


ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

_ lJ  * _ CM  1 _ 1  .  i  Ik  MM/I/- 


22b.  TELEPHONE  (Include  Area  Code) 


/-i 

A^-N-fr- 


iU 


IS 


22c  OFFICE  SYMBOL 

/v  /V _ 


OD  FORM  1473, 84  MAR 


83  APR  edition  may  be  used  until  exhausted. 

All  other  editions  are  obsolete.  _  me 

90  08  27  035 


SECURITY  CLASSIFICATION  OF  'HIS  PAGE 


Unclassified 


r 


i  > 


)•  ii';  ; 


I' 


'u 


1  '  '  1  •  M  J  "l  I  I  > 


!l  C*  »!/)■<;  f 


on 


New  Approaches  to  Beamforining,  Null  Steering,  and  Filtering 


e  i :  t  r  - 


Final  Report  11  July  1990 


Prometheus  Inc. 
103  Mansfield  Street 
Sharon,  MA  02067 


Contract  #  F  49620-88-C-0028 


Prepared  by: 


Michael  J.  Barrett,  John  J.  Benedetto,  Stephen  Boyd, 

James  S.  Byrnes,  Andre  Giroux,  Nasser  M.  Khraishi,  Kenneth  S.  Miller 
Donald  J.  Newman,  Gerald  Ostheimer,  Richard  Roy,  Harold  S.  Shapiro 


Submitted  to: 


Dr.  Arje  Nachman 

Mathematics  and  Information  Sciences 
Air  Force  Office  of  Scientific  Research 


Submitted  by: 


James  S.  Byrnes,  President 


f 


New  Approaches  to  Beamforming,  Null  Steering  and  Filtering 


Abstract 

This  final  report  consists  of  ten  sections.  The  first  is  a  report  describing  our  abilities  in  the  co-site  radio 
frequency  interference  (RFI)  problem.  The  second  is  the  paper  Heisenberg  Wavelets  and  the  Uncertainty 
Principle,  which  will  be  expanded  upon  and  submitted  to  the  appropriate  technical  journal.  The  third  is  a 
revised  version  of  the  paper  A  Multidimensional  Wiener-  Wintner  Theorem  and  Spectrum  Estimation ,  which 
will  appear  shortly  in  the  Transactions  of  the  American  Mathematical  Society.  The  fourth  is  the  expository 
paper  Problems  on  Polynomials  with  Restricted  Coefficients  Arising  from  Questions  in  Antenna  Theory,  to 
appear  shortly  in  the  Proceedings  of  the  1989  NATO  Advanced  Study  Institute  on  Fourier  Analysis  and 
its  Applications.  The  fifth  is  the  paper  An  Ideal  Omnidirectional  Transmitting  Array,  and  Optimal  Peak 
Factor  Array,  for  Less  Than  Half-wavelength  Spacing,  to  be  submitted  shortly  to  the  IEEE  Transactions  on 
Antennas  and  Propagation.  The  sixth  is  the  revised  paper  Properties  on  the  Unit  Circle  of  Polynomials  with 
Unimodular  Coefficients,  which  just  appeared  in  the  Proceedings  of  the  American  Mathematical  Society. 
The  seventh  gives  an  overview  of  awd,  which  is  our  new  method  of  choosing  shading  coefficients,  based  on 
a  convex  programming  approach.  The  eighth  is  the  latest  revision  of  the  paper  A  Specification  Language 
Approach  to  Solving  Convex  Antenna  Weight  Design  Problems,  which  will  be  submitted  shortly  to  the 
appropriate  technical  journal.  The  ninth  is  version  3.0  of  the  awd  User’s  Manual,  incorporating  several 
improvements  which  have  been  added  recently.  Several  specific  awd  examples  which  were  done  for  the 
Naval  Underwater  Systems  Center  constitute  the  tenth  and  final  section.  The  report  concludes  with  the 
Appendix,  which  contains  various  miscellaneous  information. 
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1  Introduction 

This  document  is  to  inform  the  clients  of  Prometheus  Inc.  of  our  interest 
and  abilities  in  the  co-site  radio  frequency  interference  (RFI)  problem.  This 
problem  has  arisen  due  largely  to  the  tremendous  difficulty  and  concomitant 
lack  of  experience  and  practice  in  the  simulation  of  mutual  interference  be¬ 
tween  transmitters  and  receivers  caused  by  system  nonlinearities.  A  special 
simulation  technology  has  been  developed  to  handle  this  problem. 

Given  a  particular  transceiver  (co-sited  transmitter  and  receiver)  and  set 
of  transmitter  waveforms  as  well  as  a  mathematical  description  of  nonlinear¬ 
ities  present,  the  simulation  technology  developed  predicts  the  accumulated 
RFI  in  specified  receiver  bandwidths.  This  is  obviously  useful  in  the  design 
of  transceivers  to  minimize  RFI  due  to  system  nonlinearities.  In  situations 
where  the  design  of  the  system  is  fixed  and  RFI  reduction  is  of  importance, 
the  simulation  technology  provides  estimates  of  the  dominant  sources  of  RFI. 
This  is  valuable  information  in  attempting  to  arrive  at  cost  effective  solutions 
to  the  problem.  The  actual  solutions  are,  of  course,  problem  dependent  and 
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can  incorporate  temporal  adaptive  filtering  techniques  as  well  as  sensor  array 
processing  techniques  that  take  advantage  of  spatial  diversity  of  the  various 
RFI  sources. 


2  Background 

There  are  many  practical  applications  in  which  the  ability  to  simultaneously 
transmit  and  receive  RF  signals  using  a  single  antenna  system  is  desirable. 
Ideally,  though  physically  co-located,  the  two  functions,  transmission  and 
reception,  can  be  separated  in  the  frequency  domain  so  as  to  prevent  mutual 
interference.  However,  nonlinearities  in  the  output  stages  of  the  transmitters 
and  in  the  input  stages  of  the  receivers,  and  nonlinear  mixing  effects  in  the 
platform  environment  can  lead  to  co-channel  interference.  These  nonlinear 
problems  are  difficult  for  most  engineers  to  understand  and  predict,  and  are 
particularly  difficult  for  those  engaged  in  writing  EMI  specifications. 

Co-site  interference  is  associated  primarily  with  several  nonlinear  effects. 
Undesirable  nonlinear  mixing  of  signals  from  different  sources  can  occur  in 
the  skin  of  the  platform  and  in  the  front-ends  of  receivers.  Broadband  self¬ 
noise  associated  with  switching  transients  in  the  frequency  synthesizers  of 
spread-spectrum  systems  is  mixed  with  itself  in  the  output  stages  of  trans¬ 
mitters,  since  microwave  amplifier  output  stages  can  often  be  quite  nonlinear. 
These  effects  cause  the  received  signal-to-interference  ratio  to  deteriorate  as 
the  transmitted  signal  passes  through  the  output  stage.  In  addition,  in¬ 
termodulation  products  and  noisy  harmonics  generated  in  nonlinear  mixing 
from  various  external  sources  can  appear  in  the  wrong  bands,  causing  further 
deterioration  of  the  performance  of  co-sited  systems. 

It  is  generally  accepted  that  the  current  difficulty  in  understanding  and 
coping  with  co-site  interference  between  transmitters  and  receivers  is  partly 
due  to  a  lack  of  simulation  models  for  the  calculation/prediction  of  intermod¬ 
ulation  and  harmonic  interference.  In  addition  to  the  inherent  mathematical 
complexity  of  nonlinear  systems,  the  generation  of  efficient  computer  models 
for  these  systems  is  made  even  more  difficult  by  the  very  great  bandwidths 
and  frequency  separations  involved  in  EMI  between  such  systems. 


3  Predicting  Nonlinear  Effects 

Volterra- Wiener  analysis  [1]  is  probably  the  be:t-known  technique  for  han¬ 
dling  nonlinear  systems,  although  it  has  not  lived  up  to  its  promise  as  a 
means  of  predicting  nonlinear  effects  in  communication  systems.  One  of  t1-  i 
major  difficulties  is  that  Volterra  transfer  functions  [VTFs]  (or  Volterra  ker¬ 
nels)  for  nonlinear  amplifiers  'vith  memory  grow  in  dimension  with  the  order 
of  each  term  in  a  power  series.  That  is,  the  second-order  VTF  is  a  two- 
dimensional  function  (a  complex- valued  function  on  a  two-dimensional  real 
frequency  domain),  the  third-order  VTF  is  three-dimensional,  etc..  Due  to 
the  exponentially  increasing  complexity,  Volterra  series  are  generaiiy  trun¬ 
cated  at  orders  of  around  four  or  five. 

In  general,  the  computational  complexity  involved  in  approximating  broad¬ 
band  nonlinear  systems  with  fading  memory  by  Volterra  series  is  prohibitive. 
For  memoryless  or  instantaneous  nonlinearities,  however,  the  computational 
complexity  is  greatly  reduced.  The  resulting  approximation  can  be  viewed 
as  a  Taylor  series  expansion  of  the  nonlinearity  at  the  origin.  As  a  conse¬ 
quence,  the  issue  of  convergence  of  the  series  must  be  addressed.  For  example, 
non-analytic  memoryless  nonlinearities  do  not  have  Volterra  series  represen¬ 
tations.  The  nonlinearities  must  first  be  approximated  by  analytic  functions, 
adding  another  degree  of  complexity  and  approximation. 

The  new  technique  proposed  herein  involves  expansion  of  the  instanta¬ 
neous  nonlinearities  in  a  set  of  polynomials  orthogonal  with  respect  to  the 
complex  Gaussian  measure.  It  differs  from  Volterra  expansions  most  sig¬ 
nificantly  in  that  analytic  functions  are  not  required;  tabular  I/O  charac¬ 
teristics  of  the  instantaneous  nonlinearities  can  be  handled.  When  applied 
to  a  single  amplifier,  the  new  technique  also  differs  considerably  from  the 
analysis  of  Fuenzalida  et.  al.  [2],  wherein  nonlinear  amplifier  characteris¬ 
tics  are  expanded  in  a  series  of  Bessel  functions.  But  it  has  in  common 
with  that  analysis  the  ability  to  accurately  model  both  the  amplitude  mod¬ 
ulation/amplitude  modulation  [AM/AM]  and  amplitude  modulation/phase 
modulation  [AM/PM]  characteristic  curves.  In  fact,  these  may  be  input  to 
the  simulation  in  equal  decibel  increments  rather  than  in  volts. 

The  alternative  approach  developed  by  Prometheus  Inc.  is  specifically  tai¬ 
lored  to  the  RFI  and  nonlinear  amplifier  problem.  It  incorporates  engineering 
models  designed  for  the  purpose  of  generating  intermodulation  effects  in  a 
computer  representation  which  is  suitable  for,  and  compatible  with,  existing 
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computer  simulations  of  military  communication  and  radar  systems.  The 
simulations  can  accurately  predict  the  interference  effect  of  one  system,  or 
group  of  systems,  upon  another. 

The  most  difficult  problem  in  simulating  communication  systems  is  the 
sample  rate  required.  For  instance,  a  typical  commercial  satellite  transponder 
has  a  bandwidth  of  40  MHz,  which  requires  complex  sampling  at  20  MHz  to 
prevent  aliasing  after  down-conversion  to  baseband.  Now  consider  simulating 
a  combination  of  two  such  transponders  separated  by  100  MHz.  It  would  be 
unthinkable  to  simulate  these  systems  in  the  actual  combined  bandwidths 
they  occupy,  as  this  would  require  sample  rates  in  excess  of  100  MHz.  Where 
multiple  signal  spread-spectrum  systems  are  involved,  the  overall  spread- 
spectrum  bandwidth  of  a  given  system  is  generally  too  great  to  consider 
simulating  directly.  Such  systems  are  usually  simulated  in  the  de-spread 
baseband  bandwidth  of  one  signal,  calculating  where  necessary  the  effects  of 
the  spreading  and  de-spreading  processes  on  the  other  signals  involved. 

The  new  simulation  technique  allows  the  generation  of  all  of  the  inter¬ 
mediation  products  that  fall  into  a  given  band  from  a  group  of  systems  in 
several  widely  spaced  frequency  bands.  The  simulations  can  be  carried  out 
in  sample  rates  appropriate  to  the  transmission  bandwidths  of  each  system, 
referred  to  baseband  in-phase  and  quadrature  [IP&Q]  samples. 


4  Example 

By  way  of  introduction  to  the  new  technique,  let  us  suppose  that  a  metal 
fastener,  bolt,  or  rivet  on  the  surface  of  the  communications  platform  has 
become  an  undesirable  RF  radiator  due  to  a  metallic  interface  that  acts 
like  a  half-wave  rectifier.  The  rectifier  causes  mixing  of  signals  from  two 
transmitters,  resulting  in  intermodulation  products  that  fall  into  the  band  of 
a  co- located  receiver.  An  adequate  model  for  the  purpose  at  hand  would  be  to 
sum  the  two  signals  and  pass  the  sum  through  an  instantaneous  nonlinearity, 
the  actual  characteristics  of  which  could  be  determined  experimentally. 

For  the  purposes  of  computer  simulation,  adding  the  two  signals  together 
at  RF  and  performing  the  indicated  operation  is  not  practical  because  the 
bandwidth  and  associated  sample  rates  are  astronomical.  It  is  preferable  to 
carry  out  simulations  in  the  baseband  bandwidth  of  each  interfering  system, 
and  to  generate  intermodulation  products  as  they  appear  in  the  baseband  of 
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the  system  that  is  being  interfered  with.  Thus,  data  are  required  only  in  the 
conventional  IP&Q  baseband  format  (compatible  with  existing  simulations 
of  the  systems). 

In  the  hypothetical  problem,  there  are  two  inputs  and  one  output.  The 
output  band  of  interest  is  only  a  small  fraction  of  the  entire  bandwidth  of 
intermodulation  products  resulting  from  the  nonlinear  mixing  of  the  two 
signals.  Furthermore,  the  actual  carrier  frequencies  are  not  present  in  the 
data  (since  all  data  is  referenced  to  baseband).  Therefore,  a  means  of  keeping 
track  of  the  sum  and  difference  frequencies  between  carriers  that  set  up  the 
entire  pattern  of  intermodulations  is  required. 

The  nonlinear  computer-aided  analysis  and  simulation  technique  that  was 
recently  developed  [3]  for  the  analysis  and  simulation  of  nonlinear  effects  in 
microwave  amplifiers  satisfies  these  requirements.  Any  complex-valued  (in¬ 
stantaneous)  nonlinear  function  of  a  complex  input  or  inputs  can  be  repre¬ 
sented,  and  the  carrier  frequencies  factored  out ,  so  that  mixing  between  two 
close,  widely  spaced,  or  overlapping  radio  frequency  bands  can  be  done  with 
complex  baseband  IP&Q  samples.  There  are  no  rule-of-thumb  approxima¬ 
tions  involved,  and  computations  may  be  carried  out  to  any  desired  degree 
of  accuracy.  The  technique  is  limited  only  by  the  degree  of  accuracy  of  the 
model  of  the  nonlinearity.  It  can  be  readily  demonstrated  that  all  microwave 
nonlinear i Lies  of  interest  here  are  sufficiently  instantaneous  to  satisfy  the 
fundamental  requirement.  The  technique  will  produce  output  for  only  the 
desired  band  of  the  system  that  is  being  interfered  with,  without  the  need 
for  generating  intermodulation  products  that  are  outside  the  bandwidth  of 
interest. 

The  new  technique  is  based  upon  a  nonlinear  input-output  theory  de¬ 
veloped  by  Price  [4].  Like  the  Vol terra- Wiener  analysis,  it  was  originally 
intended  for  random  inputs,  but  is  suitable  for  deterministic  inputs.  It  uses 
expansions  of  the  nonlinear  instantaneous  transfer  characteristics  in  complex 
Hermite  polynomials,  rather  than  in  a  power  series  in  the  input  like  the  VTFs 
and  Wiener  G-functionals.  A  brief  discussion  of  complex  Hermite  polyno¬ 
mials  is  given  in  the  Appendix.  Again,  the  interest  here  is  not  in  Gaussian 
random  process  inputs,  but  rather  in  actual  non-Rayleigh  envelope  signals. 
Unlike  previous  traveling  wave  tube  [TWTA]  analyses,  there  is  no  attempt 
here  to  determine  output  statistical  properties.  The  main  objective  is  to 
provide  a  detailed  simulation  of  the  actual  signals. 


5  Conclusion 


The  Prometheus  team  has  already  successfully  applied  complex  Hermite 
polynomial  expansions  to  microwave  amplifier  saturation  characteristics,  and 
demonstrated  excellent  curve  fitting  over  the  20  dB  range  of  the  Rayleigh 
decibel  density.  In  that  application,  actual  saturation  characteristic  curves 
were  fit  with  complex  Hermite  polynomials  of  orders  less  than  ten.  Close 
phase  shift  and  gain  agreement  between  the  polynomial  approximations  and 
the  corresponding  actual  AM/AM  and  AM/PM  saturation  curves  were  achieved. 
The  same  technique  can  be  used  with  more  general  nonlinearities,  such  as 
microwave  rectifiers,  and  will  accurately  predict  the  outputs  with  nothing 
other  than  simple  numerical  truncation.  As  much  accuracy  as  desired  can  be 
implemented. 


6  Discussion 

The  capability  described  in  this  document  provides  a  computationally  fea¬ 
sible  means  for  predicting  the  effects  of  various  instantaneous  nonlinearities 
on  the  performance  of  co-sited  RF  communication  systems.  It  is  a  valuable 
tool  in  the  design  and  analysis  of  such  systems.  In  the  design  ctage,  the  sus¬ 
ceptibility  of  the  performance  of  a  particular  communication  system  design 
to  various  types  of  nonlinearities  can  be  predicted.  An  area  for  future  work 
would  be  the  inclusion  of  this  capability  into  an  automated  design  facility, 
where  optimum  system  parameters  minimizing  performance  degradation  are 
sought. 

In  the  system  analysis  and  debugging  phase  of  communication  system  de¬ 
sign,  the  tool  is  capable  of  predicting,  given  several  nonlinearities,  those  that 
are  contributing  most  prominently  to  the  performance  degradation.  Conse¬ 
quently,  those  nonlinearities  that  should  be  eliminated  first  in  an  attempt 
to  improve  performance  can  be  identified.  Similarly,  the  tool  can  be  used 
to  identify  nonlinearities  that  are  present  in  a  particular  physical  system  in 
the  experimental  phase.  Inputs  can  be  designed  to  isolate  certain  types  of 
nonlinearities,  and  actually  applied  to  the  system.  The  presence  of  certain 
components  in  the  output  can  then  be  attributed  to  the  particular  nonlin¬ 
earity  in  question. 

Having  isolated  the  potential  sources  of  the  performance  degradation, 
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elimination  of  the  problem  can  then  be  addressed.  In  the  situation  where  an 
array  of  receivers  is  present  and  the  significant  IIP  nonlinearities  are  point 
source  radiators,  recently  developed  techniques  in  the  area  of  sensor  array 
signal  processing  can  be  employed  to  significantly  reduce  the  efTect  of  the 
nonlinearities.  The  basic  idea  involves  spatially  isolating  the  disturbance 
source,  and  employing  a  transformation  of  the  sensor  array  outputs  (receiver 
inputs)  that  eliminates  the  disturbance. 


A  Mathematical  Development 

In  the  simple  example  described  previously,  a  half-wave  rectifier  causes  mix¬ 
ing  of  signals  from  two  transmitters,  resulting  in  intermodulation  products 
that  fall  into  the  band  of  a  co-located  receiver.  The  appropriate  mathemat¬ 
ical  model  involves  passing  the  sum  of  the  two  signals  through  an  instanta¬ 
neous  nonlinearity,  the  actual  characteristics  of  which  could  be  determined 
experimentally. 


Si  +  i2  =  z1ejuJlt  +  z2ejlM7i  (1) 

The  tilde  represents  radio  frequency  (RF)  signals  and  zx  and  z2  are  the 
complex  baseband  signals.  Only  the  pre-envelope  or  positive-frequency  rep¬ 
resentations  of  the  RF  signals  are  required.  The  RF  signal  is  now  rectified 
by  an  instantaneous  nonlinearity: 

V  —  f{z  1  +  52)  (2) 

The  instantaneous  nonlinearity  is  assumed  to  be  a  complex-valued  func¬ 
tion,  f(z),  of  a  complex  variable.  For  example,  an  ideal  half-wave  rectifier  is 
described  by 


m 


S,  Re  S  >  0 
0,  Re  S  <  0 


(3) 


in  these  representations.  Note  that  the  function  need  not  be  analytic. 

The  basic  idea  is  to  expand  the  nonlinear  function  in  a  series  of  complex 
Ifcrmite  polynomials 


n  ,m 


(•l) 
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Complex  Mermite  polynomials  appear  very  infrequently  in  the  literature, 
and  thus  are  not  wcil  known,  even  to  mathematicians.  They  are  closely 
associated  with  the  complex  normal  density  in  the  same  way  that  ordinary 
Mermite  polynomials  arc  associated  with  the  normal  density.  In  each  case,  a 
probability  density  is  the  weight  function  of  a  set  of  orthogonal  polynomials. 
The  use  of  the  complex  Mermite  polynomials  versus  another  suitable  set  of 
orthogonal  complex  polynomials  is  somewhat  arbitrary  in  this  application, 
although  they  can  be  shown  to  have  suitable  curve-fitting  properties  for  the 
problem  at  hand,  as  they  do  for  TWTA  characteristics. 

Noting  that  the  complex  Mermite  polynomials  are  orthogonal  polynomials 
with  a  weight  function  which  is  a  (unit- variance)  complex  normal  probability 
density,  the  coefficients  in  the  expansion  may  be  determined  by: 

I<n,m  =  c(n,m)  J  exp(  —  ^\z\2)f(z)Hn  m(zm,  z)dz,  (5) 

where  c(n,  m)  is  a  normalizing  factor,  and  it  is  understood  that  the  differential 
elements  are  the  real  and  imaginary  parts  of  2  rather  than  the  complex 
variable  itself.  The  integral  is  actually  a  two  dimensional  integral  in  real 
variables. 

The  complex  Mermite  polynomials  form  a  pyramid  of  the  type  as  indi¬ 
cated  in  Figure  1  with  those  along  the  outside  border  of  the  pyramid  of  the 

1 

2  2* 


Figure  1:  Pyramid  of  Complex  Hermite  Polynomials 

form  5n  or  2*" .  The  internal  polynomials  lead  with  the  term  2n2*m  which 
is  indicated  by  that  particular  polynomial’s  position  on  the  pyramid.  Every 
term  of  a  particular  polynomial  will  retain  the  difference  between  the  num¬ 
ber  of  powers  of  2  and  2*  that  the  leading  term  has.  The  expansion  of  an 
amplitude  sensitive  nonlinearity,  like  the  TWTA,  uses  only  those  polynomi¬ 
als  along  the  diagonal  for  which  there  is  one  more  power  of  2  than  of  2". 
This  occurs  because  a  purely  amplitude  sensitive  nonlinearity  does  not  gen¬ 
erate  harmonics  [5],  but  only  intcrmodulations  in  or  near  the  original  band 


I-  8 


of  transmission.  The  half-wave  rectifier ,  on  the  other  hand,  has  a  general 
expansion  using  terms  from  all  of  the  diagonals. 

The  carrier  frequencies  can  now  be  factorcd-out,  and  the  output  can  be 
obtained  in  only  a  desired  band.  In  any  particular  Hermite  polynomial,  the 
difference  between  the  number  of  powers  of  2  and  of  its  complex  conjugate 
2*  is  fixed.  Using  (1),  (2),  and  the  binomial  theorem: 


*)  =  Ecoeff-x(zieJUl'  +  2^T(vM!  +  ^e-«Om  (6) 

X 

=  2>eir.  x  £  ("  W-y(  V"11)"-''  m 

1  iA  w 

*  £  (8) 

The  carrier  frequencies  ux  and  a>2  can  be  factored-out  with  their  various 
multiplicities,  leaving  only  the  baseband  signals  zx  and  z2,  with  their  inter¬ 
modulations,  for  the  numerical  simulation  to  contend  with.  But,  by  keeping 
track  of  the  carrier  frequencies  as  a  side  problem,  it  is  readily  evident  where 
the  intermodulations  belong  in  the  frequency  spectrum.  A  particular  term 
would  be,  for  instance,  of  the  form: 

ziz2 ~vz[qZ2m~q  exp[/t[(pu;1  +  (n  -  p)u>2  -  qu  1  -  (m  -  <7)^2]]  (9) 

The  sums  and  differences  of  the  carrier  frequencies  may  be  accounted  for, 
and  the  frequency  of  the  intermodulation  band  specified.  Only  the  desired 
band  of  intermodulation  products  need  be  simulated.  The  corresponding 
coefficients  in  the  complex  Hermite  polynomial  expansion  may  be  determined 
in  advance  of  the  simulation. 
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Section  II 

Heisenberg  Wavelets 


Heisenberg  wavelets  and  the  uncertainty  principle 

Jolm  J.  Benedetto* 

Prometheus  Inc. 

For  the  past  few  years  the  term  “wavelet”  has  indicated  a  function  rp  G  T2(5?)  for  which  the 
set  {V’m.n}  of  dilations  and  translations,  associated  with  the  affine  group,  is  an  orthonormal 
basis  of  Z/2(3?).  In  light  of  the  theory  of  wavepackets  and  the  sustained  interest  in  Weyl- 
Heisenberg  decompositions  because  of  the  importance  of  frequency  localizaton,  we  refer  to 
both  “affine  wavelets”  ip  and  “Heisenberg  wavelets”  <p  G  L2(3£),  where 

<f>m,n(t)  =  e2*'mbt<Kt-na)  (1) 

for  fixed  a,  b  >  0.  The  term  “Heisenberg”  is  used  because  of  the  underlying  importance  of 
the  Heisenberg  group  in  analyzing  (1).  Elementary  properties  of  {ipmtn}  and  {<£m,n}  are 
given  in  [Ben89,  HW89].  This  note  concerns  Heisenberg  wavelets. 

The  construction  of  Heisenberg  wavelets  is  a  traditional  goal  in  signal  processing,  e.g., 
Gabor’s  classical  paper  [Gab46].  The  idea  originated  in  quantum  mechanics  where  “over 
complete”  sets  (coherent  states)  are  associated  to  objects  such  as  the  Weyl-Heisenberg  group 
or  the  von  Neumann  lattice.  Related  devices  include  the  Wigner  distribution  (1932)  and  the 
ambiguity  function.  The  Wigner  distribution  became  a  tool  in  signal  analysis  in  the  1940’s 
(by  Ville)  and  presently  plays  an  important  role  in  a  host  of  electrical  engineering  problems 
[CM80].  The  ambiguity  function  is  a  staple  in  optics  and  high  resolution  radar. 

1  The  uncertainty  principle  inequality. 

The  uncertainty  principle  inequality  is 

ll/lll  <  4<r|l((  -  fo)/Wlhll(7  -  70)7(7)112.  (1.1) 

The  Fourier  transform  f  is  defined  as 

7(7)  =  /  f(t)e~2wit’1  dt, 

where  t  6  5?,  7  €  9?  (=3?),  and  integration  is  over  5?.  Our  notation  is  standard,  e.g.,  [SW71], 
and  (1.1)  is  the  simplest  of  a  whole  class  of  uncertainty  principle  inequalities,  e.g.,  [Ben90]. 

There  is  equality  in  (1-1)  in  the  case  /  is  the  Gaussian.  Because  of  this,  one  of  Gabor’s 
ideas  from  [Gab46]  was  to  write  discrete  expansions  of  L2(5?)  elements  in  terms  of  {<pm,n}, 
where  <j>  is  the  Gaussian.  This  part  of  his  program  is  subject  to  some  scrutiny  because  of  the 
Balian-Low  theorem:  If  {d’m.n}  is  a  frame  for  <P  G  T2(3?)  and  ab  =  1  then  either  t<f>(t)  L2(9?) 

or  7^(7)  i  L2(%)  [BHW]. 


*  The  author  is  also  Professor  of  Mathematics  at  the  University  of  Maryland,  College  Park, 
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2  Weak  uncertainty. 

Given  tp  G  L2( 3£)  for  which  \\tp\\2  =  1.  It  is  well  known  and  easy  to  verify  that 

IK*  ~  *m,n)0m,n(*)||2||(7  “  7m,»)0m,n(7)l|2  =  IK*  ~  *0 .0 )  0  (  *  )  1 1 2 1 1(7  ~  T'O  0  )  0  (  7 )  1 1 2 ,  (2.1) 

where  the  expected,  values  (*m,n,7m,n)  are 

*m,n  =  J  t\4>m,n{t)\2  dt  and  7 m.n  =  J  7|0m,n(7)|2  <*7^ 

e.g.,  [Ben90]. 

An  analogous  result  is  true  for  Heisenberg  wavelets — 

Proposition  2.1.  Given  <f>  G  L2( 3?)  for  which  ||<^||2  =  1  and  let  ab  =  1.  Assume  that 
t<p(t),  7^(7)  G  L2  and 

J  t\<p(t)\2  dt,  J  7|^(7)|2^7  €  3?. 

Define  the  variances 

n  =  IK*  ~  *m,n)^m,n(*)|l2  an<^ 

*^m,n  =  11(7  —  7m,n)^m,n  (7)112 

Then 

°<  Am,n^m,n=  *2|<A(*)P  <**  ~  (/  *|^(*)|2^  j  72|^(7>|2  <*7  ~  7^(7)|2  <*7)  j- 

3  Strong  uncertainty  and  complex  analysis. 

Proposition  2.1  allows  us  to  conclude  that 

sup  Ax2miHA p^n  <  00. 

m,n 

On  the  other  hand,  if  the  expected  values  (tm,nilfm,n)  are  replaced  by  a  fixed  pair  (*o,  70), 
then  it  is  easy  to  see  that  we  have  the  “strong”  uncertainty  property, 

sup  ||(*  ~  *o)^m,n(<)||2||(7  -  7o)^m,n(7)ll2  =  CX>.  (3.1) 

m,n 

The  terminology,  “weak”  and  “strong”  uncertainty,  is  discussed  in  [Ben90]  but  is  not  essential 
to  our  discussion  here.  (In  fact,  at  the  risk  of  making  an  apparently  confusing  remark,  (3.1) 
is  really  a  weak  form  of  strong  uncertainty  since  it  is  valid  only  because  of  translation  from 
the  origin,  not  because  of  internal  changes  within  each  <j)m,n •) 

In  this  section  we  propose  to  give  a  complicated  proof  of  (3.1)  for  dealing  with  refinements 
of  (3.1)  in  which  the  sup  is  replaced  with  various  paths  in  Z  x  Z. 

Lemma  3.1.  Suppose  f  is  analytic  in  the  tube  3?  x  (0,  b ),  /  has  nontangential  limits  a.e.  on 
3?,  and  f  =  0  on  K  C  3?  for  which  \K\  >  0.  Then  f  is  the  zero  function. 

In  this  generality,  D.  Hamilton  has  pointed  out  that  topological  arguments  a  la  Privalov’s 
theorem  can  be  used  in  the  proof.  Alternatively,  one  could  proceed  in  the  usual  way  by 
Jensen’s  theorem  after  setting  up  the  following  mapping.  Choose  C  C  K,  \C\  >  0,  such  that 
|/(2)|  <  A  on  (J  T7,  7  G  F,  where  Ty  C  3?x  [0,  6)  is  a  triangle  and  Ty\ {7}  C  5?x  (0,  6).  For  an 
appropriate  set  B  D{JTy  choose  a  bi-absolutely  continuous  map  from  B  onto  the  unit  disc; 
and  then  apply  the  Jensen’s  theorem  argument  found  in  the  books  of  Koosis  (pp.  76-77), 
Levinson  (p.74),  or  Porcelli  (p.60). 
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Lemma  3.2.  Suppose  F  €  T2(3?)  \  {0}  and  supp  F  is  contained  in  a  half-line.  Then 
F  —  f  E  // 2 ( 5? )  \  {0}  and,  in  particular,  f  is  an  analytic  function  with  3?  as  a  boundary. 

Lemma  3.2  is  the  easy  direction  of  the  Paley- Wiener  theorem,  and  also  ties  in  with  the 
work  of  F.  and  M.  Riesz,  Levinson,  and  Beurling  and  Malliavin.  In  light  of  Lemma  3.2  and 
the  following  application,  we  note  that  we  are  only  using  Lemma  3.1  in  the  case  /  €  //2(9?). 

Using  Lemmas  3.1  and  3.2,  our  complicated  proof  of  (3.1)  proceeds  as  follows.  We  first 
write 

v(m,  n)  =  ||(t  -  <o)^m,n(Oll2ll(7  -  7o)^m,n(7)ll2 

=  J \(u  —  to  +  na)<f>(u)\2  du  j  [(A  —  70  +  m6)^(A)|2  dX. 

For  case  1,  suppose  supp  <f>  is  not  contained  in  a  half-line  so  that  for  all  n  >  0 


v(m,n)  >  ( na  )2  J  \<j>(u)\2  du  J \(\  -  ~/o  +  mb)d>(\)\2  d\. 


Thus, 


Vm,  lim  v(m,  n)  =  00. 
n— *00  v  ' 

For  case  2,  suppose  supp  (f>  is  contained  in  a  half-line.  Then  for  all  m  >  0, 

v(m,n)  >  (mb)2  f  \(u  —  to  +  na)(f>(u)\2  du  j  |^(A)|2dA. 

J  7 70 

By  Lemmas  3.1  and  3.2,  |^(A)|2  d\  >  0  and  so 

Vn,  lim  vim,  n)  =  00. 
m— *00  v  ’  ' 


3.3  Remark. 

(2.1),  (3.1),  and  Proposition  2.1  can  be  compared  with  Bourgain’s  weak  uncertainty  in¬ 
equality  [Bou88],  cf.  [Ben90,  section  4].  Bourgain’s  orthonormal  basis  was  constructed  by 
perturbations  of  Heisenberg  wavelets. 

4  The  Balian  conjecture. 

The  Balian-Low  theorem,  mentioned  in  Section  1,  was  first  stated  by  Balian  [Bal81]  in 
1981.  At  the  same  time  he  posed  the  following  conjecture :  there  are  no  orthonormal  bases 
{#*}  Q  L2(3J)  with  the  property 


sup  A  if.  <  00  and  sup  A  7^  <  00, 
k  k 


where 


A t\  =  (0k(t),  t29k(t))  -  ( 0k(t ),  tdk{t))  and 

A 7 k  =  (h(l),l2h(l))  ~  (hh),30k(l)), 

cf.  Proposition  2.1. 

The  proof  of  (2.1)  shows  the  validity  of  the  conjecture  for  affine  wavelets,  and  the  Balian- 
Low  theorem  confirms  the  conjecture  for  Heisenberg  wavelets.  Bourgain’s  theorem  casts  some 
doubt  on  the  conjecture,  cf.  [D  J  J] . 
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Section  III 

A  Wiener-Wintner  Theorem 


A  MULTIDIMENSIONAL  WIENER-WINTNER  THEOREM  AND  SPECTRUM  ESTIMATION 


John  J.  Benedetto 
Prometheus  Inc. 


Abstract 


Sufficient  conditions  are  given  for  a  bounded  positive  measure  p  on 

IRd  to  be  the  power  spectrum  of  a  function  <p.  Applications  to  spectrum 

estimation  are  made  for  the  cases  in  which  a  signal  <p  is  known  or  its 

autocorrelation  P.  is  known.  In  the  first  case,  it  is  shown  that 
<P 


If  Cr) 


'dp^(y)  = 


lim 

R-+a> 


I  B(R) 


|f«<p(t)|  dt, 


BCR) 


where  P  =  p  ,  B(R)  is  the  d-dimensional  ball  of  radius  R,  and  f  ranges 
<P  <P 

through  a  prescribed  function  space.  In  the  second  case,  an  example,  which 
is  a  variant  of  the  Tomas-Stein  restriction  theorem,  is 


Vf  €  L1(IRd)nLP(lRd), 


(| 


If  (e) 


"dpd-i(0) 


1/2 


2  Md- 


^2]  [llfll 


■r 


iifii 


»]• 


'd-l 


where  1  <  p  <  2d/(d+l)  and  th?  power  spectrum  p  is  the  compactly 
supported  restriction  of  surface  measure  to  the  unit  sphere 
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1.  Introduct ion 


Wiener’s  generalized  harmonic  analysis  is  a  multi-faceted  subject.  It 
was  conceived  to  deal  with  noise,  hidden  periodicities,  and  turbulence;  it 
was  developed  in  terms  of  Tauberian  theorems  and  sophisticated  spectral 
analysis;  and  it  continues  to  play  a  role  in  prediction  theory,  ergodic 
theory,  stochastic  processes,  and,  of  course,  harmonic  analysis  [B;W).  An 
important  fact  from  the  general  theory  is  a  theorem  due  to  Wiener  and  Wintner 
[ WW, 241-242] :  each  bounded  positive  measure  p  on  the  real-line  is  the 
Fourier  transform  of  a  function  P  which  can  be  written  in  the  form, 


i.e.,  each  p  is  the  power  spectrum  of  some  signal  <p  having  autocorrela¬ 
tion  P.  In  Section  3  we  extend  this  result  to  functions  <p  on  d-dimensional 
euclidean  space  Rd;  and  we  use  the  theorem  as  background  for  proving  a 
d-dimensional  type  of  Wiener-Plancherel  formula  in  Section  5  and  a  weighted 
Fourier  transform  norm  inequality  and  restriction  theorem  in  Section  6.  The 
results  in  Sections  5  and  6  are  viewed  in  terms  of  estimating  the  power 
spectrum,  i.e.,  "spectrum  estimation". 

Section  2  is  devoted  to  notation.  Theorem  3.3  (in  Section  3)  is  the 
d-dimensional  version  of  Wiener  and  Wintner’ s  theorem.  It  should  be  pointed 
out  that  their  proof  is  cryptic,  and  our  proof  on  Rd  is  inspired  by  the 
intricate  argument  of  Bertrandias  for  R  [Be],  cf . ,  [Ba,  Chapitre  2]. 

Section  4  provides  a  geometrical  application  of  Theorem  3.3  for  u.  ..  the 

d_  1 

restriction  of  surface  measure  <r,  ,  to  the  unit  sphere  Z,  of 

d-1  d-1 

d-dimensional  euclidean  space  Rd(=!Rd). 


In  Theorem  5 . 2  (of  Section  5)  we  derive  a  formula  to  express 
J|  f (?)  I ^d  p^(y )  in  terms  of  arithmetic  means  involving  f  and  a  gi 


III2 


is  a  bounded  positive 


function  <p.  (Integration  is  over  R  and  p 


measure  on 


which  is  the  Fourier  transform  of  a  d-dimensional  version  of 


the  function  P  in  (1.1);  the  function  f  denotes  the  Fourier  transform 


*  I 


f(t)e  2Tt^  ^dt,  where  integration  is  over  Rd  and  y  €  Rd. )  Using 


Theorem  5. 2  with  functions  f  having  the  appropriate  shape,  we  can  estimate 
the  mass  of  p^  on  prescribed  regions  of  Rd  in  terms  of  the  known 
functions  f  and  <p.  This  is  a  form  of  spectrum  estimation  and  further 
remarks  are  made  in  Section  5. 

The  results  in  Section  6  are  elementary  but  of  some  interest  because  of 
recent  work  on  weighted  norm  inequalities  and  restriction  theorems.  The 

2  p 

(L  ,1/)  restriction  theorem  of  Tomas  and  Stein  asserts  that 


(1.2)  Vf  €  L1(Rd)nLP(Rd) 


ff  -  2  11/2 

.  I  | f  ( 0 )  |  d<rd_ J ( e )  <  C(p) Ilf  1 


for  any  fixed  1  <  p  <  2(d+l  )/(d+3) ,  e.g.  ,  [  T) .  LP(Rd)  and  ||f||  are  the 

P 

usual  Lebesgue  space  and  norm.  Discussions  of  the  role  of  curvature  and 
extensions  of  (1.2)  to  other  smooth  sub-manifolds  are  found  in  [CD,  103-109; 
S,  325-329],  Weighted  LP  versions  of  (1.2)  are  found  in  [ BH ] .  Proposi¬ 
tion  6. 1  is  a  p-weighted  Fourier  transform  norm  inequality  with  explicit  norm 
constant  in  terms  of  the  inverse  Fourier  transform  P  of  p;  and  Proposition 
6. 2  is  a  corollary  which  can  be  viewed  as  a  restriction  result  since  we  take 
p  =  p^  The  inequality  (1.2)  is  deep  and  Propos i t ion  6.2  is  elementary. 

In  (1.2)  the  constant  c(p)  is  not  explicit,  p  =  2(d+l)/(d+3)  is  largest 
possible,  and  the  right  hand  norm  is  }j F || ^ .  In  Proposition  6.  2  the  constant 
is  explicit  and  the  values  of  p  extend  beyond  2(d+l )/(d+3) ,  but  the  right 
hand  norm  is  l|f||^  +  II f"  1! ^ .  As  far  as  spectrum  estimation  is  concerned, 

Proposition  6. 2  can  be  used  to  give  an  upper  bound  of  the  mass  of  p  ,  on 

d- 1 

prescribed  regions  of  Rd  in  terms  of  the  known  functions  f  and  P.  This, 
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too,  is  a  form  of  spectrum  estimation  and  further  remarks  are  made  in 
Section  6.  Naturally,  in  light  of  Theorem  3. 3  and  the  generality  of 
Proposition  6.1.  our  discussion  of  pd-1  extends  to  many  other  measures. 
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2.  Notation 

Besides  the  L^-spaces,  1  <  p  <  co,  mentioned  in  Section  1.  we  shall 
deal  with  a  number  of  other  spaces.  To  introduce  them  we  let  IRd(  =  IRd)  be 
the  dual  group  of  d-dimensional  euclidean  space  Rd  and  let  X  be  a  locally 
compact  subspace  of  !Rd.  Then  CMX)  is  the  vector  space  of  complex-valued 
continuous  functions  f  :  X  — *  C  having  compact  support  supp  f  £  X.  A 


measure  p  on  X  is  a  linear  functional  defined  on  CMX)  satisfying 

lim  <p,f  .>  =  0  for  every  sequence  {f  .}  £  C  (X)  having  the  properties  that 
.  0  J  c 

J-*°° 

lim  || f  ,||  =  0  and  supp  f  .  £  K,  where  K  £  X  is  a  compact  set  independent 

J  00  J 

J-»00 


of  j  and  jj  •  ■  ■  I!  is  the  usual  sup  norm,  e.g.,  [Bo].  M ( X )  is  the  space  of 

CO 

measures  on  X  and  M+(X)  =  {p  e  M(X)  :  <p,f>  >  0  for  all  non-negative 
f  e  ( X) >  is  the  space  of  positive  measures  on  X.  Similarly,  ( X )  is 
the  subspace  of  M(X)  having  bounded  variation,  i.e.,  the  above  mentioned 
convergence  criterion  on  C  (X)  is  replaced  by  (C  (X),  ||  •  ■  •  ||  );  and 

C  C  oo 

M,  (X)  consists  of  the  positive  elements  of  M,  (X).  The  support  of 
D+  D 

p  €  M(X)  is  denoted  by  supp  p.  We  write  <p,f>  =  ffaOdpty)  and  in  the 

JX 
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A  ^  p 

case  X  =  IR  we  write  <4,f>  =  I  f(  j- )dji(y ) .  If  )i  e  ( X )  then  4  is  well 


defined  on  C,  (X)  =  {f  :  f  is  continuous  and  bou  ided  on  X}. 
b 

For  p  e  (0,oo),  Lf  (IRd)  is  the  set  of  functions  f  :  IRd 

loc 


C  for 


which  |fr  is  locally  integrable  with  respect  to  Lebesgue  measure.  If 

4  e  M+(!Rd)  then  L^(IRd)  designates  the  set  of  Borel  measurable  functions 

f  defined  4  a.e.  on  (Rd  for  which  ||f||  =  (  f|f  (7)  1  Pdpi  t  )  )1/p  <  00. 

P>  P  J 

There  is  the  standard  adjustment  for  p  =  00.  Also,  we  write  p'  =  p/(p-l). 

The  Fourier  transform  defined  in  Section  1  extends  from  L^([Rd)  to 
M^(IRd)  (and  beyond),  arid  4  is  the  inverse  Fourier  transform  of 
4  €  Mb(IRd). 


The 

characteristic 

function 

of  the 

set 

S  £  lRd 

is  %s 

and  the 

Lebesgue 

measure  of  S 

is  denoted  by 

ISi. 

We  set 

B(R)  = 

{  x  6  IRd 

Ix|  <  R> 

so  that  1  BCR)  |  =  a) .  , 

a- 1 

Rd/d 

where 

M  = 

d-1 

2nd/2/(dr(d/2) ) 

surface  area  of  Finally,  we  mention  that  we  shall  frequently  deal 

with  sets  of  the  form  B(R)\B(r),  and  that  the  boundaries  of  such  sets  are 
unimportant  for  our  results.  As  such,  on  this  point,  we  shall  not  usually  be 
concerned  about  boundaries  of  such  sets  in  our  calculations. 


3.  The  Miener-Wintner  Theorem  on  1R 

Given  4  €  Mb+(R  )  and  let  5^  be  the  Dirac  measure  supported  by  {u} 
It  is  well-known  that  there  is  a  sequence  {4^}  £  M^fsupp  4)  of  positive 
discrete  measures, 

N 

n 

4  =  /  a  .  6 

n  L-,  J .  n  u 

j  =  l 


,  a .  >0, 

w .  j,  n 

J.  n 


such  that  (a)  :  j  =  1,**-,N  }  £  supp  4  for  each  n, 

j ,  n  n 
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and  lim  p  =  p  in  the  (vague)  topology  cr( ( IR  ),  ( R  )),  e.g.,  [Bo, 

n-*oo  "d  'd 

Chapitre  III,  Section  2,  no.  4].  Actually,  (3.1)  and  the  criM^OR  ),  )) 

convergence  allow  us  to  conclude  that  lim  p  =  p  in  the  "Levy"  topology 

-d  ~d 

crd^dR  ),  Cb((R  )),  e.g.,  [  Ma,  91-93], 

For  a  given  p  e  Mb+([Rcl)  and  sequence  {p^}  -  Mb+(supp  p)  as  above,  we 
define 


v  1/2 

n  J.n 


2ni  t  ■  ti> . 

J>n 


so  that 


„  .  V  1/2  .  ..  1/2 

o  <  >  a .  <  N  sup  a . 

n  “  £—>  J* n  n  !<  c<n  J-n 

J=1  n 


Lemma  3.1.  For  each  n. 


(3.2) 


TEXT  L„, 


R-»oo  |B(R)I  Jb(R) 


uniformly  on 


Proof.  The  left  hand  side  of  (3.2)  is  p  (t)  +  lim  c(R,t,n)  (if  this  limit 

n  R-*co 


exists),  where 


,n  .  ,  1  V  1/2  1/2 

c(R,  t.n)  =  .  d-75t  r  >  a.  a, 

I  B(R)  |  L_  j ,  n  k,n 


2tu  t  •  cj 


27iix'  (u .  -u,  ) 

e  J’n  k’n  dx 


O  -t  _fd~2.  rR  d 

271  V  1/2  1/2  j.n,  ,  2  }  2  .  . 

in,  n  rr  /  a .  a,  e  w .  -w.  r  J,  A2nr  g>  .  -w,  )di 

I  B(  R)  I  /_  j ,  n  k,  n  j.n  k,n  d-2  j.n  k,n 


and  J,  „  is  the  Bessel  function  of  order  (d-2)/2.  Since  |J,  „(s)|  <  1 
d-2  d-^ 

~2~  ~2~ 

- 1  /2 

(at  the  origin)  and  J  (s)  -  0(s  ),  s — *»,  we  see  that  lim  c(R,t,n)  =  0 

R-.00 

2 


independently  of  t. 


i.  e.  d. 


We  shall  also  need  the  following  easy  fact. 


co  d 

Lemma  3.2.  Given  a  sequence  {R  }  increasing  to  infinity  and  </>  e  L  OR  ). 

n  i  oc 


Fix  t  e  IR  .  Then  for  each  fixed  p, 


I  B(R) 


<p[  t+x)<p(x)dx  =  0. 


BCR  ) 
P 


We  use  Lemma  3. 1  to  define  a  specific  sequence  { R^}  and  a  specific 


function  <p  in  the  following  way.  From  the  uniform  convergence  we  know  that 


Vn  >  1,  3A  >  A  .  such  that  Vt  e  R  and  VR  >  A  , 
n  n- 1  n 


(3.3) 


1 

I  B(R) | 


ip  (t  +x)p  (x)dx  -  P  (t)  <  - r-  . 

n  n  n  2°  1 


We  set  R  =  ( A  +1 ) ( A„+2) • • • (A  +n)  so  that  each  R  >  n!  and  the  sequences 
n  1  2  n  n 

(R  >,  <R  /R  }  =  (A  ,  +n+l}  ,  and  {R  -R  }  =  {  ( A  +1 )  •  •  •  ( A  +n)(A  Al+n)> 

n  n+1  n  n+1  n+1  n  1  n  n+1 

increase  to  infinity.  For  this  sequence  {R^}  and  for  {(/M  defined  above 
we  define  <p  on  by  setting  ^(t)  =  v^t)  for  R^<  1 1 1  <  R  ,  n  >  1, 

and  letting  <p(t)  =  0  for  |t|  <  R  Clearly,  <p  s  L™Qc(Rd)  and  va-lues 
of  <p(t)  for  |t|  =  r^  are  not  important. 


Using  this  notation,  setup,  and  the  lemmas,  we  can  state  the  Wiener- 
Wintner  construction  on  R^. 
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Theorem  3.3.  Given  p  e  M^dR  )  with  corresponding  functions  and 


(p  €  L?  (!Rd).  Assume 
loc 


n 

(3.4) 

lim  i-  /  max  ( ||  <p  .  1 1I  ,  ll<p  ||  )  max  ( ||  II  , 

n->co  n!  ^  J-l  co  J  CO  J  « 

j=2 

=  0. 

Then,  for 

each  t  e  IR^, 

(3.5) 

Plm  1 B( R) | 
R-*oo 

<p(  t+x)^>(  x)dx  =  p 

JB(R) 

(t). 

Proof,  a. 

Given  t  e  and 

c  >  0.  Because  of  our 

setup  we 

can  choose 

P2  =  Pj(t) 

for  which  R^+j  ~ 

|t|  >  Rk  +  |t|  for  all 

k  -  Pj 

and  we  can 

choose  p2  =  p2(c,t)  >  p^  such  that 


(3.6) 


Vk  >  p2,  ||ik(t)  -  p  ( t )  |  <  c. 


We  shall  find  p  =  p(e,t)  >  p2  such  that  for  all  sufficiently  large  R, 


(3.7) 


y(t+x)^>(x)dx  -  p  (t)  <  14c; 


B(R)\B(R  ) 
P 


and,  hence,  (3.5)  is  a  consequence  of  Lemma  3. 2  and  a  "lim"  argument, 
b.  For  each  p  >  p„  and  each  R  >  R  ,  write 

2  p 


(3.8) 


1 

I B  ( R ) 


B(R)\B(R  ) 
P 


y>(t+x)ip(x)dx  =  ^(b(k,R)  +  c(k,R)) 
k=P 


<p(t+x)y>(x)dx, 


B(R)\B(R  ) 
n 


where 


III- 


b(k.R)  = 


I  B(  R)  | 


<p(  t+x)<p(  x)dx, 


c(k,R)  = 


I  BCR)  | 


<p(t+x)</>(x)dx, 


BUW'BnWltl) 


and  n  =  n(R)  is  the  largest  integer  n  for  which  R^  <  R. 
In  part  c  we  shall  verify  that  for  p  > 


(3.9) 


n-1 

1 im  c(k, R)  =  0. 
R-»co 

k=p 


Parts  d  -  h  are  devoted  to  showing  that  there  are  p^  =  p^Cc.t)  >  p and 
R(c,t)  >  p^  such  that,  for  R(c,t)  >  p  >  p^, 


VR  >  R(c 


,t).  £ 


b(k, R) 


(3. 10) 


M  (t)  I  BCR)  |  X  UB(Rk+l 


1 1 1 )  |  -  IBCiyi)  +  r(p,R,t), 


where  |r(p,R,  t)|  <7 c  +  -■  -  .  Parts  i.  -  k  contain  the  proof  that 


(3. 11) 


VR  >  R(c,t), 


I B( R) | 


_  „  I B  ( R  )  | 

<p(t+x)<p(x)dx  =  m  ( t )  ( 1  -  i B ( )  +  s ( R,  t ) , 


B(R)\B(R  ) 
n 


where  |s(R,t)|  <4 c  +  — .  (3.7)  is  obtained  in  part  X  by  combining  (3.9), 

2n 

(3. 10),  and  (3.11). 


c.  If  x  e 


B(Rk+1)\B(Rk+1  -  | t | )  then  <  |x+t|  <  RR+2.  Thus, 


setting  cR  =  II v>k II ^  max(  Ilvk+1 II J .  we  have 
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n-1 


n-1 


Z  =(k-R>  s  ~H  Z  ck[<*i  -  (\.rltl)d 

i  K  ,  L  J 


k=p 


k=p 

n-1 


,  >  y  r,  -  n-iii-)11] 

Rd  A  k  ^*n  >w  J 


n  k=p 
n-1 


S  r1  Z  Ck  "ktl  Z(j](\Z 


n  k=p  j=l 

n-1 


‘[T)l 

k=P 


d  "  ,  s  iiHiii! 

k  k+1  n! 


n-1 

z 

k=l 


V 


(3.9)  follows  by  comparing  the  right  hand  side  of  this  estimate  with  our 
hypothesis  (3.4). 

d.  In  order  to  estimate  b(k,R)  we  define  the  sets 
A^t)  =  (x  6  BfR^-ltlKBdy  :  t  +  x  e  B(Rk+1  )\B(R^)>  and  A^ft)  = 
{x  e  B(Rk+1“  1 1 1  )\B(R^)  :  t  +  x  €  BOysBU^ )} .  Then,  for  k  >  pr 
-|t|)\B(R^)  =  A^( t )  \j  A^_  ^ ( t ) ,  a  disjoint  union;  and  we  have 

<Pk(t+x)#>k(x)dx 


b(k,  R)  = 


1 


I  B(R)  | 


(3. 12) 


JAk(t)uB(Rk) 


1 

I  B( R) | 


¥>k(t+x)vk(x)dx  + 


1 

|  B  ( R )  | 


<Pk_i ( t+x)^p  (x)dx. 


B(V 


A 


k-1 


(t) 


We  shall  estimate  these  three  integrals.  The  first  two  will  involve  the 

factor  p  (t)  and  the  latter  becomes  small  by  the  size  of  Ak_^(t). 

In  the  process  of  making  these  estimates  we  need  the  bound, 

n 


(3.  13) 


1 


n 


z 


K 


<  2. 
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In  fact,  the  left  hand  side  of  (3.  13)  is  bounded  by 


,  1  1  1 
1  +  -  +  — , - T-T-  +  *  '  ‘  + 


n  n(n-l) 

and  it  is  easy  to  see  that  for  each  k  >  3 

k+1  k 


n! 


(k+1 


ITT  I  J*  <  ET  Z  s 


J=0 


J=0 


5 

3’ 


hence,  (3.13)  is  valid. 

e.  Our  initial  step  in  estimating  the  first  integral  on  the  right  hand 
side  of  (3.12)  is  to  prove  that  there  is  p^  =  p^(c,t)  >  p^  such  that  for 


all  k  >  p3« 
(3. 14) 


1 

* 

^k(t+x)<pk(x)dx  - 

<Pk(t+x)ipk(x)dx 

!B(R  +  -!t 1 ) ! 

B(Rk+rltl) 

Ak(t)uB(Rk) 

The  difference  of  the  two  integrals  is  the  integral  over  A^_^(t),  and  if 
X  6  Ak  i^  then  <  |x|  <  +  |t|.  Thus,  setting  d^  =  the  left 

hand  side  of  (3. 14)  is  bounded  by 


(3. 15) 


[(Vm)d  -  r£]  * 


d  d 


0W|t|) 


k  d(Vltl)d'1’ 


where  this  last  inequality  is  a  consequence  of  the  mean  value  theorem.  By 
(3.4),  the  right  hand  side  of  (3.15)  tends  to  0  as  k  — >  co;  and  (3.14)  is 
obtained. 

By  means  of  (3.3),  (3.6),  (3.14),  and  the  triangle  inequality  we  have 


(3. 16) 


where 


I  BCRk+-l“  1 1 1  )  I 


k  >  p3-  Clearly, 


<p^( t+x)<p^(x)dx  -  p  (t) 

JA  (t)uB(R^ 

we  can  write  (3. 16)  in  the  form, 


<2  c  + 


,k+T 


IH-11 


(3. 17) 


IB(Rk+1-lt|)| 


A  (t)uB(R^) 


t+x)^k(x)dx  =  p  (t)  +  /3^_  ( t )  ( 2c  +  -g-j-j-), 


where  |/3  (t)|  <  1  and  k  >  p^. 

f.  To  estimate  the  second  integral  on  the  right  hand  side  of  (3. 12)  we 
use  (3.3),  (3.6),  and  the  triangle  inequality  to  write 


(3. 18) 


B(V 


B(V 


V>k(t+x)</>k(x)dx  -  p  (t)  <  c  +  j-, 


for  k  >  p^.  Analogous  to  part  e,  (3.18)  yields 


(3. 19) 


iBtiyi 


B(V 


^k(t+x)<pk(x)dx  =  p  (t)  +  *k(t)(c  + 


where  lyk(t)|  <  1  and  k  >  p^- 


g.  Since  Ak  l(t)  S  BCR^ 1 1 1  ) \B ( ) ,  the  third  integral  on  the  right 
hand  side  of  (3. 12)  is  bounded  by 


(3.20) 


IBdyi 


l^>k_1(t+x)<pk(x)  Idx  < 


Ak_l(t) 


K-lUVco  r  d 

-k-1-“--k  00  |(Rk+|t|)d 
^n 


-4 


h.  By  parts  e,  f,  and  g,  we  can  write 


I  b(k’R)  '  TbTrTT  Z  IB(Rlc,1-ltm(MV(U  *  Pk(t)(2c  *  -i^)) 

k=p  k=p 


(3.21) 


TbTrTT  Z  *  rk(t)(c  +  -j^) 


n-1  p 

TbTrTT  X 

l,— ^  * 


Vk-1  (t+x)<pk(x)dx- 


k=p  Ak-1 ( t ) 
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To  estimate  the  last  term  in  (3.21)  we  use  (3.20)  and  compute 


n-1 

Jd  Z  "viUV.  [(vtl)d  - Rd|] 

n  k=p 


d|t| 


n-1 

I 

k=P 


:r1c+  i  t  i  ) 


d-l 


<  d  1 1 1 2 
R 

n 


d-1  V-1  Rv  H-1 

,  n 

k=P 


By  (3.4),  this  term  is  less  than  e  for  all  large  R.  Consequently,  (3.21) 
allows  us  to  obtain  (3.10).  It  is  at  this  point  that  we  invoke  (3.13). 


i.  We  now  estimate  the  last  term  in  (3.8).  Analogous  to  part  d  we 

define  the  sets  AR  ^(t)  =  {x  €  B(R)\B(Rn)  :  t  +  x  e  B(Rk+1  HBlR^)}  for 

k  =  n  -1,  n,  n  +  1.  Then,  for  n  >  B(R)\B(Rn)  = 

A  (t)uA  (t)uA_.  .  (t),  a  disjoint  union.  (Actually,  we  have 

r\,n“i  n,  n  a,  n+i 

n  >  p^. )  Therefore,  the  last  term  of  (3.8)  is 

<p  (t+x)^  (x)dx, 
n  n 


IbTr) 


¥>(t+x)^(x)dx  = 


1 


I  B  ( R )  I 


JB(R)\B(R  ) 
n 


jAd  (t)uB(R  ) 
R,  n  n 


(3.22) 


|B(R  )| 

» 

1 

n 

<p  (t+x)p  (x)dx 
n  n 

1 B(R) | 

1  B( R  )| 

U 

B(R  )  J 

I  B(  R) 


■* 

<p  ,  (t+x)m  (x)dx  + 
n-1  n 

<p  ,  (t+x)</>  (x)dx 
n+1  n 

UAn _ _  (t)  J 

An  _ , i ( t )  J 

j.  The  last  term  of  (3.22)  is  estimated  as  in  part  g.  The  inclusions 

Ad  ,(t)  £  B(R  + | t I )\B(R  )  and  AD  ,(t)  £  B(R)\B(R  -|t|)  are  valid, 

R,n-1  n  n  R, n+1  n+1 

but  it  should  be  observed  that  A  ( t )  can  be  the  empty  set,  e.g. ,  if 

K,  n+ 1 


III-  13 


-  1 1 1  .  In  any  case,  this  last  term  of  (3.22)  is  dominated  by 


R  <  R 


n+ 1 


:3'23)  3  {'V.U'n'.  [<Vltl)d  -  Rn>]  *  -  (V.-|tl)d]} 


where  “lRd  -  (R  -|t|)^]“  is  defined  to  be  0  if  R  <  R  .  -  |t|.  Using 
n  n+ 1 

the  mean  value  theorem  again,  (3.23)  is  bounded  by 


■R_  + 1 1  |^d-l  d(R-(R  -It  I  ) ) 

R  "T’n-l"co"lf'n"ooi  R  |  +  R  11  V-lUV® 

n  ^  n 


,  i .  i  /-R  +  1 1 1  < 

41  'Vl'J'n'.pH - ] 

n  >*  n  * 


Our  hypothesis  (3.4)  allows  us  to  conclude  that  this  quantity  is  less  than  c 

for  all  large  R.  (For  the  second  term,  recall  that  R  <  R  , .) 

n+1 

k.  Our  estimation  of  the  first  two  integrals  on  the  right  hand  side  of 
(3.22)  follows  steps  e  and  f.  For  example,  analogous  to  (3.14),  we  have 


1 

IBCR) | 


¥>n(t+x)y>n(x)dx  - 


B(R) 


r 

_(t] 


<p  (t +x)(»  (x)dx 
n  n 


,r  ,.)uB(R) 

R,  n  n 


IBCR)  |  ^n11.*  [|AR,n-lU)l  +  |AR,n+l(t)  [}] 


and  the  right  hand  side  is  less  than  c  for  all  large  R  by  the  calculation 
used  in  step  j.  Proceeding  in  this  way,  as  in  parts  e  and  f,  we  obtain 


I B  ( R ) 


?>n(t+x)<pn(x)dx  =  p  (t)  + 


Vnlt)(2c 


,n+l ' 


Ad  (t)uB(R  ) 
R,  n  n 


and 


where 
(3. 22) , 


<p  (t+x)ip  (x)dx  =  p  (t)  +  y  (t.)(e  +  — 5_), 
n  n  R,n  n+1 

B(R  ) 
n 

1/3  (t)|,  |  y  (t)|  <  1.  Combining  these  facts  with  part  j  and 

K,  n  K, n 

we  have  verified  (3.11). 


1 

IB(Rn) | 
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1.  For  p  >  large  enough  and  for  all  large  R,  n(R)  >  p,  the  right 

hand  side  of  (3.8)  is 


(3.24)  p'(t)  + 


p  ( t ) 


I  B(R) 


n-1 


Z  (,B(Rk*l“ltni  -  IB,V']  ■  l8<Rn 


)l 


L  k=p 


plus  an  error  term  bounded  by  13e.  This  is  a  consequence  of  (3.9),  (3.10), 
and  (3.11).  We  rewrite  the  second  term  of  (3.24)  as 


pv(t) 

I B(R) | 


[( 


I  B( R  )|  - 
n 


+  •  •  •  + 


|B(Rn-|t|)|J  +  j^|B(Rn  l)|  -  |B(Rn  l-|t|)|j 
[lB(Rp+1)l  -  IB(RpM-ltl)l]  .  IB(Rp)l], 


which,  in  absolute  value,  is  bounded  by 


Ii£-friEK~  (vltl,d]  tRp] 

K  k=p 


<  lMV(t)| 


d  1 1 1 
R 


rR  -vd 


and  which,  in  turn,  is  less  than  c  for  all  large  R.  Thus  (3.7)  is  veri¬ 
fied,  and,  by  the  observation  in  part  a,  the  theorem  is  proved. 

q.  e.  d. 


2 

Remark  3.4.  The  hypothesis  (3.4)  can  be  weakened.  For  instance,  llv^ll^  -  CR^ 
is  sufficient.  There  is  also  left  unanswered  the  problems  of  constructing  <p 

A  ^ 

for  arbitrary  p  e  M^+(!R  )  and  for  characterizing  those  p  for  which 
<p  €  L^dR^).  The  former  problem  is  answered  positively  for  all  p  in  case 
d  =  1;  the  latter  is  not  yet  solved  in  case  d  =  1. 
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4.  Example 


A  ^ 

Let  =  Ml?  I  -  D  6  M^+(!R  )  be  the  measure  corresponding  to  a 

mass  distributed  homogeneously  over  E,  , ;  the  total  mass  is  taken  as  o>  . 

a-i  d-l 

so  that  the  compactly  supported  measure  pd  is  the  restriction  of  surface 

measure  cr^  to  E^_ ^ .  The  transform  pd_^(t)  is  easily  computed  in  terms 

of  the  approximants  —  — -)  to  p,  , ,  where  *(r)  =  x,  »(r)  for 

t  t  'a-l  [-1/2,  1/2) 

r  €  (R.  In  fact,  we  have 


p,  ,  (t)  =  lim 
x-*0  J 


1  ,  lrl-1,  2nit-y 

-  Xi~~ — )  e  d? 


(4.  1) 


=  lim  - 
x->0  T 


d- 


■(I. 


2irit*p0  ,  ,  -  , 

e  H  d<r  ,  (0)  | dp 


'd-l 


-(**) 

=  2ir|t|  2  Jd_2(2ir|t| ), 


since 


27rit-p0  , 
e  K  dcrdl(0) 


=  2tt(  1 1  |p) 


.(**) 

1  2  ' 


'd-l 


Jd_2(27i|t|p). 

~2~ 


2  cl 

Theorem  4.1.  Given  d  >  2  and  the  measure  u.  ..  There  is  w  e  L.  (IR  ) 
-  d-l  r  loc 

such  that 


Vt  e  (R 


1  im 
R->co 


1 

I B  ( R )  | 


<p(t+x)<p(x)dx  =  p  (t) . 

a- 1 

JB(R) 


Proof.  The  result  will  follow  from  Theorem  3.3  once  we  construct  (u  >  and 
-  -  n 

{ <Pn }■  as  in  Section  3  where  satisfies  (3.4). 
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To  this  end  we  first  cover  E,  ,  in  the  following  way.  For  each  n  we 

d- 1 

choose  about  n^  *  balls  B.  ,  j  =  1,--*,N  <N,  n^\  each  having 

j,n  °  n  d 

diameter  bounded  by  D,/n  and  whose  union  covers  Z,  For  d  =  2,  this  is 

d  d-1 

doable  since  the  circle  can  be  divided  into  n  arcs  of  equal  length  2n/n, 

and,  hence,  n  balls  of  diameter  D^/n  can  be  chosen.  For  d  =  3,  we  first 

consider  the  equator  and  the  polar  great  circle.  We  choose  n  equi-spaced 

points  on  both  circles.  For  convenience  let  n  be  even  so  that  the  polar 

points  are  marked  on  the  polar  great  circle.  Then  spin  the  polar  great 

circle  stopping  at  each  marked  point  on  the  equator.  At  each  stop  we  mark 

points  on  Z^  determined  by  the  rotated  points  of  the  polar  great  circle. 

2 

In  this  way  we  obtain  Nn  <  N^n  points  on  E^.  For  any  marked  point  p  on 

the  equator  we  choose  a  ball  B(p,r)  where  the  radius  r  is  large  enough  to 

include  two  marked  points  of  the  equator  on  either  side  of  p.  Balls  of  this 

size  centered  at  the  other  marked  points  on  Z2  yield  the  desired  set 

{B.  }  (since  B(p,r)  covers  the  largest  area  on  E  between  a  block  of  6 

J  >  C* 

neighboring  points).  For  d  >  3  we  use  this  process  by  considering  d  -  1 

cl—  1 

"great  circles"  Z,  „  to  find  the  n  points  and  isometric  balls  B. 
s  d-2  r  J,n 

centered  at  these  points.  Analytically,  this  can  be  accomplished  by 


imbedding  Z^  ^  in  the  unit  cube,  covering  the  surface  of  the  unit  cube  by 
cl  (J— •  i 

2  n  balls  of  diameter  D^/n  (equidistributed  on  each  face  of  the  cube), 
and  linearly  sliding  each  ball  to  the  origin  where  it  is  designated  Bj  n 
when  its  center  intersects  Z^_^.  Clearly  (geometrically),  these  B  n 


cover  Z 


d-r 


We  can  now  construct  (u  }.  For  each  fixed  n  we  set  B.  (Z.  . )  = 

n  j.n  d-1 

B.  nZ,  J  =  ;  we  then  obtain  a  partition  of  E,  ,  by  defining 

j,n  d-1  n  d-1 

°,  _  =  B  (Z  )  and 
1 , n  1 , n  d-1 
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.  =  B.  (I.  .)  n  fV  K  Jl  •  j  =  2,  •  •  •  ,  N  . 

j.n  j.n  d-1  1  1  k.n  d-1  J  n 


There  is  K^,  independent  of  n  and  j  =  !,•••  N  ,  such  that 


(4.2) 


Md-l(°j,n)  "  d-1' 

^  n 


Besides  the  geometric  proof  of  (4.2)  we  can  proceed  analytically  as  follows. 

Take  g  =  ^  ,  and  define  grf  on  lRd  by  setting  ( 0 )  =  0, 

j.n1  d-r 

gd(p9)  =  g(0)  for  0  <  p  <  1  and  0  e  £d_^,  and  g^fy)  =  0  for  |y|  >  i. 
By  the  definition  of  p 


(4.3) 


g(e)do’d_1(0)  =  g  d(r)dy; 


and  (4.2)  is  obtained  by  calculating  the  right  hand  side  of  (4.3)  using 
spherical  coordinates. 

Now  that  we  have  a  partition  (0.  :  j  =  1,***,N  <  N,nd  }  of  , 

^  j.n  °  n  d  d-1 

satisfying  (4.2)  and  having  diameters  bounded  by  Dd/n,  we  choose 

u.  e  0.  (when  0.  *  <b)  and  define 

j,n  j.n  j.n 


i  =  7(1,  ,(0.  )  8 

n  /  .  d-1  j.n  u. 


<p  ( t )  =  y  p  ( 0  .  ) 

n  ZL  d-1  j.n 


2nit • w  . 

1/2  e  J’n. 


It  remains  to  prove  that  lim  p  =  p  in  the  Levy  topology  (so  that  we  have 

n-*co 

(3.6))  and  that  (3.4)  is  satisfied. 


Ill- 


^  d  ~  d 

The  proof  that  lim  |i  =  p  in  the  Levy  topology,  i.e.,  crlM^UR  ),  C^( [R  )) 
n-»oo 

proceeds  as  follows.  Given  g  e  C^(K^);  since  g  is  uniformly  continuous 


on  Z,  ,  there  is  (r  }  decreasing  to  0  so  that  if  0,0'  e  Z 

d-1  n  o-i 


have 


Euclidean  distance  less  than  D,/n  then  |g(0)  -  g(0')|  <  r  .  Therefore, 

d  n 

for  each  j, 


X0  ( 0  )g(0)do'd_1  (0 )  =  g(w  n) 
J  *  ^ 


*o  (e)dVi(e) 

J.  n 


"d-1 


"d-1 


+  (3  .  r  !  p  ,  ,  I  ( 0  .  ) , 

J, n  n  ^d-1  j, n 


|£.  |  <  1.  Summing  over  j  we  obtain  the  desired  convergence  since 

J  *  ^ 


g(3r)dpn(r)  I  <  wd  l  rn- 


To  verify  (3.4)  we  first  use  (4.2)  to  observe  that 

(£l, 

I*.  *  "a  2  • 


Thus ,  we  have 


n2k  n*‘ 


X  "“'"ViL/’VJ  £  4d  Z  kd  l  £  N 


2k  (n+1 ) 


d  d  n! 


k=2 


k=l 


The  right  hand  side  tends  to  0  and  so  (3.4)  is  satisfied. 


a.  e.  d. 


Remark  4.2.  The  geometrical  construct  ion  in  Theorem  4 .  1  was  made  with  (3.4) 
in  mind.  There  are  other  methods  of  partitioning  Z^_ ^  so  that  the  results 
of  Section  3  can  be  invoked.  For  example,  let  d  =  3  and  consider  the 
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"simplex"  determined  by  the  six  points  (±1,0,0),  t0,±l,0),  and  (0,0,11).  The 
trace  of  the  boundary  of  this  octahedron  on  (with  the  "eraser-end  of  the 

pencil"  at  the  origin)  determines  8  subsurfaces  of  Z^_  ^  having  equal  areas. 
At  the  next  step  we  make  the  barycentric  triangulation  of  each  face  and  pro¬ 
ceed  to  establish  48  subsurfaces  having  equal  areas.  In  this  way  we  produce 
a  partition  growing  exponentially  in  cardinality,  but  compensated  for  by  much 
finer  estimates  than  (4.2).  In  particular,  at  the  klh  step  the  partition 

{0.  }  has  the  properties  that  j  =  l,---,2^6^  1  and  p  (0.  ,  )  =  w  /(2^6k  *). 

J ,  k  c.  J ,  K 

3  k-1  1/2 

Thus  II®,  II  <  2  6  )  and,  hence,  we  have  the  estimate  (for  verifying 

k  oo  2 


(3.4)), 


max  ( ||  <p 


k-1  oo  ’ 


>V'J  - 


n! 


k=2 


k=0 


which  tends  to  0.  For  this  partitioning,  our  verification  of  convergence  in 

the  Levy  topology  also  applies  since  the  diameters  of  the  0.  tend  to  0 

J »  K 

as  k  — )  oo. 


5.  Generalized  harmonic  analysis  and  spectrum  estimation 


2  d 

Definition  5.1.  Given  ®  e  L,  (R  )  and  define 
-  loc 

VR  >  °’  P<p,  R  =  |B(R)T  (lp*B(R)  M<f*B(R)  } 


so  that  P  _  e  L*  (Rd)  £  M(Rd).  Suppose  that  there  is  a  continuous 
ip,  R  loc 

positive  definite  function  P  for  which  lim  P  „  =  P  in  the  cr(M(Rd), 

*  R->®  *’R  V 

C^fR^))  topology.  Then  P^  e  L°°(Rd)  is  the  autocorrelation  of  <p  and 

P  =  u  is  the  power  spectrum  of  <t>. 

<P  <P  - - 
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Given  the  vague  convergence,  lim  P  _  =  P  ,  the  positive  definiteness 

of  P  can  be  easily  verified. 

<P 

Suppose  the  data  characterizing  a  given  signal  < ip  is  known.  In  the 
following  result,  f  can  be  thought  of  as  a  properly  shaped  window  function 
so  that  the  left  side  of  (5.1)  represents  the  power  of  <p  in  the  region 
supp  f .  Formula  (5.1)  provides  a  method  for  computing  this  power  in  terms  of 
the  known  functions  f  and  <p.  In  practice,  then,  numerical  estimates  of 
the  right  hand  side  of  (5.1)  lead  to  a  spectrum  estimation  algorithm. 


2  d 

Theorem  5.2.  Given  <p  e  L,  (DR  )  with  autocorrelation  P  and  power 

loc  <p 

spectrum  p  .  Assume  there  is  an  increasing  function  i(R)  on  (0,co)  for 

2 

which  sup  | ( x )  |  <  i(R)  and  lim  i(R)  /R  =  0.  Then 
I  x  |  <R  R-»oo 


(5.1)  Vf  c  C  (Rd),  |f(r)|2d|i  (y)  =  lim  Tplp-r 

C  *  R-»oo  |B(R) 


|f*<p(t)  Tdt. 


(A  simple  estimate  shows  that  f»<p  €  L2  (Rd)  for  f  e  C  (IRd)  and 

loc  c 

^  LL(Rd>-> 

Proof.  a.  By  Bochner’s  theorem  we  have  e  Mb+(Rd);  and,  by  the  cr(M(Rd; 

C  (Rd))  convergence  lim  P  _  =  P  and  the  Parseval  relation,  we  have 
C  R->oo  *’R  * 


(5.2)  Vf  e  Cc(Rd),  lim  |b?r)|  1  (^B(  R) } ~  ( y)  1 2  1  ^  7 }  1 2dy  =  ITty)  1 2dM  ( y ! 

R-*co  J  ^ 


Parseval  can  be  applied  since  (f»f)P  e  L*(Rd). 

<P 

Using  Plancherel’s  theorem  we  see  that  the  left  hand  side  of  (5.2)  is 


(5.3) 


I T*  ( <p^B(  r)  H  t- )  I  dt; 
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in  particular  this  limit  exists  for  each  f  e  Cc(Rd).  We  complete  the  proof 
once  we  show  that  (5.3)  equals 


(5.4) 


1  im 
R-*oo 


1 

I  B  ( R )  | 


I  f  *<p(  t )  |2dt; 

B(R) 


2  d 

in  particular,  we  prove  that  this  limits  exists,  noting  that  f»</:  e  L7  (Ra). 

loc 

2  d 

b.  Using  the  triangle  inequality  for  L  (R  )  and  the  generalized 
Minkowski  inequality,  we  have  the  estimate 


[l  B(R)  I  1  f*  (^*b(R)  ^ ( ]  1  dt]  [jbT 


I  B( R) | 


.  1/2 


i(f*«»)(t)zB(R)(t)rdt 


i 


I B( R) | 


1/2 


lf*(</>ZB(R)Ht)  ~  (f*<p)  (t)*Dfnn  (t)  |  dt 


B  ( R ) 


1/2 


(5.5) 


IB(R) 


1/2 


f(x) 


-1  1 2  -1  1/2 

(liP*B(R))(t”x)  "  ^B(R)  (t)¥>(t~x)JdX|  dt 


1 


1 B  ( R )  | 


1/2  imf(x)  <raB(R))lt-x)  -  :tB(R)(t,'B<t-x) 


2  -.1/2 
dt  dx 


1 


I B( R) | 


1/2 


If(x)  | 


supp  f 


(^B(R))(t_x) 


*B(R)UMt-x) 


2  -]  1/2 
dt  dx. 


for  each  f  e  C^tR  )■  Consequently,  the  proof  of  the  result  is  complete  once 
we  show  that  the  right  hand  side  of  (5.5)  tends  to  0  as  R  — >  c. 

c.  Let  #>(R,t-x)  =  )  (t-x)  -  *B(R)  ( t  )0(  t-x) .  Then  </r(R,t-x)  is 

0  if  t  €  B(R)  and  t  <=  x  +  B(R)  or  if  t  t  B ( R )  and  t  -  x  e  B(R). 

Also,  y(R,t-x)  =  <p(t-x)  if  t  <t  B(R)  and  t  e  x  +  B(R),  and  ip(R.t-x)  = 

-<p(  t  -x )  if  t  €  B(R)  and  t  -  x  g  B(R).  Thus,  the  right  hand  side  of 
(5.5)  is 
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(5.6) 


|f(x)  I  |  v>(  t-x  )  I  dt  + 


o  V72 

|  <p(  t-x  )  |  ^dt  dx, 


supp  f 


where  X  _={t:t«B(R)  and  t  €  x  +  B(R)}  and  Y  ={t  :  t  €  BCR) 
x,  R  x,  K 

and  t  x  +  B  ( R ) } . 

If  R^.  =  sup  {  | x  1  :  x  e  supp  f }  ,  then  straightforward  calculations  show 
that  X  _.  £  B(R+R  )\B(R)  and  Y  £  B(R)\B(R-R  )  for  x  e  supp  f  and 

X  t  X\  I  X  |  H  I 

large  R.  Consequently,  we  can  make  the  estimate 


sup 

xesupp  f 


|<p(t-x)|  dt  +  | ( t — x )  |  dt 


(5.7) 


- 

<  sup  (|B(R+R  )l  -  | BCR) | )  sup 
xesupp  f  teB(R 

+  flB(R)  |  -  I  B(R-R  )  |1  sup  \s 

L  1  JteBCR)\B(R-Rf) 


I  <p(  t-x)  I ' 


R+R  )\B(R) 


i<p(t-x)  r 


<  a(<p,  R,  f )  (  I  B(  R+Rf )  |  -  |  B(R-Rf )  I  )  , 


where 

a(ip,R,f)  =  sup 


i<p(t-x) r 


xesupp  f'-  teB(R+R  )\B(R) 


teB(R)\B(R-R  ) 


Ip(t-x) r 


By  the  mean  value  theorem  we  have  the  estimate, 


(5.8) 


B(R+R  )|  -  |B(R-R  ) 


1/2 

f2dR„1 

1/2 

fR+R  ] 

< 

f 

f 

R 

V-  - 

R 

Combining  (5.6),  (5.7),  and  (5.8),  we  see  that  the  right  hand  side  oi 
(5.5)  is  bounded  by 
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d-1 

,  2dR  .  i/2rR+R  ,  2 

Isupp  f|  Hfll^  — j  —  a((p,R,f); 


and  this  will  tend  to  0,  thus  completing  the  proof,  once  we  show  that 
1/2 

a(</\R,f)/R  tends  to  0  as  R  — >  co. 

If  t  €  B(R+R^)\B(R)  and  x  e  supp  f  then  |t-x|  <  R  +  2R^,  and  if 
t  €  B(R)\B(R-R^J  and  x  e  supp  f  then  1 1  — x |  <  R  +  R^,.  Therefore, 
a(<p,R,f)  <  V2  i(R+2R^),  and  so 


v"2  i(R+2Rf)  ,  R+2Rf  .1/2 

£  (R+2Rf)'/2  I--5-- 


which  tends  to  0  by  hypothesis. 


Corollary  5. 3.  Given  the  hypotheses  of  Theorem  5.2.  If  f  e  LP(Rd)  and 

p 7  d  a 

<p  €  I,  (IR  ),  1  <  p  <  co,  then  the  tempered  distribution  f  is  a  well- 

2 

defined  element  of  L  (IR  )  and 

% 

(5.9)  If  Bo  „  ^  IlfILIML,. 

2'%  p  p 

Proof.  If  f  e  C^dR^)  then  (5.9)  is  a  consequence  of  (5.1)  and  Holder’s 

inequality.  (5.9)  extends  to  all  of  LP(Rd)  since  Cc(Rd)  =  LP(Rd), 

1  <  p  <  00;  as  such,  (5.9)  simultaneously  defines  f  and  gives  a  quantita¬ 
tive  norm  bound  (of  course,  f  is  known  to  exist  as  a  tempered  distribution 
since  the  elements  of  LP(R  )  are  tempered). 

a.  e.  d. 
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oo  d 

Remark  5.4.  a.  One  can  formulate  a  version  of  Theorem  5.  2  for  L  (!R  ) 


2  d 

instead  of  L7  (IR  ).  In  this  case  the  autocorrelation  can  be  defined  in 
loc 

terms  of  the  weak  *  <r(Lro(Rd),  L*(IRd))  topology.  As  such,  the  analogue  of 
Theorem  5.2  does  not  require  any  hypothesis  involving  i(R);  and  the  proof 
is  much  simpler  since  part  c  (of  the  proof)  can  be  replaced  by  a  simple 
estimate  and  application  of  Lebesgue  dominated  convergence,  cf.  ,  (Me;  B, 
89-92]. 

2  d 

b.  For  <p  e  L,  (IR  ),  the  positive  definite  approximant  P  of 
loc  <p,  R 

Section  5  should  be  compared  with  the  approximant 


Q<p,R(t)  |  B(  R) 


p(t+x)^>(x)dx 


B(R) 


of  Section  3.  First,  each  P  „  and  Q  R  >  0,  is  continuous,  and  the 

-  <P,R  ¥>,R 

bounds 


|FV,R(U  1  "  |  B( R) 


l?>(x)|  dx 


B(R) 


and 


IQ<p,R(t)l  "  | B ( R) 


‘B)l  ff, 


l<p(x)|  dxj 
B(R)  JB(R) 


H 


,  1/2 


T_t(/)(x)|  dx 


are  valid  for  all  t  e  IR  .  Second,  if  <p  satisfies  the  growth  condition  (in 


terms  of  i(R))  given  in  Theorem  5.2  then  {P  „  :  R  >  0}  and  (Q  „  :  R  >  0} 

-  R  (p,  R 

have  the  expected  similar  behavior  at  infinity.  More  precisely,  we  can  show 

that  if  lim  P  0  =  P  (resp.  ,  lim  Q  D  =  P  )  in  the  cr(M(Rd),  C  (!Rd)) 

R->co  '"•R  »  R,.  <"-R  * 

topology  then  lim  Q  =  P  (resp.,  lim  P  0  =  P  )  in  the  <r(M(Rd),  C  (Rd)) 
R-*oo  *'R  *  R-»oo  *>’R  * 

topology,  and,  in  either  case,  P^  is  positive  definite.  There  is  a 
corresponding  pointwise  result. 
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6.  An  elementary  restriction  theorem 


We  begin  with  a  straightforward  Fourier  transform  weighted  norm 
inequality  having  a  measure  weight. 


Proposition  6.1.  Given  p  e  M^+(R^)  for  which  p  =  P  e  (R^),  p  € 
Then 


(6.  1) 


Vf  e  L1  (IRd)nLp([Rd) ,  l|f|L  ,  ^ 


IIPIli^dlflLimL)^2  ^  ( ill  Pll*'2)(  Ill'll,  +  II f II  ) . 

p  1  p  ^  P  i  P 


where  (L1  ( !Rd ) nLP ( IRd ) ,  ||  •  *  - 1| .  +  ||  •  •  •  ||  )  is  a  Banach  space. 

1  P 


Proof.  The  Parseval  relation. 


(6.2) 


Vf  €  L 1  ( IRd ) . 


|f(r) I  duty)  = 


f *f ( t )P( t )dt , 


1  H 

is  valid  since  f*f(t)P(t)  e  L  (IR  ).  The  right  hand  side  of  (6.2)  is 


f 


f(x)  f (x-t)P(t)dt  dx  < 


|f(x)|  ||f  ||p  ||P||p/  dx. 


This  and  the  arithmetic-geometric  means  inequality  yield  the  result. 


The  following  is  our  elementary  restriction  theorem  with  explicit 
constant . 

Proposition  6.2.  Given  d  >  2  and  1  <  p  <  2d/(d+l).  Then  lip,  .11  , 
- -  d- 1  p 


[  1  ,  00  ]  . 


o.  e.  d. 


norm 


<  00 


(6.3) 


If (0)  I  dcr,  (0) 
a-1 


1/2 


+  Ilf  II 


PJ 


'd-1 


Proof.  Using  (4.1)  we  first  show  that  p^_^  e  (R0*)  for  1  <  p  <  2d/(d+l). 


-1/2, 


v  ,, 


'Vi11?'  =  Vi(2rr)1 


S  - >03, 

e.g.  , 

Lemma 

3.  1 ,  but, 

v  >  -1. 

Thus, 

we  can 

estimate 

( c'~2  % 

d-1  k  2 
r  r 


J  (2itr) 
d-2 


dr 


as  follows.  The  first  integral  is  finite  (for  any  p  >  1)  by  the  stated 
asymptotic  property  of  J^.  The  second  integral  is  finite  if  d  -  (dp'/2)  + 
(p'/2)  <  0  and  this  follows  if  1  <  p  <  2d/(d+l)  and  d  >  2.  The  case 
p  =  1  must  be  treated  separately  but  is  trivial. 

The  result  follows  from  Proposition  6.1. 


q.  e.  d. 


Remark  6.3.  a.  In  light  of  Theorem  4.1.  the  surface  measure  p^  ^  is  the 
power  spectrum  of  a  signal  <p,  e.g..  Definition  5.  1  and  Remark  5.4b.  Since 
the  norm  constant  in  Proposition  6. 2  is  explicit  and  computable  by  (4.1),  we 
see  that  (6.3)  provides  a  means  of  estimating  an  upper  bound  for  the  power  of 
<p  in  the  region  supp  f,  even  though  we  don’t  have  precise  knowledge  of  <p 
itself. 

b.  In  Proposition  6. 1  where  p  e  M^dR^)  is  given,  if  p  =  1  then  P 

is  always  an  element  of  L°°(R^ ) ,  and  if  p  =  2  then  P  e  L2(Rd)  so  that 
1  ^  d  2  *s  d 

p  <=  L  (R  )nL  (R  ).  In  the  case  p  =  p^  it  is  clear  that  p  =  ®  can  not 
be  used  in  Proposi tion  6.1  since  p^  i  L^(R^). 
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Section  IV 

Antenna-Polynomial  Problems 


Problems  on  Polynomials  with  Restricted  Coefficients 
Arising  from  Questions  in  Antenna  Array  Theory 

J.S.  Byrnes 


We  present  mathematical  formulations  of  some  classical  problems  in  antenna  design.  For 
references  concerning  these  problems,  see  the  paper  of  Newman  and  Giroux  in  this  volume. 

To  begin,  we  give  precise  mathematical  interpretations  to  two  standard  electrical  engi¬ 
neering  definitions,  that  of  maximum  sidelobe  level  and  beamwidth,  which  occur  in  the 
analysis  and  synthesis  of  a  line  array  of  equally-spaced,  identical,  omnidirectional  antenna 
elements.  As  is  well  known,  the  model  for  such  an  array  is  a  polynomial,  where  its  degree 
is  one  less  than  the  number  of  elements,  and  the  coefficients  are  the  weights,  or  shading 
coefficients,  of  the  array. 

Throughout  this  discussion,  let  z  —  e'e  lie  on  the  unit  circle  T,  let  P(z)  denote  a  poly¬ 
nomial  with  coefficients  a;,  at  least  two  of  which  are  nonzero,  and  let  Go  be  the  smallest 
nonnegative  value  of  6  where  the  maximum  of  )P(z)|  occurs.  Now  let 


Si  =  min{  6  :  6  <  0 0  and  |P(e‘s)|  is  increasing  on  (6,  do) }, 
82  =  max{  6  :  6  >  Go  and  |.P(e'*)|  is  decreasing  on  ( Go ,  £) }, 
crp  =  max  |/5(e,*)l,  and 

0€[0o  —  ,^2) 

Pp  =  62  —  S 1 . 


The  graph  of  (one  full  period  of)  |.P(e‘*)|  is  the  (modified)  beam  pattern,  ap  is  the 
(modified)  maximum  sidelobe  level,  and  Pp  is  the  (modified)  beamwidth  of  the  array.  We 
use  the  adjective  “modified”  because  we  have  omitted  the  change  of  variable,  normalization, 
and  expression  in  terms  of  decibels  that  one  finds  in  the  engineering  literature,  but  these 
differences  are  of  no  import  here. 


All  problems  presented  will  be  of  the  following  form: 

Determine  the  existence  and  construction  of  a  P  which  lies  in  a  certain  class,  whose  beam 
pattern  satisfies  certain  properties,  and  which  is  optimum  (or  nearly  so)  with  respect  to  one 
of  its  parameters.  These  classes  are: 


U  =  {  P  :  |ay|  =  1 }  (the  unimodular  polynomials), 
W0  =  {P:|ai|=lorO}, 

£  =  {P  :«i  =  ±i), 

£0  =  {  P  :  =  ±1,0}, 


T>k  =  {  P  :  max  |  — —  |  <  I\  }  (this  maximum  is  the  dynamic  range  of  the  array), 
Om 

)K,o  =  {P-  max  |-^-|<A’}1 
C  =  any  one  of  the  above  classes. 


Although  a  reasonable  bound  on  the  dynamic  range  is  ordinarily  important  (an  old  rule 
of  thumb  is  that  “K  =  2  and  everyone  is  happy,  K  =  10  and  some  are  happy,  I\  =  100  and 
nobody  is  happy”),  which  would  seem  to  render  zero  coefficients  fatal,  ♦  his  is  not  usually 
the  case.  Since  a  zero  coefficient  simply  means  that  that  element  is  tu.ned  off  (cither  by 
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choice  or  because  it  is  not  functioning  properly),  this  does  not  affect  the  dynamic  range 
in  a  meaningful  way.  Furthermore,  it  appears  that  in  many  cases  allowing  zero  coefficients 
enables  a  significant  reduction  in  the  degree  of  P  (see  problem  4  below). 

The  time  has  come  to  state  specific  problems: 


1.  Given  a  finite  (possibly  empty)  subset  S  of  F,  </  =  e'#°  £  S,  a  >  0,  0  >  0,  and  a  class  C, 
find  P  G  C  of  minimum  degree  (if  it  exists  at  all)  such  that  P{z)  =  0  for  all  c  £  S,  o  p  <  <r, 
and  0p  <  0. 


2.  Given  c  >  0  and  a  class  C,  find  P  €  C  of  minimum  degree  such  that 
||P(^)|-|P(Z2)||<e  for  all  z,,z2€r 

(a  so-called  u/tra-/Iat  polynomial). 

3.  Given  S  =  {e1®1 ,  elil , . . . ,  el9m)  nonempty,  c  >  0,  6  >  0,  and  a  class  C,  find  P  6  C  of 
minimum  degree  such  that  P(z)  =  0  for  all  z  £  S  and 

||P(2l)|  -  |P(z2)||  <  c  for  all  zt  =  e,T| ,  z2  = 

with  the  property  that  (ri,  t-2)  G  ( 0j  —  6, 0j  +  6)  is  empty  for  1  <  j  <  m  (a  notch  filter). 

It  is  appropriate  to  observe  that,  although  optimization  of  the  degree  is  called  for  in 
problems  1-3,  near-optimal  solutions  are  definitely  of  interest.  For  example,  the  problem  of 
most  practical  interest  is  1,  in  the  case  where  the  cardinality  of  S  is  very  small  (even  0,  1, 
or  2),  and  this  appears  mathematically  interesting  as  well.  Also  along  these  lines,  we  have: 


4.  Prove  or  disprove  the  conjecture:  The  P  €  U  of  minimum  degree  with  an  n-fold  l  H 
z  —  1  is 

n  <>-*’'>. 

rn  =  0 

In  fact,  even  a  proof  that  the  degree  of  such  a  P  must  be  exponential  in  n  would  be  of 
interest.  Note  that  if  the  class  U  is  replaced  by  Uo  (actually  even  £o)>  then  there  is  a  P  of 
degree  less  than  n3  satisfying  the  required  property. 

Analogous  problems,  when  one  of  the  other  parameters  aside  from  the  degree  is  optimized, 
are  of  interest  as  well.  In  this  case  there  is  an  additional  class  of  polynomials  to  be  considered: 

Mn  =  {  P  .  degree  P  —  n  and  J  ay  J ,  0  <  j  <  n,  are  preassigned  }. 

Thus,  another  typical  problem  is: 


5.  Given  q  =  e',°,  0  >  0,  the  degree  n  (often  between  10  and  40,  although  values  can  be 
as  much  as  a  few  thousand),  and  the  class  C  of  P,  find  P  such  that  0p  <  0  and  ap  is 
minimized. 

The  simple  job  of  stating  other  varieties  of  these  problems  is  left  to  the  reader.  Their 
solution,  a  much  more  daunting  task,  also  awaits  the  reader’s  efforts. 
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Section  V 

An  Ideal  Peak  Factor  Array 


An  ideal  omnidirectional  transmitting  array, 
and  optimal  peak  factor  array, 
for  less  than  half-wavelength  spacing 

James  S.  Byrnes  f 
Prometheus  Inc. 

21  Arnold  Avenue 
Newport,  RI  02840 


Abstract 

We  give  an  explicit  construction  of  an  ideal  omnidirectional  transmitting  array,  and  optimal  peak  factor 
array,  for  the  case  of  a  linear  array  of  identical  omnidirectional  elements  with  uniform  spacing  of  less  than 
half-wavelength.  The  construction  is  based  upon  the  Byrnes  polynomials,  introduced  by  the  author  in  1977. 

Introduction 

A  consideration  which  is  often  important  in  the  synthesis  of  an  antenna  pattern  is  the  peak  factor,  which 
is  the  ratio  of  the  peak  to  average  power  of  the  array.  For  a  line  array  of  equally  spaced  omnidirectional 
elements,  the  mathematical  model  for  the  pattern  is  a  polynomial  whose  coefficients  are  the  weights,  or 
shading  coefficients,  of  the  array.  The  quantity  to  be  synthesized  in  this  case  is  the  magnitude  of  P(z)  on 
the  unit  circle  C  =  {z  :  |z|  =  1},  and  the  peak  factor  is  now  the  ratio  of  the  sup  norm  to  the  L2  norm  of  P. 

Clearly  the  peak  factor  is  bounded  below  by  1,  and  the  classical  problem  is  to  make  it  as  close  to  1  as 
possible.  In  addition  to  arising  in  peak  power  limited  transmitting  and  other  aspects  of  antenna  design,  the 
identical  question  occurs  in  digital  filtering.  Furthermore,  the  same  polynomial  problem  has  been  the  object 
of  intense  study  by  mathematical  analysts  for  more  than  fifty  years,  including  such  notable  mathematicians 
as  P.  Erdos  [3],  G.H.  Hardy,  J.P.  Kahane  [4],  T.  Korner  [5],  J.E.  Littlewood  [6],  and  D.J.  Newman  [7]. 

A  parallel  concern,  to  engineers  and  mathematicians  alike,  is  the  question  of  synthesizing  various  antenna 
patterns  (i.e.,  constructing  polynomials  with  specified  moduli  on  the  unit  circle)  when  certain  restrictions 
are  placed  upon  the  coefficients.  When  transmitting,  for  example,  one  usually  wants  to  maximize  the  total 
power  output  of  the  array,  which  is  achieved  when  each  individual  element  broadcasts  at  full  power.  Thus, 
in  this  case,  the  coefficients  must  have  the  same  magnitude.  In  addition,  the  usual  peak  power  limitation 
means  that  the  output  power  must  be  as  close  to  contact  as  possible  in  all  directions,  i.e.,  |P(z)|  should  be 
as  close  to  constant  as  possible  on  C.  The  fundamental  result  in  this  specific  area  is  that  of  Kahane  [4]: 

Theorem.  There  is  an  absolute  constant  c  >  0  such  that  for  each  positive  integer  n,  there  is  a  polynomial 
P(z)  of  degree  n  with  unimodular  coefficients  which  satisfies 

|P(z)|  =  y/ii  +  E  for  all  z  6  C,  where  (El  <  cn™  log n.  (1) 

Since  any  such  P  has  n  +  1  coefficients  of  modulus  1,  so  that  its  L7,  norm  is  obviously  \/n  +  1,  this  shows 
that  the  ideal  peak  factor  of  1  is  indeed  achievable  asymptotically,  even  when  the  coefficients  are  required 
to  have  the  same  magnitude.  Note  that  Kahane’s  theorem  shows  that  the  ratio  of  the  maximum  modulus  of 
P  to  the  minimum  modulus  of  P  can  be  asymptotically  I  as  well,  a  stronger  result  than  the  peak  factor  one 
just  mentioned.  As  an  aside,  we  mention  that  the  corresponding  problem  when  the  coefficients  are  required 
to  be  ±1  is  one  of  the  important  unsolved  problems  in  this  area  of  mathematical  analysis. 

Although  Kahane’s  beautiful  result  was  greeted  with  much  enthusiasm  by  the  mathematical  community, 
it  has  not  found  application  in  the  engineering  problems  discussed  above.  This  is  because  his  proof  is 
an  existence  one:  he  shows  probabilistically  that  such  polynomials  must  occur.  Thus  far  such  functions 
have  not  been  actually  constructed.  This  too  is  an  important  outstanding  question,  in  both  mathematics 

t  The  author  is  also  Professor  of  Mathematics  at  the  University  of  Massachusetts  at  Boston 
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and  engineering.  The  purpose  of  this  note  is  to  give  an  explicit  construction  of  an  ideal  omnidirectional 
transmitting  array,  and  optimal  peak  factor  array,  for  the  case  of  a  linear  array  of  identical  omnidirectional 
elements  with  uniform  spacing  of  less  than  half-wavelength.  In  certain  circumstances  such  close  spacing  can 
cause  a  notable  increase  in  mutual  coupling.  This  physical  problem  will  not  be  addressed  here. 


Results 

As  is  well  known,  for  a  linear  array  of  n  equally  spaced  identical  omnidirectional  elements  the  array 
factor  is  the  polynomial 


P(z)  =  *  =  e 

i=o 


2  xiu 


d  sin  6 


u  = 


(2) 


Here,  the  ay  are  the  shading  coefficients,  d  is  the  array  spacing,  A  is  the  wavelength,  and  0  is  the  angle 
of  incidence  (0  =  0  is  broadside). 

Note  that  for  d  =  A/2,  generally  considered  the  ideal  case,  2iru  goes  from  —n  to  ir  as  sin  0  goes  from  —1 
to  1,  so  that  z  traverses  C  completely.  However,  if  d  <  A/2,  then  a  portion  of  C  will  be  omitted  by  z.  This 
is  precisely  the  property  that  we  now  exploit. 

In  fact,  the  construction  follows  immediately  from  earlier  work  of  the  author  [2],  The  crucial  property 
of  the  polynomials  introduced  there,  and  later  called  Byrnes  polynomials  [1,  5],  is  given  by  part  (if)  of  the 
basic  theorem  in  [2]: 

Theorem.  For  n  a  positive  integer,  let 


n— 1 n-1 

P(z)  =  ^2  exp(2irijkn~1)P+kn ,  z  =  e2w,u . 

k—Q j=0 


(3) 


Then  for  any  e,  n  1  <  e  <  ^ ,  we  have  \P(z)\  =  n  +  E  for  all  u,  e  <  |u|  < 
where  |E|  <  1  +27r_1  +5(7rc)_1. 

To  clarify  the  significance  of  this  result  in  the  present  context,  suppose  for  simplicity  that  d  =  A/4,  so 
that  z  is  now  traversing  exactly  1/2  of  C  as  sin0  goes  from  -1  to  1.  The  theorem  immediately  yields  the 
correct  choice  of  coefficients,  gotten  by  merely  changing  z  to  —z  in  (3).  Hence,  the  optimal  polynomial  is 


n-ln-l  . 

P(z)  =  J2Y.(~iy+knex  P(2*«j*n“V+t"  z  =  e2xiu ,  u=^j-  .  (4) 

*=07=0  4 

This  gives  |P(z)|  =  n  +  E,  |£|  <  1  +22jr-1  <  9  for  all  6.  Note  that  the  degree  of  P  is  n2  —  1,  so  that  its 
L2  norm  is  n  and  it  represents  an  array  with  n2  elements.  Furthermore,  the  error  E  is  uniformly  bounded 
independent  of  n,  so  that  we  indeed  have  an  optimal  peak  factor  array.  In  addition,  the  coefficients  of  P 
are  unimodular,  so  that  it  also  represents  an  ideal  omnidirectional  transmitting  array.  These  properties  arc 
illustrated  in  Figures  l.a-l.e,  showing  polar  plots  of  |F(z)|  in  Equation  (4)  for  various  choices  of  n  and  for 
—  7r  <  9  <  v.  Observe  that  for  small  n  the  actual  array  patterns  exhibit  an  error  that  is  considerably  smaller 
than  that  obtained  for  the  general  case  considered  in  [2].  Finally,  the  spread  (i.e.,  the  difference  between  the 
maximum  and  minimum,  expressed  in  dB’s)  of  |P(z)|  as  a  function  of  the  number  n2  of  elements  is  shown 
in  Figure  2.  This  graph  clearly  illustrates  the  desired  flatness  of  |P(2)|,  as  described  above. 
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Section  VI 

Polynomials  with  Unimodular  Coefficients 


Properties  on  the  Unit  Circle 
of  Polynomials  with  Unimodular  Coefficients 

By  Donald  J.  Newman  and  Andre  Girouxt 
Prometheus  Inc. 

21  Arnold  Ave. 

Newport,  RI  02840 


The  following  problem  is  posed  in  [1]:  given  the  magnitude  of  the  coefficients  of  a  polynomial  P,  a  finite 
subset  S  of  the  unit  circle  C,  and  a  point  p  on  C  distinct  from  those  in  S,  choose  the  phases  of  the  coefficients 
so  that  P(z)  =  0  for  all  z  in  S,  the  maximum  on  C  of  |P(z)|  occurs  at  z  =  p,  and  the  maximum  of  |P(z)| 
on  a  subset  of  C  excluding  an  appropriate  interval  around  p  (the  “beamwidth”)  is  as  small  as  possible. 
As  explained  in  [1],  this  problem  arises  naturally  in  linear  antenna  theory,  filter  theory,  and  other  classical 
electrical  engineering  applications. 

Our  main  result  is  encompassed  in  the  following  theorem: 

Theorem.  Given  o-i ,  a?, . . . ,  an ,  /?  on  C,  /?  ^  a,-,  there  exists  a  polynomial  P  with  coefficients  of  modulus 
one  such  that  the  zeros  of  P  on  C  are  precisely  the  points  a  1 , 02, . . . ,  a„  and  that 

||P||  =  max{|P(z)|:z€C}  =  |P(/?)|. 

(The  points  aj  are  not  assumed  to  be  distinct:  multiple  zeros  are  allowed.) 

We  shall  prove  the  theorem  in  two  steps,  first  constructing  a  polynomial  with  coefficients  of  modulus  one 
with  the  prescribed  zeros  (lemma  1)  and  then  modifying  it  so  that  it  also  has  the  prescribed  maximum 
(lemma  4). 

Let  U  denote  the  (non-linear)  class  of  polynomials  with  coefficients  of  modulus  one.  We  introduce  in  U  the 
operation  of  “encapsulation”  which  we  denote  by 

Definition.  Let  M  be  the  degree  of  P.  Then 

(P®Q)(z)  =  P(z)Q(zm+1). 

Explicitly,  then,  if  P(z)  =  X^tio0*2*  anc*  Q(z)  —  X^lo  bjZ* ,  ^cn  P  ®  Q(z)  =  C|2«,  where 

cj(M+\)+k  =  atbj  for  0  <  k  <  M  and  0  <  j  <  N. 

From  this  expression,  it  is  clear  that  ®  is  associative  and  that  U  is  indeed  closed  under  it.  Also,  if  Q  does 
not  vanish  on  C,  then  the  zeros  of„P  ®  Q  are  precisely  those  of  P. 

Lemma  1  .  Given  Oj ,  02,  • . . ,  otn  on  C,  there  exists  P„  in  U  which  vanishes  at  the  a;-  and  nowhere  else  on 
C. 

Proof.  Let  u  =  exp2<n'/3.  We  shall  prove  a  little  more,  namely  that  such  a  polynomial  can  be  found 
with  constant  term  one  and  leading  term  — w.  We  use  induction  on  n.  Let 

Qi(z)  =  1  -  z/<*i 


and  let 

P\(z)  =  Q\(z)  if  oiw  =  1 

and 

Pi(z )  =  Qi(z)  ®  (1  +  z  +  OjWZ2) 

=  1  —  (1/0  l)z  +  Z2  —  (1/oj  )z3  +  Q]  usz4  — wz5 

f  The  authors  are  also  with  Temple  University  and  the  University  of  Montreal  respectively. 
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otherwise  (the  condition  aiu  ^  1  also  ensures  that  the  polynomial  1  +  z  +  aiwz2  ^  0  on  C).  Assume  now 

that  R(z)  =  1  H - wzN  is  a  polynomial  in  U  of  which  the  zeros  on  C  are  precisely  the  points  fij  =  Qj/an 

for  1  <  j  <  n  —  1 .  Set 

S(z)  =  R(z)(l  -  zN  +  z2N  -  z3N+1) 

=  1  H - f  uz4N+1. 

To  verify  that  S  is  in  U,  we  need  only  check  the  modulus  of  the  coefficients  of  zN  and  z2N .  But,  the 

coefficient  of  zN  is  — w  —  1  =  u2  and  that  of  z2N  is  u  +  1  =  —  w2.  We  now  claim  that 

1  -  zN  +  z2N  -  z3N+l 

does  not  vanish  c  n  C  except  for  a  simple  zero  at  z  =  1.  Indeed,  for  the  sum  of  four  points  on  C  to  vanish, 
it  is  necessary  and  sufficient  that  they  cancel  in  pairs.  There  are  thus  three  possibilities: 

a)  1  —  zN  =  0  and  z2N  —  z3iV+1  =  0;  this  gives  zN  =  1  and  zN+1  =  1,  hence  z  —  1; 

b)  1  +  z2N  =  0,—  zN  —  z3N+1  =  0;  this  implies  z2N  =  — l,z2^+1  =  —1,  which  is  impossible; 

c)  1  —  z3N+l  =  0,  —  zN  +  z2!V  =  0;  in  this  case,  z3/^+1  =  1  and  zN  =  1,  therefore  z2N+l  =  zJV+1  =  zN  =  1, 

hence  z  —  1  again. 

Setting 

Qn(z)  =  S(z/an)  =  1  +  •  •  •  +  (u/a4nN+1)ziN+1 
we  obtain  a  polynomial  in  U  vanishing  on  C  precisely  at  the  points  oi, c*2»  •  •  • ,  an.  Then 

pn(z)  =  Qn{z)  if  a^+1  =  -l 


and 

Pn(z)  =  Q„(z)  ®  (1  +  z  —  a^+1z2)  otherwise. 

This  completes  the  induction  and  the  proof  of  the  lemma. 

Lemma  2.  If  n  >  3,  the  polynomial  pn(z)  =  l  +  z-)-z2  +  --  -  +  zn~l  —  zn  +  zn+1  + - K  z2n~l  does  not 

vanish  on  C  and  max{|p„(z)| :  z  €  C}  is  attained  if  and  only  if  z  =  1. 

Proof.  Since 

z2n _ 1 

pn{z)  =  ~-±-2zn 

we  see  that 

pn(e,'2*)e-'(3"-1)*  =  Sin2—  -  2e‘* 
sin  9 

so  that  p„(e'29)  =  0  implies  that  e'9  is  real  which  implies  in  turn  that  pn(e'2S)  =  Pn(l)  =  2rz  —  2,  a 
contradiction. 

Our  statement  about  the  maximum  modulus  property  of  pn(z)  follows  from  the  fact  that  pn2{z)  has  positive 
coefficients. 

Lemma  3.  If  n  >  0,  the  polynomial  qn(z)  =  l  +  z  +  z2  +  --  -  +  zn  —  zn+1  does  not  vanish  on  C. 

Proof.  One  has  (1  —  z)qn{z)  =  1  —  2zn+1  +  zn+2,  which,  by  the  triangle  inequality,  can  vanish  only  if 
l,zn  +  1  and  zn  +  :  are  in  the  same  direction;  that  is,  here,  z  =  1.  But,  obviously,  9„(1)  ^  0. 

Lemma  4.  If  P  is  any  polynomial  such  that  P{0)  ji  0,  then  there  exists  a  polynomial  T  in  U  which  does 
not  vanish  on  C  and  is  such  that  P<S>T  attains  its  maximum  modulus  on  C  at  z  =  0. 
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«? 


Proof.  We  can  assume  that  0=1.  let  T  =  qm  ®  r  ®  p„  where  qm  and  p„  are  as  in  lemmas  3  and  2 
respectively,  and  r(z)  =  1  +  e'ez  —  z 2  with  9  ^  ±ir/2:  r(z)  does  not  vanish  on  C  and  Imr'(l)/r(l)  =  2sin  9. 
Consider  first  P *  =  P  <S>  qm  ®  r-  If  M  is  the  degree  of  P ,  then 


P*'(l)  _  P'(  1) 
P*(l)  ~  P(l) 


+  (M  +  1) 


^(1) 

?m(l) 


r'fl) 

+  (M  +  l)(m  +  2)--^i 


and  therefore 


Im 


™,Imrn+0+(M+1Km+2)Imrn- 


p'(  i) 

=  Im  !■  +  (M  +  l)(m  +  2)2sin  9. 


P(  1) 


A  proper  choice  of  m  and  6  will  cancel  this  expression,  making  zero  a  stationary  point  of  the  function 
|P*(e'*)|2.  Since  this  function  does  not  vanish  at  9  =  0,  this  stationary  point  can  be  transformed  into  the 
maximum  point  over  all  9  by  multiplying  P*  by  a  delta-like  function,  in  our  case  by  encapsulating  it  with 
pn  for  suitably  large  n.  Indeed,  one  has 


»  /  ,  sin  tx9  _  .  „ 

w*  >1  =4-4ii^75cosS/2  + 


By  computing  the  second  derivative,  it  is  easily  seen  that  pn{9)  falls  off  like  e~cnB* ,  and  this  shows  that  it 
is  an  adequate  delta  function  for  our  purposes. 

This  completes  the  proof. 


Reference. 

[1]  James  S.  Byrnes,  Donald  J.Newman,  “Null  Steering  Employing  Polynomials  with  Restricted  Coeffi- 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
OVERVIEW  OF  PRESENTATION 

OVERVIEW  OF  PRESENTATION 

•  Sensor  array  signal  processing  —  problem  formulation 

•  Various  objectives  and  approaches  —  real-time  versus  off-line 

•  Problem  complexity  and  realistic  simplifying  assumptions 

•  Antenna  array  weight  design 

•  AWD  -  the  program 

•  Examples 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
SENSOR  ARRAY  PROCESSING 


General  Problem  Characteristics 


•  Inputs  axe  discrete  samples  of  a  continuous  function  of  one  or  more 
independent  variables 

—  digital  filtering  of  analog  signals 

—  phased  array  radars 

•  The  function  being  sampled  generally  consists  of  a  superposition 
of  multiple  signals  of  interest  and  noise 

—  sinusoids  in  analog  signals  -  time  series  analysis 

—  plane  waves  impinging  on  an  array  -  sensor  array  process¬ 
ing 

•  Signals  may  be  functions  of  parameters  of  interest 

—  frequencies  of  sinusoids 

—  directions-of- arrival  (DO As)  of  plane  waves 

•  Reconstruction  of  signals  as  functions  of  time  may  be  of  primary 
interest 
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SENSOR  ARRAY  PROCESSING 
OBJECTIVES  AND  APPROACHES 


Objectives 

•  Obtain  parameter  estimates  in  real-time 

—  track  time- varying  frequencies  of  sinusoids  in  noise 

—  track  moving  targets  -  active/passive  radar  tracking 

•  Obtain  optimal  parameter  estimates  off-line 

—  trajectory  reconstruction  -  error  analysis 

•  Reconstruct  one  or  more  signals  individually  in  multiple  signal 
environment  in  real-time 

—  establish  reliable  communication  links 

—  passive  listening  in  multiple  source  environments 

Approaches 


•  Off-line 

—  many  approaches  possible 

—  choice  depends  on  computational  cost  and  available  re¬ 
sources  ( cf ’  maximum-likelihood) 

•  Real-time 

\ 

—  adaptive  array  processing 

—  pseudo -adaptive  array  processing 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
PROBLEM  FORMULATION 


Problem  Statement 


•  GIVEN 

—  sensor  array  of  known  fixed  geometry, 

—  statistics  (second-order)  of  noise, 

—  parameters  of  signals  of  inteiest  ( e.g .  DOA,  center  fre¬ 
quency) 

—  parameters  of  unwanted  signals  (e.g.  jammers) 

•  FIND 

—  a  linear  combination  (weight  vector/linear  functional)  of 
sensor  element  outputs 

that 

•  MINIMIZES/MAXIMIZES 

•  —  a  cost/merit  function  of  the  combined  array  output 

•  SUBJECT  TO 

—  various  design  specifications. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
SIMPLIFYING  ASSUMPTIONS 


Simplifying;  Assumptions 


•  Medium  of  transmission  is  assumed  to  be  isotropic  and  non- 
dispersive. 

•  Array  consists  of  sensor  elements  with  fixed  relative  locations. 

•  Signal  sources  are  assumed  to  be  farfield  and  lie  in  the  plane  of 
the  array. 

•  Sensor  elements  are  omni-directional. 

•  Elements  have  constant  gain  inside  a  predetermined  (finite)  fre¬ 
quency  band  of  operation  and  negligible  gain  outside  such  band. 

•  Noise  present  in  every  sensor  output  is 

—  additive,  zero-mean,  and  stationary, 

—  uncorrelated  with  the  signals  of  interest  (SOI). 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
PROBLEM  FORMULATION 


Combined  Array  Output 


n  i 
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PROBLEM  FORMULATION 


ANTENNA  ARRAY  WEIGHT  DESIGN 
PROBLEM  FORMULATION 


Basic  Definitions 

•  Complex  weight  vector  (to  be  designed) 

W\ 

w2 

w  — 

wN 

•  The  noise-free  sensor  output  vector 

si 
S2 

Sn 

•  The  noise  vector 

n\ 
n2 

nN 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
PROBLEM  FORMULATION 


Basic  Equations 


•  Noise-free  output 


•  Noise  output 


•  Total  output 


•  Total  noise  power 


S(0,  cj)  =  w*s(@,  lj) 


Sn  =  w*n 


z(0,u)  =  S(0,u)  +  Sn 


where 


Pn  =  E  <?i\wi[  = 

i=i 


E  =  diag((ii,  (72, ,  ctn)  =  E(nn*) 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
PROBLEM  FORMULATION 


Sensor  Element  Output 


•  The  ith  sensor  element  output  can  be  written  as 

Si{6 ,  lj)  =  exp(j»27r/(xz-  cos  6  +  y*  sin  9)) 

•  /  is  the  normalized  (dimensionless,  /actual //nominal)  signal  fre¬ 
quency, 

•  6  is  the  incident  angle  (degrees), 

•  Xi  is  the  z-location  (in  units  of  Anom;nai)  of  element  z,  and 

•  yi  is  the  y-location  (in  units  of  Anomjnai)  of  element  i. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
REAL-WORLD  PROBLEMS 


Real-World  Problem  Characteristics 


•  Signals  of 

—  unknown  number, 

—  unknown  directions  of  arrival  (DOAs), 

—  unknown  frequencies. 

•  Modeling  errors 

—  unmodeled  signal  environment  (e.g.  reflectors,  echo,  etc.), 

—  inter-element  effects, 

—  anisotropic  media, 

—  nonlinearities. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
HANDLING  REAL-WORLD  DESIGN 


On-Line  Adaptive  Methods 


•  Maximizes  (minimizes)  a  certain  merit  (penalty)  function  (e.g. 
SNR)  of  the  combined  array  output. 

•  Rely  on  some  knowledge  of  SOI, 

—  DOA,  frequency, 

—  functional  form. 

•  Difficult  to  take  advantage  of  apriori  knowledge  of  interfering  sig¬ 
nal  frequencies,  directions-of-arrival,  etc. 

•  Designed  to  be  computationally  efficient  and  very  fast. 

•  Basically  places  pattern  nulls  in  interference  DOA’s.  Thus,  is  quite 
sensitive  to  interference/signal  correlation. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
HANDLING  REAL-WORLD  DESIGN 


Deterministic  Methods 


•  Assumes  knowledge  of  frequencies  and  directions  of  arrival  of  sig¬ 
nals  of  interest  and  interference. 

•  Computationally  demanding 

—  slow  compared  to  adaptive  methods, 

—  unsuitable  for  on-line  compuation. 

—  pre-storage  of  off-line  results  is  necessary. 

•  Merit  functions,  pattern  nulls,  pattern  constraints  and  array  gains 
are  arbitrary. 

•  Sensitivity  to  data  can  be  reduced  as  desired  by  clever  choice  of 
design  constraints  and/or  the  merit  function. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
HANDLING  REAL-WORLD  DESIGN 


Combining  Adaptive/Deterministic  Methods 


•  Adaptively  estimate  signal/interference  parameters  and  obtain  an 
initial  set  of  weights. 

•  Select  an  appropriate  merit  (penalty)  function  and  a  desired  set 
of  desing  constraints,  and 

—  run  a  suitable  deterministic  scheme  ( e.g .  our  scheme), 

—  obtain  optimal  weight  vector, 

—  such  weights  produce  desired  output. 

•  One  suggestion  to  use  results  of  deterministic  schemes  on-line  is 
to 


—  devise  a  combination  possible  scenarios, 

—  use  off-line  computations  to  find  corresponding  optimal 
weight  vectors, 

—  store  resulting  weight  vectors  in  memory, 

—  using  on-line  data,  determine  most  appropriate  weight 
vector  to  use. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
DETERMINISTIC  ANTENNA  WEIGHT 

DESIGN 


Closed  Form  Methods 


•  Usually  old  methods  (40’s  and  50’s). 

•  Easy  to  implement  (hand  calculations,  charts,  . . . ) 

•  Restrictive  in  nature. 

•  Fast. 

•  Example:  Dolph’s  method. 

Iterative  Methods 

•  Newer  methods.  Uses  currently  available  computing  power. 

•  Implementation  requires  intensive  computations. 


•  Less  restrictive. 

•  Much  slower. 

•  Example:  our  method. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
A  HISTORY  OF  DETERMINISTIC 

METHODS 


Closed  Form  Methods 


•  Stone  (U.  S.  Patents  1,643,323  and  1,715,433): 

—  A  binomial  expansion  representaion. 

—  Design  of  patterns  by  choice  of  coefficients. 

•  Schelkunoff  (1943): 

—  Polynomial  representation  of  arrays  of  equally  spaced  an¬ 
tenna  elements. 

—  Design  of  patterns  by  cleverly  placing  zeroes  of  the  repre¬ 
senting  polynomial. 

•  Dolph  (1946): 

—  Design  using  Tchebychev  polynomials. 

—  Obtained  optimal  tradeoff  between  sidelobe  gain  and 
mainlobe  width. 

—  First  to  talk  about  solutions  that  are  optimal  in  some 
sense. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
A  HISTORY  OF  ITERATIVE  METHODS 


Non— convex  Design 

•  Cheng  and  Tseng  (1965,1967): 

—  Optimization  of  Hermitian  ratios. 

—  Unconstrained  case:  generalized  eigenvalue  problem. 

—  Constrained  case:  iterative  solutions. 

—  Method  actually  known  before.  Others  used  same  meth¬ 
ods.  For  example,  Butler,  et  a/.,  (1971,1972). 

•  Waren,  et  al, (1966, 1967): 

—  General  mathematical  programming  approach. 

—  Optimization  of  magnitude  and  phase  of  weights. 

—  Design  of  array  geometry. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
A  HISTORY  OF  ITERATIVE  METHODS 


Convex  Design 


•  McMahon,  et  al,  (1972): 

—  Optimize  real  and  imaginary  parts  of  weights  taken  sepa¬ 
rately. 

—  Translate  problem  into  linear  program. 

—  Fake  a  “tuned”  objective  funtion. 

•  Wilson  (1976): 

—  Similar  to  above  approach. 

—  Uses  a  Tchebychev-norm  approximation  technique. 

•  Evans  and  Fortmann  (1975):  a  functional  programming  approach. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 

CONVEXITY 


Basic  Definitions 


Definition  1  (Convex  Set)  (]  C  V  is  said  to  be  a  convex  set  if 
fix  +  (1  —  fi)y  E  fi  whenever  x,  y  E  fi  and  0  <  fi  <  1. 

Definition  2  (Convex  Function)  i4  function  /  :  V"  — ►  is  a 
convex  function  on  a  convex  set  Cl  C.  V  if  f(fix  +  (1  —  /i)y)  < 
fif(x)  +  (1  —  /i)/(y)  /or  a//  £,y  E  fi  and  0  <  fi  <  1. 

Definition  3  (Convex  Program) 

minimize  <f>(u>) 
subject  to  cu  E  fi 

where  <1 >  :  F  — >  3ft  zs  a  convex  function  and  fi  25  a  convex  set. 


* 


V 1 1  -  2  0 


ANTENNA  ARRAY  WEIGHT  DESIGN 

CONVEXITY 


Basic  Facts 


Fact  1  Given  the  convex  program  described  above,  and  given  that 
the  set  ft  is  non-empty  and  closed,  then  the  set  of  points,  Ho  Q  ft? 
that  solves  the  convex  program  is  a  non-empty  convex  set. 

Fact  2  Let  f\  :  V  — ►  and  •  V  — ►  be  a  convex  and  concave 
functions  (respectively)  on  a  convex  set  ft  E  V,  then  the  sets 

•  {w  G  ft|/i(u>)  <  b},  and 

•  {we  ft  1/2(0;)  >  b} 

are  convex  for  any  b  e  3?- 
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ANTENNA  ARRAY  WEIGHT  DESIGN 

CONVEXITY 


Basic  Theorems 


Theorem  1  Let  {/i,  /2, . . . ,  fn},  be  a  set  of  real  valued  convex  func¬ 
tions  on  a  convex  set  Cl  C  V ,  then  the  following  functions  are  convex 
on  Cl. 

•  Fi(x)  =  max  /,-(x). 

1 

•  F2(x)  =  £”=1  ctifi(x),  where  az-  >  0. 

Theorem  2  Let  f  :  S  x  Cl  —>  5R,  be  a  real  valued  function,  where 
Cl  cy  is  a  convex  set  and  E  C  U  is  any  subset  ofU.  If  the  function 
fa(u>)  =  f(a,u)  is  convex  on  Cl  for  any  fixed  a  £  E,  then  the  function 
F  :  Cl  — >  defined  by 

F(u)  =  sup  /(cr ,  uj) 

<7€E 

is  convex  on  Cl. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 

CONVEXITY 


Some  Convex  Functions  of  Weight  Vector 


•  Total  noise  power,  Pn  =  w*T,w. 

•  Real  and  imaginary  parts  of  a  certain  (scaled)  element  weight, 
awi. 

•  Magnitude  of  a  certain  (scaled)  element  weight,  aw{. 

•  Real  and  imaginary  parts  of  scaled  noise-free  pattern  function, 
aS(6,uj). 

•  Scaled  array  noise-free  pattern  gain,  |a5(0,a;)|. 

•  max  1 5(0,  a;)  I  for  any  fixed  9. 

u>i<u<u>2 

•  max  |S(0,u;)|  for  any  fixed  u. 

6\<6<.62 

•  max  |5(0,  o;)|  are  also  convex  for  any  6 1,  #2?  and  ujo- 

Oi<O<02,ui<uj<u)2 

•  Any  positive  combination  of  the  above  functions. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 

CONVEXITY 


Why  Convexity 


•  Efficient  algorithms  for  solving  convex  programs  exist.  That  is, 
such  algorithms 

—  are  fast, 

—  are  numerically  stable, 

—  can  detect  whether  an  optimal  solution  exists  or  not. 

•  Locally  optimal  solutions  are  globally  optimal. 

•  Additional  information,  such  as  the  Lagrange  multipliers  (dual 
prices)  are  extremely  useful  in  judging  a  design. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
GENERAL  CONVEX  WEIGHT  DESIGN 


General  Outline 


•  Reduce  weight  design  problems  to  convex  programs. 

•  Approximate  general  convex  programs  by  ones  for  which  efficient 
algorithms  exist. 

•  Solve  approximate  convex  programs. 

•  Report  optimal  weights  and  shadow  prices  (Lagrange  Multipliers). 
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-H 


(feasibility, 

Lagrange  multipliers,  etc. .  .  ) 


Simulation 

I 


_ J 


VII-26 


AWDSL  Program 


Solution  t(/0pt 


Figure  1:  Our  implementation  of  awd,  a  specification  language  compiler  and  solver. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
GENERAL  CONVEX  WEIGHT  DESIGN 


Questions  and  Answers 


Question:  What  antenna  weight  design  problems  can  be  reduced  to  con¬ 
vex  programs? 

Answer:  A  very  general  class  of  weight  design  problems  can  be  cast  as 
convex  programs. 

Question:  Can  large  problems  be  translated  to  convex  programs  in  real 
time? 

Answer:  Automating  the  translation  step  from  a  design  problem  to  a 
convex  program  using  computers  rather  than  having  it  manually 
done,  makes  it  possible  to  handle  large  problems. 

Question:  How  closely  does  the  approximate  program  represent  the  orig¬ 
inal  problem? 

Answer:  If  the  translation  is  accurately  done,  the  approximate  program 
can  be  made  arbitrarily  close  to  the  original  problem. 
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ANTENNA  ARRAY  WEIGHT  DESIGN 
GENERAL  CONVEX  WEIGHT  DESIGN 


Questions  and  Answers 


Question:  Can  the  approximate  program  be  solved  effectively 
and  are  such  solutions  reliable? 

Answer:  If  the  approximate  program  is  a  linear  or  quadratic  pro¬ 
gram,  then  there  are  fast,  reliable  and  effective  algorithms 
to  find  solutions. 

Question:  Can  the  mathematical  program  solution  be  related  to 
the  original  design  problem? 

Answer:  An  inverse  translation  is  necessary  to  relate  the  opti¬ 
mization  results  to  the  original  problem. 

Question:  Are  there  any  extra  advantages  that  warrant  the  ex¬ 
tra  effort? 

Answer:  Additional  information  available  at  the  end  of  the  op¬ 
timization  (e.g.  Lagrange  multipliers)  is  as  important  as 
the  actual  solutions. 

Also,  if  such  method  fails,  then  no  other  method  can  suc¬ 
ceed  in  solving  the  problem.  Non-convex  iterative  methods 
cannot  assert  such  a  powerful  claim. 

Question:  Where  does  the  proposed  method  fit  in  the  typical 
engineering  design  cycle? 

Answer:  To  be  discussed  later. 
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A  Proposed  Approach  To  Convex  Design 


•  Introduce  a  specialized  specification  language,  called 
Antenna  Weight  Design  Specification  Language 
or  AWDSL,  to  facilitate  the  specification  of  a  convex 
weight  design  problem. 

•  Introduce  a  specification  language  compiler  to 

—  automate  translating  a  design  problem  into  a  con¬ 
vex  mathematical  program,  and 

—  link  the  solutions  and  whatever  extra  information 
resulting  from  optimization  to  the  original  design 
problem. 

•  Choose  a  class  of  convex  optimization  problems  which  can 
arbitrarily  approximate  the  original  design  problem  and 
for  which  effective  algorithms  exist. 

•  To  facilitate  implementing  optimization  techniques,  think 
of  a  complex  weight  vector  as  a  pair  of  real  vectors. 

w  =  Wsr  +  jw% 

where 

are  the  real  and  imaginary  components  of  the  complex 
weight  vector. 
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Our  Proposed  AWDSL  Structure 


array_description  { 

description  of  array  geometry 
and  element  characteristics 
•  •  • 

> 

minimize  { 

design  objectives 
•  •  • 

} 

subject_to  { 

design  constraints 
•  •  • 

> 
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The  Array  Description  Section 


array_description{ 

/*  number  of  antenna  elements  in  array  */ 
n_elements  =  number  of  elements ; 

/*  lower  and  upper  bounds  on  frequency  band  * 
lower_freq  =  lower  bound; 
upper_freq  =  upper  bound; 

/*  specify  the  x  locations  */ 

element_x_loc  (0  =  x-loc  of  element  i; 

/*  specify  the  y  locations  */ 

element_y_loc(0  =  y-loc  of  element  i; 

/*  specify  the  element  noise  powers  */ 
element_pwr(z)  =  power  of  element  i; 

/*  specify  the  initial  weights  */ 

element_w(0  =  {real  part  ,  imaginary  part)  ; 

} 
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The  Objective  Section 


An  objective  term  is  of  the  form 
expr  *  g; 

where  expr  is  any  positive  expression  and  g  is  one 
of  the  following  convex  functionals 

•  noise_pwr, 

•  max_mag_S(0i,  02,  i), 

•  max_mag_S(0i,  @2,  j)(/),  or 

•  max_mag_S(0iJ  02,  j)(/i,  /2,  &). 
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The  Constraints  Section 


Constraints  are  in  the  form  of  functional  inequali¬ 
ties  such  as 

g  <=  expr  ; 
expr  >=  g  ; 
g  >=  expr  ; 
expr  <=  g  ; 
g  ==  expr  ; 
expr  ==  g  ; 

exprl  <=  g  <=  expr2  ; 
exprl  >=  g  >=  expr2  ; 
abs(g)  <=  esrpr  ; 
expr >=  abs( g) ; 
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The  Constraints  Section/continued 


The  g  is  one  of  the  following  functionals 

•  re  (scLfctr ), 

•  im  (  scLfctr  *W  (i) ), 

•  mag  (scLfctr  *W  (z)  ) , 

•  null(0 ,  /), 

•  re  ( scLfctr *S  (  0 ,  /)  ) , 

•  im(scL/c£r*S(0,  /)), 

•  max_mag_S(0i,  j), 

•  max_mag_S(0i,  #2>  j)(/),  and 

•  max_magJ5(0i,  &2,  j)  C/i ,  /2,  fc). 
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Producing  and  reporting  awd  results 


Optimization  usually  ends  with 

•  found  strong  minimum:  the  unique  feasible 
optimal  weight  vector  (w)  was  found. 

•  found  weak  minimum:  non-unique  feasible 
weight  vector  (w)  that  minimizes  the  objec¬ 
tive  function  was  found. 

•  solution  appears  to 

be  unbounded:  happens  only  if  the  weight 
vector  (w)  in  the  AWDSL  program  does  not 
have  explicit  bounds.  It  is  generally  a  symp¬ 
tom  of  an  ill-formulated  design  problem. 

•  no  feasible  point  found:  there  is  no 
weight  vector  (w)  that  satisfy  the  con¬ 
straints.  Some  (or  all)  of  the  constraints  in 
the  constraints  section  must  be  relaxed. 
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Output  Listing 


The  listing  for  design  constraints  will  be  six- 
column  lines  containing  the  following  information 


constraint  value  lo_bnd  up_bnd  lambda  status . 

These  columns  are  as  follows. 

Column  1:  description  of  functional,  e.g. 

•  re(scLfctr*\Ki) ) ,  or  may  be 

•  maxjaag_S(0i,  02 ,  j)(/i,  /2,  k) . 
Column  2:  value  of  functional. 

Column  3:  value  of  lower  bound  on  functional. 
Column  4:  value  of  upper  bound  on  functional. 

Column  5:  value  of  Lagrange  multiplier  for  func¬ 
tional. 

Column  6:  This  field  is  empty  if  functional  is 
within  its  bounds.  Otherwise,  it  will  contain 
one  of: 

•  lb 

•  ub 

•  violates  lb 

•  violates  ub 
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Sample  Ouput  Listing 

The  following  segment  in  the  constraints  section  of 
the  source  file  filename 

re(S(0) )  >=  1.0; 
mag(S(0) )  <=  1.5; 
null(270,0 . 9) ; 

may  produce  the  following  lsiting 


re(S(0) ) 
mag(S(0) ) 
null (270, 0.9) 


1  1  le+10  1.04  lb 

1.082  0  1.5  0 

0  0  0  0 


Similarly, 

max_mag_S(90,270 ,50) (1 . 1)  <=  0.1; 
may  produce  the  following  lines  in  filename,  out 


max_mag_S (90, 270,50)  0.1035  0  0.1  1.87  ub 

Above  constraint  attains  bounds  at: 
1.1  90  0.1035  1.24 

1.1  180  0.0999  0.63 
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Abstract 

In  this  paper  we  introduce  a  new  computer-aided  approach  to  the  design  of  weight  vectors 
for  antenna  arrays  of  omnivariant  elements  responding  to  signals  in  the  far-ficld  of  the  array. 
The  problem  of  concern  is  the  deterministic  problem,  where  frequencies  and  directions  of  arrival 
for  sources  are  assumed  to  be  known. 

The  approach  relies  on  the  ability  to  cast  many  typical  design  problems  as  convex  program¬ 
ming  problems  with  the  real  and  imaginary  parts  of  the  complex  weight  vector  as  the  decision 
variables.  In  this  approach  we  introduce  an  antenna  weight  design  specification  language  com¬ 
piler  as  a  tool  for  translating  design  objectives  and  specifications  from  everyday  engineering 
language  into  a  convex  mathematical  program.  Resulting  convex  program  is  then  passed  to  a 
standard  convex  program  solver  and  results  arc  reported  in  an  easily  interpretable  form. 


1  Introduction 

1.1  Design  of  Directivity  Patterns 

The  problem  of  finding  array  clement  weights  so  as  to  shape  certain  directivity  patterns 
for  the  sensor  array  is  an  old  one  and  the  approaches  used  are  quite  diverse.  In  general,  we 
can  divide  the  methods  used  into  two  main  categories,  closed  form  and  iterative  techniques. 
On  the  one  hand,  closed  form  techniques  are  easy  to  utilize,  yet  the  class  of  problems  for 
which  such  methods  can  be  used  arc  usually  restricted.  Thus,  the  computational  efficiency 

'The  authors  are  currently  with  the  Information  Systems  Laboratory.  Department  of  Electrical  Engineer¬ 
ing,  Stanford  University,  Stanford,  CA  94305 
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is  offset  by  the  limited  applicability.  Oil  the  other  hand,  iterative  methods  usually  attain 
the  less  restricted  applicability  by  sacrificing  the  computational  case  offered  by  closed  form 
solutions. 

Historically,  dosed  form  methods  were  the  first  to  appear  in  literature.  The  ability  to  use 
tables,  charts,  and  hand  calculations  were,  in  our  opinion,  the  real  reason  for  this.  For  a  trans¬ 
mitting  array  of  equally  spaced  elements,  Stone1  used  a  binomial  expansion  representation  of 
the  array  pattern  with  the  exciting  currents  as  the  binomial  coefficients.  Schelkunoff  [Sch43] 
used  the  ability  to  represent  arrays  of  equally  spaced  antenna  elements  as  polynomials  to  de¬ 
cide  on  the  appropriate  currents  to  shape  directivity  patterns  of  transmitting  antenna  arrays 
by  intelligent  choice  of  zeroes  of  the  representing  polynomial.2  Dolph  [Dol46,  Rib-17,  Pri53] 
used  Tchebychev  polynomials  to  find  an  exceedingly  easy  solution  for  the  optimal  tradeoff 
between  sidelobe  gain  and  mainlobe  width.  Dolph  initially  assumed  a  linear  array  of  equally 
spaced  narrowband  transmitting  elements  symmetrically  placed  around  the  origin  so  that 
the  combined  array  output  can  be  assumed  totally  real,  and  that  the  exciting  currents  are 
in-phase.  Different  attempts  were  made  afterwards  to  relax  some  or  all  of  these  assumptions. 

The  significance  of  Dolph’s  work  is  that  he  was  the  first  to  talk  about  a  solution  to  the 
pattern  shaping  problem  that  is  optimal  in  some  sense.  This  in  itself  differentiates  his  work 
from  earlier  ones  which  gave  a  (and  not  the)  solution  to  the  problem.  What  followed  was  a 
search  for  solutions  that  are  optimal  in  one  sense  or  another.  For  instance,  one  of  the  methods 
proposed  was  that  of  representing  array  directivity  as  a  i-atio  of  two  “Hermitian”  forms 
and  attempting  to  maximize  this  ratio  by  use  of  generalized  eigenvalues  and  eigenvectors 
arguments  [CT65,  TC67].  Attempts  to  handle  extra  constraints  on  the  pattern  can  be  found 
in  [SB71,  VB72,  Che71].  The  last  two  papers  have  a  good  review  of  the  literature  (Western 
and  Eastern)  on  this  approach.  Such  Hermitian  ratios  are  not  convex,  and  thus  solutions 
obtained  for  the  constrained  case  are  generally  local.  That  is,  it  is  almost  always  the  case 
that  there  is  a  better  solution  that  the  method  did  not  locate. 

Waren  et.  al.  [LSW66,  WLS67],  suggested  the  use  of  mathematical  programming  for 
optimizing  different  aspects  of  antenna  design  (including  geometry).  The  problems  that  were 
posed,  whether  weight  design  or  otherwise,  were  generally  non-convex.  Lack  of  convexity 
meant  that  produced  solutions  were  typically  '‘local, ”  and  that  the  failure  to  produce  a 
set  of  parameters  that  satisfy  the  design  specifications  did  not  answer  the  question  about 
the  solvability  of  such  a  problem.  That  is,  the  method  failed  to  answer  the  question  ol 
whether  there  exists  any  set  of  parameters  for  which  the  design  specifications  can  be  satisfied. 
Furthermore,  convexity  is  an  underlying  assumption  for  nearly  all  efficient  optimization 
algorithms.  Thus,  specifying  non-convex  optimization  problems  limits  the  ability  to  utilize 
most  of  the  existing  efficient  optimization  procedures. 

Of  the  different  methods  utilizing  convex  programming  for  antenna  array  weight  design, 

‘U.  S.  Patents  1,643,323  and  1,715,433 

■’Using  the  principle  of  reciprocity  [Kra88],  exciting  currents  for  transmitting  antennas  can  he  thought  ol 
as  element  weights  for  receiving  antennas. 
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we  choose  to  mention  those  of  [MHM72,  EF/o,  Wii76,  \Yil7S].  Evans  and  Fortmann  [EF75] 
studied  the  related  problem  of  optimizing  patterns  for  antennas  with  continuous  line-source 
excitation.  A  functional  analysis  argument  was  used  to  solve  the  problem.  Although  this 
problem  is  removed  from  our  discussion,  the  link  can  be  seen  from  the  fact  that  an  array  can 
be  thought  of  as  a  “sampled  aperture”  (see  for  example  [Ste76] ) .  Of  particular  interest  to 
us  are  the  works  of  Wilson  [Wil76,  Wil7S]  and  McMahon,  et.  ai,  [MHM72].  In  particular, 
the  optimization  method  we  utilize  is  similar  to  that  proposed  in  [MHM72], 

In  [MHM72],  the  constraints  on  the  magnitude  of  the  pattern  were  translated  to  con¬ 
straints  on  the  real  and  imaginary  parts  of  the  pattern  function  (with  a  possible  3dI3  error). 
The  real  and  imaginary  parts  of  the  weights  were  taken  as  the  design  variables.  Thus,  con¬ 
straints  on  the  pattern  were  translated  to  linear  constraints  on  the  design  variables.  In  order 
to  obtain  a  linear  programming  solution  of  the  above  problem,  the  authors  “created”  a  linear 
objective  function  especially  tuned  so  as  to  reduce  the  mainlobe  beamwidth.  One  interesting 
aspect  of  the  above  approach  is  that  the  authors,  although  their  concluding  discussion  shows 
that  they  appreciated  the  significance  and  generality  of  their  approach,  in  our  opinion  failed 
recognize  a  major  advantage  that  their  approach  had  over  earlier  ones.  Since  the  problem 
solved  was,  due  to  discretization  of  the  pattern  and  the  3dB  errors,  less  restrictive  than  the 
original  problem  (where  the  magnitude  of  the  pattern  strictly  lies  within  the  designated  lim¬ 
its  for  all  angle  points),  and  since  the  original  problem  was  a  convex  one,  it  is  obvious  that 
if  the  method  discussed  above  failed  then  there  are  no  element  weights  that  will  satisfy  the 
original  constraints.  Thus,  if  the  designer,  using  such  a  method,  failed  in  obtaining  solutions 
satisfying  the  constraints,  the  problem  itself  as  posed  should  be  re-examined  rather  than 
attempting  different  approaches  to  solving  the  problem. 

Wilson,  [Wil76,  Wil7S],  utilized  a  similar  method.  In  his  approach,  Wilson  described  the 
array  pattern  as  a  function  of  the  phase  and  magnitude  of  element  weights.  By  imposing 
symmetry  of  the  array  around  the  origin,  the  phase  of  the  weights  was  assumed  to  be  constant 
and  the  pattern  was  assumed  to  be  totally  real  (not  unlike  Dolph’s  assumptions).  Thus, 
the  only  variables  remaining  were  magnitudes  of  the  weights.  An  approximation  scheme, 
[BY66],  was  then  used  to  find  the  variables  (magnitudes  of  weights)  that  will  minimize  the 
Tchebychev  (infinity)  error  between  the  actual  response  and  a  desired  response  at  selected 
angle  points.  This  approximation  scheme  is  based  on  linear  programming  techniques.  The 
interesting  note  about  this  method  is  that  had  Wilson  assumed  the  decision  variables  to 
be  the  real  and  ima0inary  parts  of  the  weights,  rather  than  the  magnitude  and  phase,  the 
symmetry  constraints  imposed  on  the  problem  could  have  been  easily  discarded  of,  and  the 
method  would  have  become  similar  to  that  of  [MIIM72]. 

1.2  A  Proposed  Approach  to  Optimal  Directive  Pattern  Shaping 

In  this  paper  we  propose  a  new  approach  to  optimal  directive  pattern  shaping.  We 
only  concern  ourselves  with  design  constraints  and  objectives  that  arc  proven  to  be  convex 
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functionals  of  the  weight  vector.  Non-convex  problems  will  not  be  considered.  This  is  not  to 
say  that  such  problems  are  not  important.  For  example,  the  problem  of  finding  the  optimal 
geometry  to  attain  a  certain  pattern  is  a  very  interesting  one.  Another  interesting  non- 
convex  design  problem  is  considered  in  [LJ8S] .  In  that  work  the  authors  study  the  problem 
of  minimizing  the  number  of  active  elements  (or  maximizing  the  number  of  zero  elements) 
in  the  antenna  array  while  satisfying  certain  design  constraints. 

Unfortunately,  the  current  status  of  mathematical  programming  does  not  allow  handling 
such  problems  efficiently.  Moreover,  solutions  produced  with  existing  algorithms  will,  at 
best,  be  local  solutions.  That  is,  there  is  no  way  of  knowing  how  far  such  a  solution  is 
from  the  global  optimum  one.  Furthermore,  failure  of  such  methods  to  produce  an  initial 
feasible  point  (if  the  constraints  that  define  the  set  of  feasible  designs  are  not  convex)  does 
not  convey  much  information  about  the  existence  of  a  feasible  point  (i.e.  one  which  satisfy 
design  constraints). 

In  addition,  experience  from  other  fields,  such  as  economics,  indicates  that  the  extra 
information,  e.g.  the  Lagrange  multipliers  or  shadow  prices  as  they  are  often  referred  to, 
available  from  a  typical  convex  program  solver,  is  in  many  instances  as  important  as  the 
actual  solution  itself.  These  shadow  prices  reflect  to  the  designer  the  tradeoff  between  the 
“tightness”  of  the  design  constraints  and  the  value  of  the  design  objective,  thus,  giving  the 
designer  an  opportunity  to  decide  whether  keeping  the  constraints  as  tight  as  they  arc  is 
“worth  it.” 

We  also  propose  the  use  of  a  specialized  design  language  compiler,  referred  to  as  the 
antenna  weight  design  specification  language,  as  a  mechanism  for  translating  the  design 
problem  from  engineering  specifications  to  the  corresponding  convex  mathematical  program. 
The  significance  of  this  approach  is  that  it  automates  the  important  step  of  translating 
antenna  weight  design  problems  into  convex  mathematical  programming  ones.  The  antenna 
designer  need  not  be  an  expert  in  mathematical  programming.  Thus,  the  designer  should  not 
carry  the  burden  of  translating  the  engineering  specifications,  which  may  resemble  a  large 
number  of  objectives  and  constraints,  into  a  convex  mathematical  program.  The  designer 
should  also  expect  results  of  the  design  process  to  “make  sense”  from  an  engineering  point 
of  view.  Although  mathematical  programming  results  from  a  typical  convex  program  solver 
contain  enormous  information,  this  information  need  not  make  sense  to  a  mathematical 
programming  novice. 

This  specification  language  approach,  combined  with  convex  programming,  has  been 
successfully  used  recently  as  a  computer  aided  design  tool  in  designing  controllers  [BBB*SS, 
Khr90]  for  linear  time-invariant  dynamic  systems.  In  that  work,  the  ability  to  parameterize 
all  stabilizing  controllers  of  a  dynamic  system  combined  with  the  lad  that  many  specifica¬ 
tions  on  the  control  system  translate  into  convex  constraints  on  these  free  parameters,  were 
exploited  to  represent  the  design  problem  as  a  convex  program. 
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2  The  Antenna  Weight  Design  Problem 

2.1  Underlying  Assumptions 

Consider  a  planar  array  of  N  sensor  elements  with  fixed  relative  locations.  Elements  are 
assumed  to  have  negligible  gain  outside  a  predetermined  (finite)  frequency  band  of  operation. 
This  can  be  attained  in  practice  by  passing  the  sensor  output  through  an  “ideal"  band -pass 
filter  (BPF)  before  processing  the  signal.  Medium  of  transmission  is  assumed  to  be  isotropic 
and  non-dispersive.  Signal  sources  are  assumed  to  be  farfield  and  lie  in  the  plane  of  the 
array  as  shown  in  Figure  1.  Furthermore,  it  is  assumed  that  each  sensor  in  the  array  will 
have  an  additive  zero-mean  stationary  noise  component  at  its  output.  This  noise  can  be  due 
to  thermal,  amplifier,  and/or  spatial  noise.  This  noise  is  assumed  to  be  uncorrelated  with 
the  signal  component.  Without  loss  of  generality,  it  is  also  assumed  that  noise  components 
in  different  elements  are  uncorrelated.3 


y 


Letting  s,(0,u;)  be  the  output  of  the  ilh  element  with  respect  to  a  certain  signal  arriv¬ 
ing  from  the  0  direction  with  angular  frequency  u>,  a  fundamental  assumption  is  that  the 

JA  simple  whitening  argument  can  be  used  to  justify  this  for  correlated  noise 


VIII- 


5 


(1) 


combined  array  output  due  to  such  signal  can  be  written  as 

,v 

y{0,u)  =  DtDi(*(0fW)  +  7i,). 

t  =  l 

where  w ,  is  referred  to  as  the  element  weight  of  the  ith  element,  w,  is  its  complex  conjugate, 
and  77,  is  the  corresponding  additive  noise  for  such  an  element.  This  is  equivalent  to  saying 
that  element  weights  afFect  the  combined  array  output  in  a  linear  fashion.  Letting  tv,  s(0,tv), 
and  7i  be  complex  vectors  in  CN ,  with  the  itk  element  of  these  vectors  as  u>,-,  s,(0,u),  and  n,, 
respectively,  then  another  way  to  state  the  combined  array  response  is 

y(0,  w)  =  wm(s(0,<v)  +  n)  (2) 

where  the  operation  represents  taking  the  conjugate  transpose  (Hermitian)  of  the  vector 
tv.  1  his  tv  is  usually  referred  to  as  the  weight  vector.  Sometimes  we  shall  refer  to  the 
combined  array  noise-free  response  as 

S(0,iv)  =  U)’s(0,u)).  (.'}) 

As  for  the  noise  component,  under  the  assumptions  above,  the  total  noise  power  i.i  the 

output  can  be  thought  of  as 

K  =  J2at\wi\2'  (4) 

i=l 

where, 

£  =  diag(<7x,  cr2, . . .  ,cr.v)  =  £(nn*).  (5) 

i.e.  the  cr’s  are  non-negative  real  numbers. 

The  relationship  between  sensor  element  attributes  and  the  combined  array  output  is 
shown  in  Figure  2  for  a  linear  array  of  N  elements.  In  this  figure  we  have,  on  purpose,  dis¬ 
tinguished  between  the  noise-free  response  of  each  sensor  element,  s,  in  the  above  discussion, 
and  the  additive  noise,  nt. 

2.2  Problem  Statement 

Assume  that  we  have  an  antenna  array  as  described  above  and  that  the  array  is  receiving 
signals  with  potential  interference  or  jamming  Using  adaptive  antenna  design  techniques, 
a  certain  criterion  ab<  ut  the  quality  of  reception  of  a  desired  signal  versus  reception  of  the 
jamming,  such  as  SNR,  can  be  optimized.  Most  such  methods  read  to  jamming  by  placing 
a  pattern  null  (that  is  a  zero  of  the  directive  pattern  function)  at  the  directions  in  which 
jamming  is  expected.  Although  these  solutions  are  optimal  as  far  as  the  optimization  crite¬ 
rion  is  concerned,  they  are  quite  sensitive  to  changes  ii.  incident  angles  and  other  deviations 
from  underlying  assumption ;. 
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For  the  deterministic  weight  design  problem,  it  is  assumed  that  desired  signals  are  arriving 
at  known  angles  of  incidence  with  prespecified  frequencies.  It  is  also  assumed  that  jamming 
signals  are  arriving  at  predetermined  angles  and  frequencies.  The  weight  design  problem  can 
be  stated  as  follows. 

Find  an  antenna  array  (complex)  weight  vector  that  extremizes  a  given  design 
objective  while  satisfying  certain  design  constraints  on  the  directive  pattern  func¬ 
tion  of  the  array. 

3  The  Convex  Mathematical  Programming  Approach 

3.1  Basic  Definitions 

In  this  section,  some  of  the  basic  definitions  and  terminology  used  in  the  remainder  ol 
this  paper  will  be  reviewed.  For  a  detailed  discussion,  the  reader  is  advised  to  refer  to  any 
standard  textbook  on  the  subject  (e.g.  Luenberger  [Lue69,  LueS4],  Bazaraa  and  Shetty 
[BS79],  or  Rockafeller  [Roc70]). 

In  the  ongoing  argument,  V  is  assumed  to  be  any  linear  space  on  the  real  field  3?. 

Definition  1  (Convex  Set)  ft  C  V  is  said  to  be  a  convex  set  if  /tx  +  (l  —  /i)y  6  ft  whenever 
x,  y  6  ft  and  0  <  y.  <  1. 

This  is  to  say  that  if  two  points  are  in  a  convex  set,  then  any  point  between  them  on  the 
line  connecting  the  two  points  is  also  in  the  set.  It  is  usually  assumed  that  the  empty  set  is 
a  convex  set. 

Definition  2  (Convex  Function)  A  function  f  :  V  —>  is  a  convex  function  on  a  convex 
set  ft  C  V  if  f  (fix  +  (1  —  y)y)  <  yf(x)  +  (1  —  y)f{y)  for  all  x,y  6  ft  and  0  <  /t  <  1. 

This  is  to  say  that  the  line  connecting  two  points  on  the  graph  of  a  convex  function  lies 
above  the  graph  of  the  function  between  these  two  points. 

A  convex  mathematical  program  is  the  problem  of  finding  an  element  in  a  given  convex 
set  that  minimizes  the  value  of  a  certain  convex  function  on  the  set.  This  can  be  stated  as 

minimize  $(u;)  ^ 

subject  to  uj  6  ft 

where  $  :  V  — ♦  3?  is  a  convex  function  on  the  convex  set  ft  C  \  .  This  convex  function,  <!>,  is 
usually  referred  to  as  the  objective  function  while  the  set  ft  is  often  referred  to  as  the  feasible 
set.  ft  is  usually  assumed  to  be  a  closed  set. 

Many  of  the  results  on  convexity  will  not  be  mentioned  here,  but  it  should  be  noted  that 
an  intersection  of  convex  sets  is  convex,  and  that  any  positively  weighted  linear  combination 
of  convex  functions  is  a  convex  function. 
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3.2  Why  Convex  Programming 

A  general  mathematical  programming  problem  is  the  same  as  that  described  in  6,  except 
for  (one  of)  the  convexity  requirements  on  the  objective  function  and  the  set.  Although  this 
provides  a  more  general  setup  that  allows  formulating  much  more  general  design  problems, 
such  as  those  mentioned  in  the  introduction,  yet,  it  is  desired  to  maintain  convexity  in  the 
mathematical  programming  formulation  of  the  problem  for  the  reasons  discussed  below. 

The  most  important  feature  that  a  convex  program  has  is  that  local  optimal  solutions 
are  global  optimal  solutions.  That  is,  no  local  solutions  exist.  Thus,  the  notion  of  global 
optimality  for  a  convex  program  is  well  defined.  That  is,  it  can  be  easily  verified  whether  a 
solution  is  the  optimal  one  or  not  by  using  some  of  the  well-known  local  optimality  tests, 
such  as  the  Kuhn-Tucker  conditions  [Lue84].  This  is  not  true  for  a  general  mathematical 
program  since,  generally  speaking,  there  is  no  quick  test  for  global  optimality  of  a  solution 
since  solutions  for  a  general  mathematical  programming  problem  can  be  local.  Furthermore, 
local  optimal  solutions  for  a  general  non-convex  program  may  be  arbitrarily  far  from  the 
global  optimal  one.  Thus,  it  may  very  well  be  the  case  that  starting  the  same  algorithm 
from  different  initial  points  will  produce  different  results. 

Another  advantage  that  convex  programming  has  over  general  mathematical  program¬ 
ming  is  the  abundance  of  convex  programming  algorithms  that  are  effective  and  robust.  By 
effectiveness  it  is  meant  that  given  a  convex  programming  problem  of  a  certain  type,  an 
algorithm  can  be  found  that  will  correctly  determine  whether  there  are  any  feasible  points. 
If  there  are  any,  the  optimal  one  will  be  found.  On  the  other  hand,  robustness  means  that 
an  algorithm  is  numerically  stable  and  will  handle  a  wide  spectrum  of  problems  of  a  certain 
type  ( e.g .  linear  programs)  while  producing  results  that  are  correct  within  a  prespecified 
numerical  tolerance.  Furthermore,  for  most  classes  of  convex  programming  problems,  algo¬ 
rithms  have  been  developed  with  well  estimated  rates  of  convergence.  Thus,  such  algorithms 
are  guaranteed  to  produce  results  within  a  prespecified  precision  in  a  precalculable  amount 
of  time. 

One  more  useful  feature  of  convex  programs  is  that  solution  sets  are  convex.  Thus,  even 
for  non-differentiable  convex  optimization  problems  [Kel60,  Sho85,  AkgS4],  there  are  many 
methods  that  will  try  to  find  the  solution  set  by  continuously  refining  a  convex  set  that 
contains  the  solution  set.  Many  of  these  methods  have  well  established  (although  generally 
poor)  rates  of  convergence.  Many  of  these  methods  are  still  undergoing  intensive  research 
which  may  have  a  promise  for  the  future. 

Finally,  most  algorithms  which  solve  convex  programs  will  produce  what  is  known  as 
the  Lagrange  multipliers  or  shadow  prices  which  resemble  a  solution  for  a  dual  optimiza¬ 
tion  problem.  Although  the  main  use  of  these  multipliers  within  an  algorithm  is  to  verify 
the  optimality  of  the  solution  found,  they  have  another  interesting  practical  interpretation. 
These  multipliers  reflect  the  inherent  tradeoff  among  objectives  and  constraints.  That  is, 
the  multipliers  (shadow  prices)  tend  to  quantify  the  answer  to  the  question 
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“how  much  does  maintaining  a  constraint  as  tight  as  it  currently  is  contribute 
to  the  value  of  the  objective  function?” 

A  large  positive  multiplier  indicates  that  relaxing  the  lower  bound  on  the  related  constraint 
may  significantly  reduce  the  value  of  the  objective  function,  while  a  large  negative  multiplier 
indicates  that  the  upper  bound  should  be  examined. 

This  is  quite  an  important  feature  since  in  most  engineering  applications  and  otherwise, 
the  designer  does  not  usually  have  an  a  priori  knowledge  of  the  exact  tradeoffs  among  certain 
design  constraints  and  design  objectives.  Quantifying  these  tradeoffs  is  a  valuable  asset  to 
the  designer. 

3.3  The  Antenna  Weight  Design  Problem  as  a  Convex  Program 

The  ability  to  write  a  weight  design  problem  as  a  convex  mathematical  program  is  now 
investigated.  In  particular,  we  note  that  there  is  a  one-to-one  correspondence  between  vectors 
of  CN  and  vectors  in  ?R2N .  For  a  vector  w  £  CN  can  be  thought  of  as  w#  +  jw$,  where 
w%,  €  3?^  and  j  =  \/—l.  Similarly,  we  can  think  of  a  vector  W  £  3?27V  as  WT  =  [u;^,  w^\T ■ 
That  is,  the  first  N  elements  of  a  vector  in  3£2A/  can  be  thought  of  as  corresponding  to  the 
real  part  of  a  complex  vector  in  CN ,  while  the  latter  N  elements  can  be  thought  of  as 
corresponding  to  the  imaginary  part  of  the  complex  vector.  This  remark  permits  the  study 
of  the  weight  design  problem  as  a  problem  on  3?2JV,  where  N  is  the  number  of  antenna 
elements,  rather  than  a  problem  on  CN . 

In  the  sequel,  it  will  be  assumed  that  the  weight  vector  w  £  CN  is  written  as  w  =  w^+jiUi j, 
where  w#,  £  3?w  are  called  the  real  and  imaginary  components  of  the  weight  vector. 

Instead  of  explicitly  considering  iyr  =  [iugj,  w$]T  as  a  vector  in  3?2Af,  the  pair  and  w 3 
will  be  considered. 

The  relation  between  w^,w<3  and  W  is  a  simple  linear  transformation  of  the  form 

w#  =  [/yvlOivJfF,  and 
^  =  [On\In]W, 

where  I yv  is  the  N  x  N  identity  matrix,  while  0^  is  the  N  x  N  zero  matrix.  Thus,  a  linear 
function  on  u>»  and  uk,  corresponds  to  a  linear  function  on  W ,  and  a  convex  function  on 
and  iuq  corresponds  to  a  convex  function  on  W . 

With  the  above  remarks  in  mind,  some  convex  functions  on  the  real  and  imaginary  parts 
of  the  weight  vector  are  discussed  next.  Such  functions  may  appear  as  an  objective  or  be 
used  to  construct  constraints  by  setting  upper  and/or  lower  bounds.  Here,  it  should  be  noted 
that,  to  maintain  convexity  of  the  set  satisfying  the  constraints,  lower  bounds  can  be  set 
only  for  linear  functionals. 

•  Total  noise  power,  Pn  =  cr,|iy,j2  =  ru'S w,  can  be  rewritten  Pn  — 
or  Pn  =  ai{wx  +  This  is  obviously  a  convex  function  since  the  individual  ct’s 
are  non-negative. 
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•  The  real  and  imaginary  parts  of  a  certain  element  weight,  wt,  are  (obviously)  linear 
functions  of  the  real  and  imaginary  parts  of  the  weight  vector. 

•  The  magnitude  of  a  certain  element  weight,  re,-,  is  (obviously)  convex  in  the  real  and 
imaginary  parts  of  the  weight  vector. 

•  Any  norm  of  the  weight  vector,  is  (obviously)  convex  in  the  real  and  imaginary  parts 
of  the  weight  vector. 

•  The  real  and  imaginary  parts  of  the  noise-free  pattern  function,  5(0 ,  /),  are  linear  in  the 
real  and  imaginary  parts  of  the  weight  vector.  To  see  this,  note  that  5(0,  /)  =  w  s(0,  /); 
letting  s(0, /)  =  •Sj{(0,/)  ■+■  js§(0,  /),  it  becomes  obvious  that  the  real  part  of  5  is 
w%S3i  —  WqSq,  while  the  imaginary  part  of  5  is  +  w%s%;,  both  of  which  are  linear. 

•  The  array  noise-free  pattern  gain,  |5(0,/)|,  is  convex  in  the  real  and  imaginary  parts 
of  the  weight  vector.  This  can  be  easily  seen  from  the  fact  (the  previous  point)  that 
the  real  and  imaginary  parts  of  5(0,/)  are  linear  in  the  real  and  imaginary  parts  of 
the  weight  vector. 

•  Also,  sup  |5(0,/)|  for  any  fixed  0,  sup  |5(0, /)|  for  any  fixed /,  and  sup  |.c>(0. 

/,</</2  ®i<0<02 

are  also  convex  for  any  0i,  02,  /i,  and  /2. 

•  Finally,  any  positively  weighted  linear  combination  of  the  above  functions  is  also  con¬ 
vex. 

4  The  Specification  Language  Approach 

It  was  shown  above  that  many  typical  antenna  weight  design  problems  for  a  large  class 
of  antennas  can  be  posed  as  convex  mathematical  programs.  The  theoretical  value  of  this 
approach  was  discussed.  The  practicality  of  such  a  design  method  remains  limited,  as  far 
as  the  antenna  weight  designer  is  concerned,  unless  such  a  designer  can  specify  the  actual 
design  problem  as  a  convex  program  and  interpret  the  results  efficiently.  A  computer  aided 
design  approach  to  translating  the  antenna  weight  design  problem  discussed  above  into  a 
mathematical  program  is  a  most  natural  way  for  solving  this  dilemma. 

To  do  this,  we  propose  the  use  of  an  antenna  weight  design  specification  language  (A  WDSL). 
The  purpose  of  such  a  language  is  to  enable  the  designer  to  fully  specify  the  design  prob¬ 
lem  in  everyday  engineering  language  and  not  as  a  convex  mathematical  program.  Thus, 
the  designer  can  set  up  the  design  problem  in  a  natural  way  by,  say,  setting  pattern  nulls, 
specifying  limits  on  the  magnitude  of  the  array  response,  or  selecting  any  ol  the  convex 
functionals  mentioned  above  as  a  design  objective.  A  specification  language  compiler  then 
translates  the  design  problem  into  a  convex  mathematical  program.  1  lie  designer  need  not 
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know  the  resulting  mathematical  program.  This  is  not  unlike  the  fact  that  a  good  computer 
programmer  need  not  know  anything  about  the  machine  language  for  the  computer  that  is 
being  used  without  sacrificing  the  ability  to  produce  correct  and  efficient  computer  programs. 
This  convex  program  can  then  be  solved  using  any  convex  program  solver.  Afterwards,  the 
results  are  reported  to  the  designer  in  a  form  similar  to  what  was  used  to  specify  the  prob¬ 
lem.  This  general  idea  is  outlined  in  Figure  3.  Using  the  above  outline,  a  description  of  the 


Figure  3:  Outline  of  how  to  implement  a  specification  language  approach. 

specific  approach  utilized  in  this  paper,  called  awd  [KluSS],  is  illustrated  in  Figure  4 

Figure  5,  shows  how  this  specification  language  approach  fits  into  the  typical  engineering 
design  cycle.  The  designer  starts  to  design  element  weights  for  a  sensor  array  so  as  to 
satisfy  certain  engineering  design  specifications.  The  designer  may  not  know  a  priori  whether 
the  design  constraints  are  attainable,  or  it  may  be  the  case  that  the  current  specifications 
are  based  on  earlier  iterations  (the  dashed  line  in  the  Figure).  Using  the  array  and  design 
specifications,  the  designer  sets  up  an  AWDSL  program  and  inputs  it  to  the  language  compiler 
resulting  in  a  convex  program.  Additional  information  about  the  problem  may  be  available 
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too.  The  convex  program  is  then  passed  to  a  convex  program  solver.  The  designed  weights 
(u)opt)  and  certain  additional  information  (e.g.  Lagrange  multipliers)  will  be  available  at  the 
end  of  this  phase.  If  the  problem  is  feasible,  simulating  the  results  will  show  whether  they 
are  implementable  or  not.  Results  of  the  simulation,  combined  with  the  information  about 
problem  feasibility  and  the  Lagrange  multipliers,  may  serve  as  a  tool  for  refining  the  design 
and/or  array  specifications. 

Infeasibility  of  the  problem  (under  convexity)  gives  a  clear  answer  to  whether  the  design 
specifications  can  be  attained  or  not  by  any  set  of  weights.  If  the  problem  is  infeasible, 
the  answer  is  an  unqualified  “no”  as  was  discussed  earlier.  This  in  itself  is  a  valuable 
product  of  the  method.  For  if  the  problem  is  infeasible,  the  designer  should  decide,  from 
the  additional  information  provided  by  the  convex  program  solver,  which  design  constraints 
were  the  infeasible  ones.  Then  these  “violated”  constraints  should  either  be  “relaxed,”  if  this 
is  acceptable,  or  a  different  array  be  used  if  it  is  vital  to  maintain  the  original  constraints. 
By  a  different  array  we  mean  that  the  geometry  and/or  gain  patterns  of  individual  sensors 
are  changed.  It  is  up  to  the  designer  to  decide  which  action  to  take. 

On  the  other  hand,  even  with  a  feasible  problem,  the  Lagrange  multipliers  reveal  to  the 
designer  the  tradeoff  between  attaining  certain  design  constraints  and  the  increase  of  the 
objective  value.  It  is  up  to  the  designer  to  decide  whether  maintaining  constraints  as  tight 
as  they  are  is  worth  the  corresponding  increase  in  the  objective  value.  If  the  answer  is  no, 
the  constraints  that  “cost”  most  to  maintain  should  be  relaxed. 

From  the  above  discussion,  it  can  be  seen  that  the  specification  language  approach  re¬ 
sembles  a  significant  step  towards  making  the  convex  programming  approach  to  antenna 
weight  design  a  practical  engineering  tool.  Using  this  approach,  posing  the  problem  and 
interpreting  the  results  is  a  task  manageable  by  the  typical  antenna  designer. 

Of  the  many  possible  structures  of  an  antenna  weight  design  specification  language 
(AWDSL),  we  propose  the  following.  The  AWDSL  program  should  contain  three  major 
sections, 

•  array  description, 

•  objective  function ,  and 

•  design  constraints. 

The  array  description  section  is  to  provide  all  the  information  necessary  to  evaluate  the 
objective  and  constraint  functionals.  For  instance,  the  number  of  array  elements  and  their 
geometry,  operational  frequency  band,  signal  and  noise  covariance  matrices,  and  element 
gain  patterns  should  be  stated  in  this  section  explicitly  or  by  giving  names  of  files  that 
contain  such  information. 

The  objective  section  should  contain  design  objectives.  These  should  be  convex  lunction- 
als  to  be  minimized.  Positively  weighted  linear  combinations  of  convex  functionals  and  the 


VIII- 


M 


Solution  wopl 


Additional 

Information 


(Feasibility, 

Lagrange  Multipliers,  etc...) 


I 

Simulation 


_ J 


Figure  5:  Overall  structure  of  design  procedure. 


maximum  possible  value  of  a  set  of  functionals  are  themselves  convex,  as  was  stated  above, 
and  thus  may  be  allowed. 

Finally,  the  constraints  section  should  contain  design  constraints.  Lower  and  upper 
bounds  on  affine  functionals  should  be  allowed.  Non-afiine  convex  functionals  should  have 
upper  bounds  only. 

In  the  following  we  describe  our  implementation  of  awd  [KhrSS],  a  compiler  for  a  small 
but  effective  subset  of  AWDSL. 


5  An  Implementation  of  a  Specification  Language 

5.1  Issues  Concerning  Implementation 

First  of  all  we  discuss  some  of  the  issues  involved  in  choosing  and  implementing  the 
AWDSL  subset.  In  particular,  the  following  closely  related  issues  are  involved  in  such  an 
implementation. 

5.1.1  Underlying  assumptions 

In  order  to  simplify  the  problem  from  the  computational  point  of  view,  the  following 
assumptions  are  made  in  conjunction  with  those  specified  earlier.  The  compiler  is  to  handle 
only  planar  arrays  of  omnidirectional  elements  with  fixed  relative  locations  and  negligible 
mutual  interactions.  Elements  are  assumed  to  have  constant  complex  gain  over  a  predeter¬ 
mined  (finite)  frequency  band  of  operation.  Signal  sources  are  assumed  to  be  farfield  and  lie 
in  the  plane  of  the  array. 

It  was  decided  that  the  finite  operational  frequency  band  of  interest  be  selected  by  setting 
lower  (/lower)  and  upper  (/upPe r)  frequencies  in  the  array  description  section.  Antenna  element 
responses  are  assumed  to  be  constant  over  the  selected  band.  Thus,  there  is  no  need  for 
storing  or  computing  element  patterns.  It  was  also  decided  that  all  frequencies  be  specified 
in  units  of  a  certain  nominal  frequency,  /nominal-  Although  choice  of  /nominal  is  arbitrary, 
a  good  choice  is  the  center  frequency  of  the  antenna  element  operational  frequency  band. 
A  frequency  of  interest  (/)  in  physical  units  (Hertz)  is  input  as  a  normalized  frequency, 
/  =  ///nominal-  In  other  words,  a  frequency  of  0.5  specifies  a  frequency  of  half  the  nominal 
frequency. 

As  usual,  locations  (in  rectangular  coordinates)  of  array  elements  must  be  given  in  units 
of  nominal i  the  wavelength  corresponding  to  the  nominal  operation  frequency  (/nominal)  aiul 
the  underlying  speed  of  propagation  (c)  in  the  medium  (Allolninai  =  c//n„,„inai)- 

It  is  also  assumed  that  the  noise-free  directive  pattern  function  of  concern  to  the  designer 
is  the  array  space  factor  [RWDS4,  Sch43j.  This  is  a  function  of  Angle  of  Incidence,  0 , 
(normalized)  signal  frequency,  /,  complex  weights,  w,' s,  and  element  locations.  Under  the 
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above  assumptions,  this  function  is  given  by  equation  3,  where  in  this  case 


s,(0,  /)  —  exp(j’27r  f(x{  cos  0  +  yx  sin  0 )) 


(') 


where 

/  is  the  normalized  (dimensionless)  signal  frequency, 

0  is  the  Angle  of  Incidence,  and 

Xi  and  i /,  are  the  x  and  ^-locations  (in  units  of  A„on,jna|)  of  element  i. 

5.1.2  Language  capabilities 

Not  all  convex  functionals  of  the  real  and  imaginary  parts  of  the  weight  vector  can  be 
incorporated  within  the  language  without  paying  the  price  from  a  computational  point  of 
view.  This  has  to  do  with  the  fact  that  most  efficient  optimization  codes  exploit  a  certain 
problem  structure,  such  as  linearity  or  differentiability.  This,  unfortunately,  limits  ability  to 
handle  all  types  of  convex  objectives  and  constraints  discussed  above.  Our  implementation 
handles  design  problems  that  can  effectively  be  approximated  by  linear  or  quadratic  pro¬ 
grams.  That  is,  linear  constraints  with  linear  or  quadratic  objectives.  This  choice  has  to  do 
with  the  fact  that  linear  and  quadratic  convex  programming  are  intensively  researched  fields 
with  good  collection  of  robust  and  effective  algorithms. 

The  following  objective  functionals  fit  the  above  description;  that  is,  they  can  be  effec¬ 
tively  approximated  by  linear  or  quadratic  functionals  (as  will  be  discussed  in  a  following 
section). 

•  max  max  |5(0,/)|,  maximum  value  that  the  magnitude  of  the  array  space  factor 
!\<J <h  S\<.6<62 

assumes  between  angles  0X  and  02  for  any  (normalized)  frequency  between  f\  and  _/2. 

•  total  noise  power  in  the  combined  array  output. 

The  same  can  be  said  about  the  following  functionals  that  can  be  used  to  construct  design 
constraints. 

•  max  max  15(0,  f)\,  same  as  above. 

Jl<f<h^<0<02 

•  real,  imaginary,  or  magnitude  of  a5(0,/),  where  a  is  some  complex  weighting  factor. 
0  is  some  desired  Angle  of  Incidence,  and  /  is  a  desired  (normalized)  frequency. 

•  real,  imaginary,  or  magnitude  of  aw,,  where  a  is  some  complex  weighting  factor,  and 
iw,  is  the  complex  weight  factor  for  element  i. 

Another  design  constraint  can  be  imposed  by  setting  S(0,f)  to  zero  at  a  given  Angle  of 
Incidence,  0,  for  some  (normalized)  frequency,  /.  This  is  termed  placing  a  null  at  0  and  J . 
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5.1.3  Approximations 

A  magnitude  constraint,  such  as  |i<j,|  <  b  or  (.9(0,  /)(  <  b,  can  be  viewed  as  a  disc  in 
the  complex  plane  of  radius  b.  This  disc  can  be  arbitrarily  approximated  by  a  polygon. 
In  particular,  such  a  disc  may  be  approximated  by  an  octagon  to  within  ±0.34dl3  error. 
This  can  be  accomplished,  say  for  |S|  <  b  and  S  =  +  jSa,  by  replacing  the  magnitude 

constraint  by  the  following  four  inequalities 
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As  for  sector  conditions,  such  as  max|Sj  over  a  certain  region  in  AOI  or  frequency, 
the  region  will  be  sampled  first,  then  the  magnitude  functional  at  each  sampling  point  will 
be  approximated  as  above.  The  resolution  of  this  sampling  procedure  can  be  left  to  the 
program  or  be  entered  by  the  designer  as  an  optional  argument.  Thus,  a  constraint  of  the 
form  max|S|  <  b  over  a  certain  region,  is  approximated  as  |Sj  <  b  for  a  set  of  points 
uniformly  spaced  over  the  region.  Each  of  these  |5|  <  b  constraints  will  be  approximated  by 
four  inequalities  as  above. 

On  the  other  hand,  max  |5|  appearing  as  an  objective  in  the  form 

minimize  max  |5‘| 


over  a  certain  region  in  AOI  and/or  frequency,  will  be  handled  by  introducing  a  new  variable, 
say  z,  and  replacing  the  problem  by  the  equivalent  problem 

minimize  ~ 

subject  to  max|S|  <  s. 

The  latter  constraint  can  be  handled  as  above. 

5.1.4  Implementation  of  language  compiler 

The  question  of  how  to  implement  the  language  compiler,  from  the  programming  point  ol 
view,  is  a  fundamental  one.  In  our  implementation  of  awd,  standard  UNIX,  [K PS-1,  SF85], 
utilities  were  used  to  facilitate  this  formidable  task.  In  particular,  lex  [LSS-l],  a  lexical 
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analyzer,  and  yacc,  [Joh84],  a  compiler  compiler,  were  used  to  construct  the  compiler  part 
of  awd. 

As  shown  in  Figure  6,  a  lexical  analyzer,  built  using  lex,  will  detect  language  tokens  and 
pass  them  to  a  parser  (compiler),  constructed  using  yacc.  This  parser  recognizes  certain 
sequences  of  tokens.  If  a  sequence  is  complete,  a  certain  action  is  carried  out.  These  actions 
include  creating  and  updating  a  certain  internal  structure  called  the  design  problem.  This 
design  problem  is  actually  independent  of  the  type  of  optimization  code  used.  Syntax  errors 
are  detected  at  this  phase. 

As  a  second  phase,  the  design  problem  is  translated,  using  a  translator  written  in  C 
language,  into  a  convex  quadratic  or  linear  program,  as  appropriate.  This  convex  program, 
as  expected,  is  fine  tuned  to  the  specific  convex  solver  used.  This  is  why  we  have  included 
the  convex  solver  as  an  extra  third  phase  within  awd  (refer  to  Figure  6).  It  should  be  noted 
that  this  last  phase  is  not  a  part  of  the  language.  Dependence  of  the  compiler  output  (the 
convex  program)  on  the  convex  solver  makes  incorporating  the  solver  with  the  compiler  a 
sensible  thing  to  do. 

We  have  also  made  use  of  the  C-preprocessor  (cpp)  [KR7S],  which  is  a  UNIX  utility 
that,  among  other  things,  allows  defining  constants  using  a  “^define”  facility,  including 
external  files  within  the  source  file  using  an  “^include”  facility,  and  permits  C  language' 
style  comments  (he.  ignore  text  within  "/*"  and  "*/”).  The  source  program  is  passed  to 
cpp  before  passing  it  to  lex. 

5.1.5  Choice  of  a  convex  program  solver 

LSSOL,  [GHM*S6],  a  convex  program  solver  for  linear,  quadratic  and  least  squares  opti¬ 
mization,  from  the  Systems  Optimization  Laboratory  at  Stanford  University,  was  chosen  as 
the  convex  programming  code.  The  code  is  quite  staole  from  the  numerical  point  of  view, 
and  is  quite  efficient  on  problems  of  the  type  encountered  in  the  case  at  hand.  It  is  the 
case,  as  was  discussed  above,  that  most  of  the  convex  functionals  encountered  in  the  current 
implementation  of  awd  are  linear,  quadratic,  or  can  be  arbitrarily  (as  was  shown  earlier) 
approximated  by  a  set  of  quadratic  or  linear  functionals. 

Furthermore,  LSSOL  utilizes  the  active  set  method  [GMWS1].  Thus,  LSSOL’s  perfor¬ 
mance  is  enhanced  with  the  addition  of  equality  constraint.,.  Such  equality  constraints  (e.g. 
pattern  nulls)  are  quite  common  in  antenna  weight  design  practice.  In  addition,  LSSOL 
does  not  exploit  sparsity  of  matrices.  This  is  actually  an  advantage  since  the  mathematical 
programs  that  result  from  translating  the  design  problem,  typically,  have  full  matrices. 

5.2  Description  of  Awd 

Next,  we  give  a  brief  description  of  the  implemented  subset  itself.  In  par* icular,  syntax 
is  similar  to  the  C  language.  Statements  are  separated  by  semi-colons  (';’)-  and  groups  ol 
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statements  (multiple  statements)  are  enclosed  by  braces  (‘{’,  The  input  is  free  format. 
That  is,  any  number  of  spaces,  newlines,  and  tabs  can  replace  a  space. 

Looping  is  allowed  and  takes  the  following  form, 

for  var  =  exprl  to  cxpr2  do  strut 

for  var  =  exprl  to  exprS  step  cxpr3  do  stmt 


Numerical  constants  are  the  same  as  in  C.  Certain  constants  such  as  7r  (T 1-1...)  ami  c  (2. 7 IS...) 
are  defined.  Variables  (which  can  take  on  only  scalar  floating-point  values)  and  variable 
names  are  allowed  to  assume  single  lower  case  letters  only.  Furthermore,  most  mathematical 
functions  (cos,  sin,  etc.)  are  defined. 

Next,  we  shall  discuss  how  each  of  the  three  major  sections  of  an  AWDSL,  array  descrip¬ 
tion,  constraint,  and  objective  sections,  was  tailored  to  fit  our  implementation  of  awd.  For 
a  sample  awd  program,  the  reader  is  advised  to  refer  to  Example  6.1. 

5.2.1  The  array  description  section 

This  section  contains  the  basic  information  about  the  array  structure  and  element  behav¬ 
ior  that  allows  us  to  compute  any  permissible  functional.  In  view  of  the  basic  assumptions 
underlying  awd,  and  capabilities  of  the  implemented  language,  the  following  information  is 
enough. 

•  The  number  of  sensor  elements  in  the  array. 

•  The  upper  and  lower  bounds  on  the  normalized  operational  frequency  band. 

•  The  x  and  y  locations  (in  rectangular  coordinates)  of  the  array  elements  in  units  of 
nominal  wavelength,  Anominai. 

•  The  noise  power  of  each  element. 

Optionally,  an  initial  estimate  of  the  solution  (the  weights)  may  be  given  in  this  section. 
The  default  starting  weights  are  zeroes. 

5.2.2  The  objective  section 

The  objective  function  is  formed  by  summing  all  objective  terms  appearing  in  the  objective 
section.  An  objective  term  is  of  the  form: 

ex pr  *  g; 
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The  optional  (non-negative)  scalar  expression  expr  allows  relative  weighting  of  different 
functionals.  As  was  discussed  section  5.1.2,  an  objective  functional,  g ,  is  either  the  total 
noise  power  in  the  array  output,  referred  to  as  noise_pwr  in  awd,  or  the  maximum  value 
that  the  magnitude  of  the  pattern  function  attains  over  a  certain  region,  referred  to  as 
maxjnag_S. 

5.2.3  The  constraints  section 

This  section  is  used  to  specify  design  constraints.  As  was  mentioned  in  section  5.1.2,  a 
design  constraint  is  either  a  null  (null),  or  constraint  functional  with  certain  bounds  on  it. 
Constraint  functionals  can  be  the  real  part,  imaginary  part,  or  the  magnitude  (re,  im,  or 
mag)  of  an  element  weight  (W)  or  the  pattern  function  (S),  or  may  be  the  maximum  value 
that  the  magnitude  of  the  pattern  function  attains  over  a  certain  region  (max_mag_S). 
Functional  inequalities  are  allowed  to  take  the  forms 

g  <=  expr  ; 

<7  >=  expr  ; 
g  ==  expr  ; 
exprl  <=  g  <=  expr 2  ; 
abs( g)  <=  expr  ; 

where  g  is  one  of  the  above  functionals  (other  than  those  involvingmag,  null,  and  max_mag_S), 
and  expr ,  exprl  and  expr 2  represent  any  scalar  expressions. 

The  last  form 

abs  (  g )  <=  expr  ; 

is  equivalent  to 


-expr  <=  g  <=  expr  ; 

The  special  functionals  involving  mag,  and  max_mag_S,  can  only  be  used  in  the  form 
g  <=  expr  ; 

In  section  5.1.3,  we  discussed  how  these  constraints  are  approximated  by  linear  ones. 

5.2.4  Producing  and  reporting  awd  results 

Once  the  design  problem  has  been  successfully  compiled  into  a  mathematical  program, 
the  mathematical  program  is  solved  using  LSSOL  [GHM*S6].  LSSOL  usually  ends  in  one  of 
the  following  results  that  are  reported  by  awd, 
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•  found  strong  minimum:  the  unique  feasible  optimal  weight  values  ( w )  were  found. 

•  found  weak  minimum:  non-unique  feasible  weight  values  (10)  that  minimize  the  objec¬ 
tive  function  were  found. 

•  solution  appears  to  be  unbounded:  happens  usually  if  the  weight  values  (u?)  in 
the  AYVDSL  program  do  not  have  explicit  bounds.  It  is  generally  a  symptom  of  an 
ill-formulated  design  problem. 

•  no  feasible  point  found:  there  are  no  weight  values  (w)  that  satisfy  the  constraints. 
Some  (or  all)  of  the  constraints  in  the  constraints  section  must  be  relaxed. 

At  the  conclusion,  the  array  description,  along  with  the  final  designed  weights,  the  state 
of  all  constraints  and  the  final  objective  function  value  are  produced  and  saved  in  a  certain 
file.  Design  constraints  are  printed  in  that  file  in  the  same  order  in  which  they  were  listed 
in  the  source  file. 

The  listing  for  design  constraints  consists  of  six-column  lines  of  the  form 
constraint  value  lo_bnd  up_bnd  lambda  status. 

These  columns  are  as  follows. 

Column  1:  description  of  what  the  functional  is,  such  as  its  name,  e.g.  max_mag_S. 
Column  2:  the  value  of  the  functional. 

Column  3:  the  value  of  the  lower  bound  on  the  functional. 

Column  4:  the  value  of  the  upper  bound  on  the  functional. 

Column  5:  the  value  of  the  Lagrange  multiplier  for  the  functional. 

Column  6:  This  field  is  empty  if  the  functional  is  within  its  bounds,  or  the  upper  and  lower 
bounds  are  equal  (a  nonviolated  equality  constraint).  Otherwise,  it  will  contain  one  of: 

•  lb 

•  ub 

•  violates  lb 

•  violates  ub 

according  to  w'hether  the  functional  is  at  its  lower  or  upper  bound,  or  exceeds  its  lower 
or  upper  bound,  respectively. 

Example  6.1  contains  a  partial  listing  taken  from  an  actual  awd  run. 
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6  Examples 

6.1  A  Linear  Array  of  20-Elements 

In  this  section  it  is  desired  to  find  a  weight  vector  that  will  satisfy  the  constraints  exhibited 
in  Figure-7  for  a  linear  array  of  20  sensors  equally  spaced  along  the  y-axis  with  0.5A  spacing. 
The  objective  is  to  minimize  the  total  noise  power  in  the  received  signal.  Element  power  is 
assumed  to  be  unity  for  each  element. 


-100  -SO  -60  -40  -20  0  20  40  60  80  100 


Figure  7:  Specifications  for  array  pattern  function  design. 

A  possible  AWDSL  source  file  for  such  a  problem  looks  as  follows 

/  * 

A  linear  array  of  20  elements  spaced  at  0.5  wavelength  apart  along  the  y-axis. 
*/ 


/*  basic  definitions: 

*/ 

ftdefine  TARGET 

-15 

^define  NULL 

-55 

#def ine  BEAM.WIDTH 

10 

#def ine  NULL.WIDTH 

20 

((define  NULL_GAIN 

-40 
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#def ine  SIDE.GAIN 


-20 


array_description{ 

n_elements  =  20;  /*  the  number  of  antenna  elements  in  the  array  */ 

for  i  =  0  to  n_elements-l  do  { 

element_loc_y(i)  =  i*0.5  -4.75; 
element_loc_x(i)  =  0.0; 
element_pwr(i)  =  1.0; 

> 

> 


subject.to 

/*  target  angle  TARGET  as  defined  above  */ 
re (S (TARGET))  >=  1; 

/*  a  null  centered  at  NULL  degrees  and  of  width  NULL_WIDTH  */ 

max_mag_S (NULL-NULL_WIDTH ,  NULL+NULL_WIDTH ,  100)<=  inv_dB(NULL_GAlN) ; 

/♦no  sidelobe  has  gain  higher  than  SIDE_GAIN  dB  */ 

max_mag_S(-90,  TARGET-BEAM.WIDTH)  <=  inv_dB(SIDE_GAIN) ; 
max_mag_S (TARGET+BEAM_WIDTH  ,  90)  <=  inv_dB(SIDE_GAIN) ; 

} 

Awd  was  run  on  the  above  source  file  and  a  unique  optimal  weight  vector  was  obtained. 

The  following  plots  show  the  resulting  awd-optimai  pattern  and  how  it  relates  to  design 

constraints. 

A  partial  listing  of  the  awd  generated  output  looks  as  follows 
awd  version  2.0 


element  x-loc  y-loc  noise_pwr  complex  weight 
0  0  -4.75  1  0.02029  -0.0261412 

1  0  -4.25  1  0.02168  -0.0306461 

2  0  -3.75  1  0.0405  0.00488037 
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Magnitude  of  Array  Space  Factor  vs.  Angle 


Angle  of  Incidence  in  Degrees 


Figure  S:  Awd  optimal  design  for  Example  6.1.  Plot  of  magnitude  of  pattern  function  in 
dB  vs.  AOI. 


18  0  4.25  1  0.02168  0.0306461 

19  0  4.75  1  0.02029  0.0261412 

constraint  value  lo_bnd  up_bnd  lambda  status 

re(S(-15) )  1  1  le+10  0.1135  lb 


max_mag_S (-90,-25, 72)  0.09612  0  0.1  -0.01418  ub 

Above  constraint  attains  bounds  at: 

1  -25  0.09612  -0.01418 

max_mag_S (-5 ,90,106)  0.096120  0.1  -0.02477  ub 

Above  constraint  attains  bounds  at : 

1  -5  0.09612  -0.01396 

1  -0.5189  0.09612  -0.008679 

1  4.858  0.09612  -0.002128 

found  strong  minimum  in  998  iterations 
objective  function  value  =  0.0538599 
total  noise  power  =  0.0538599 
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Magnitude  of  Array  Space  Factor  vs  Angle 


Figure  9:  Awd  optimal  design  for  Example  6.1.  Polar  plot  of  magnitude  of  pattern  function 
in  dB. 


6.2  A  Parabolic  Array  of  '20-Elements 

In  this  section  it  is  desired  to  find  a  weight  vector  that  will  satisfy  the  constraints  exhibited 
in  Figure-7  for  the  parabolic  array  of  20  sensors  shown  in  Figure-10.  The  objective  is  still 
minimizing  the  total  noise  power  in  the  received  signal. 

In  addition  to  the  constraints  illustrated  in  the  previous  example,  a  constraint  on  the 
“back-of-array”  gain  was  included.  That  is,  the  magnitude  of  the  pattern  function  for  signals 
arriving  with  Angles  Of  Incidence  between  90°  and  270°  was  constrained.  This  constraint 
has  to  do  with  the  asymmetry  of  the  parabolic  array  as  opposed  to  the  linear  array.  The 
source  file  is  essentially  similar  to  that  for  the  linear  array  above,  and,  thus,  will  not  be 
listed.  It  is  yet  interesting  to  note  that  the  geometry  of  the  array  was  specified  by  placing 
the  statements 

for  i  =  0  to  n_elements-l  do  { 

element_loc_y(i)  =  i*0.5  -  4.75; 

element_loc_x(i)  =  (4. 75'2-element  loc  y(i)‘2)/4.75~2; 

> 

in  the  array-description  section,  and  that  the  extra  “back-of-array”  constraints  were  added 
by  simply  adding  the  line 
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Plot  of  Normalized  Element  Locations  for  Antenna  Array 


Figure  10:  A  Parabolic  Array  of  20-Elements. 


max_mag_S(90 ,  270)  <=  inv_dB(BACK_GAIN) ; 

to  the  constraints  section  of  the  source  file.  BACK-GAIN  was  defined  to  be  —  lOdB. 

Running  awd  on  the  source  file  for  this  problem,  a  unique  optimal  weight  vector  was 
obtained.  The  following  plots  show  the  resulting  awd-optimal  pattern  and  how  it  relates  to 
design  constraints. 

6.3  8-Element  Wideband  Linear  Array 

In  this  section,  awd  is  used  to  design  element  weights  for  a  linear  array  of  S  elements 
responding  to  nnr.owband  signals  with  known  center  frequencies.  The  elements  are  equally 
spaced  along  the  r-axis  with  a  1.8127  wavelength  spacing.  Using  ESPRIT  [RHR87],  it  was 
determined  that  there  were  two  incident  signals.  One,  considered  the  target  signal,  is  arriving 
at  an  angle  of  —9.39°  relative  to  broadside,  while  the  other,  the  interfering  signal,  is  arriving 
at  1 .26°. 

T  he  operational  frequency  band  of  interest  for  the  array  under  consideration  is  frequencies 
up  to  1.26kHz.  If  we  consider  this  frequency  as  the  nominal  frequency,  the  target  signal  will 
have  a  normalized  frequency  of  .75.  Thus,  it  was  determined  that  the  array  pattern  should 
have  at  least  unity  gain  for  signals  arriving  from  the  target  DOA  with  the  target  frequency. 
Furthermore,  it  was  determined  that  unwanted  signals  due  to  nonlinear  effects  in  the  sensors 
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Figure  11:  Awd  optimal  design  for  Example  6.2.  Plot  of  magnitude  of  pattern  function  in 
dB  vs.  AOI. 

were  present  at  frequencies  significantly  less  than  the  target  frequency.  Thus,  it  was  required 
that  any  signal  arriving  at  the  DOA  with  a  normalized  frequency  of  less  that  0.5  should 
encounter  at  least  — 20dB  suppression. 

As  for  the  interfering  signal,  it  was  determined  that,  regardless  of  its  frequency,  it  should 
encounter  — 40dB  suppression.  In  addition,  a  certain  band  of  possible  directions  of  arrival 
around  the  computed  interference  DOA  should  also  be  suppressed.  This  band  is  to  be  5° 
wide.  That  is  to  say,  for  any  signal  arriving  at  1.26°  ±  5°  relative  to  broadside,  with  any 
frequency  within  the  operational  frequency  band,  the  array  pattern  gain  should  be  less  than 
— 40dB. 

With  the  above  information  in  mind,  an  awd  source  program  was  written,  and  passed 
to  awd.  A  comparison  between  the  results  obtained  by  awd  and  those  of  ESPRIT  is  shown 
in  Figures  13,  14,  and  15. 
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Array  Gain  for  ESPRIT  and  AW  0  Optimal  Copy  Vectors 


Figure  13:  A  comparison  of  combined  array  pattern  obtained  using  the  ESPRIT  optimal 
weight  vector  and  that  obtained  using  Awd  for  signals  arriving  with  the  target  frequency  ol 
0.75. 
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Array  Gain  for  AWO  Optimal  Copy  Vector 


1.5 


Figure  14:  A  3-dimensional  plot  of  the  magnitude  of  the  pattern  function,  using  the  Awd 
optimal  weight  vector,  as  a  function  of  angle  and  normalized  frequency. 
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Figure  15:  A  3-dimensional  plot  of  the  magnitude  of  the  pattern  function,  using  the  ESPRIT 
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1  INTRODUCTION 


Awd  is  a  computer  aided  design  program  for  designing  antenna  element  complex  weights 
for  arbitrary  arrays  of  omnidirectional  elements  with  fixed  relative  locations  and  negligible 
mutual  interactions.  Elements  are  assumed  to  have  constant  complex  gain  over  a  prede¬ 
termined  (finite)  frequency  band  of  operation.  Currently,  signal  sources  are  assumed  to  be 
farfield.  An  array  configuration  for  the  case  when  all  elements  and  signals  lie  in  a  plane  is 
shown  in  Figure  1. 


y 


Figure  1:  Antenna  array  under  consideration. 


Awd  takes  as  its  input  an  antenna  weight  design  problem  written  in  a  specialized  an¬ 
tenna  weight  design  specification  language.  Awd  then  compiles  this  problem  into  a  convex 
mathematical  program.  The  resulting  mathematical  program  is  solved  using  a  convex  pro¬ 
gram  solver.  The  output  of  awd  is  a  set  of  complex  numbers,  one  per  antenna  element, 
that  constitutes  the  weight  vector.  These  complex  weights  determine  the  relative  scaling  and 
phase  shift  by  which  the  output  of  each  antenna  element  contributes  to  the  combined  array 
output.  This  combined  array  output  is  a  scalar  that  is  the  inner  product  of  the  awd  designed 
weight  vector  and  the  vector  output  of  the  array. 
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2  DESCRIPTION 


Awd  is  invoked  as  follows: 
awd  [-1]  filename 

where  filename  is  a  file  that  contains  a  specific  design  problem  specified  using  an  antenna 
weight  design  specification  language  that  is  discussed  below.  The  design  problem  is  then 
compiled  into  a  mathematical  program.  This  is  the  compilation  phase.  Note  that  in  the 
resulting  mathematical  program,  the  real  and  imaginary  parts  of  an  element  complex  weight 
are  considered  as  two  independent  real  optimization  variables. 

If  no  errors  occur  in  the  compilation,  an  optimization  routine  is  called  to  solve  the 
mathematical  program.  This  is  the  optimization  phase.  The  first  step  in  this  phase  is 
to  determine  the  feasibility  of  the  mathematical  program.  If  the  program  is  feasible,  optimal 
values  of  the  design  weights  are  found.  The  optimization  code  used  in  the  current  version  of 
awd  is  LSSOL  [GHM*86],  a  library  of  convex  programming  routines  for  linear,  quadratic, 
and  least  squares  optimization. 

Upon  completion,  the  antenna  array  specification,  the  weights  designed,  and  the  values 
of  the  constraint  functionals  and  their  upper  and  lower  bounds  are  written  into  the  file 
filename. out.  The  values  of  the  design  variables  are  written  into  the  file  filename,  w  if  they 
are  optimal,  or  into  the  file  filename. last  .w  if  the  optimization  failed. 

If  the  -1  option  is  specified,  the  log  information  from  the  convex  program  solver  LSSOL 
(the  optimization  code  used)  is  written  to  the  file  filename,  lslist.  Extra  options  for  LSSOL 
may  be  specified  in  the  file^aame.lsoptn. 

3  PREPARING  TO  USE  AWD 

The  following  information  concerning  the  design  problem  is  required  to  set  up  the  input  to 
awd: 

•  a  finite  operational  frequency  range, 

•  pre-steering  information, 

•  and,  antenna  element  locations. 

First  of  ail,  a  finite  operational  frequency  band  of  interest  is  selected  by  setting  lower  (/i0«er) 
and  upper  (/upper)  frequencies.  Currently,  antenna  element  responses  are  assumed  to  be 
constant  over  the  selected  band.  Note  that  for  narrowband  array  design  problems,  the  lower 
and  upper  frequencies  can  be  equal.  All  frequencies  are  specified  in  units  of  a  certain  nominal 
frequency,  /nom inai-  A  good  choice  of  /nomjnai  is  the  center  frequency  of  the  antenna  element 
operational  frequency  band.  A  frequency  of  interest  (/)  in  physical  units  (Hertz)  is  input  to 
awd  as  a  normalized  frequency,  /  =  ///nomjnai.  In  other  words,  a  frequency  of  0.5  specifies 
a  frequency  of  half  the  nominal  frequency. 
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If  the  array  is  pre-steered  to  a  certain  reference  direction  (say,  azimuthal  angle  0Te f,  and 
elevational  angle  <f>Te()  for  signals  of  a  certain  normalized  frequency  (say  /re f),  this  information 
should  be  made  available  to  awd. 

Finally,  locations  (in  rectangular  coordinates)  of  the  array  elements  must  be  given.  The 
positions  are  to  be  specified  in  units  of  Anominaj,  the  wavelength  corresponding  to  the  nominal 
operation  frequency  (/nom jnai)  and  the  underlying  speed  of  propagation  (c)  in  the  medium 

(^nominal  ^//nominal)  • 

It  is  assumed  that  the  designer  is  concerned  with  the  array  space  factor1  [RWD84,  Sch43j. 
This  is  a  function  of  azimuthal  angle,  0 ,  elevation  angle,  <f> ,  (normalized)  signal  frequency, 
/,  2  complex  weights,  W,-’s,  and  (normalized)  element  location  vectors,  P,  =  [xt-  y,  z,]T.  The 
function  is  given  by: 


n— 1 


S(0, 4>J)  =  'E  wt  exp  j(k  -  krc[)TP„ 
.'=0 


(1) 


where  n  is  the  number  of  elements  in  the  array, 


kTe[  —  27r/ref 


sin  d>nf  cos  0Ttj 
sin  <j>ref  sin  0Tt} 
cos  <f>nJ 


(2) 


and 


k  =  2ir  f 


sin  <j>  cos  0 
sin  <j>  sin  0 
cos  <f) 


(3) 


If  no  pre-steering  information  are  given,  kn(  will  be  assumed  to  be  a  zero  vector  (which  is 
equivalent  to  assuming  a  zero  reference  frequency). 

The  user  must  also  specify  the  total  noise  power  output  pt  of  each  element  in  the  array 
over  the  operational  frequency  band.  Each  p,  can  be  thought  of  as  the  integral  of  the  noise 
power  density  of  the  ith  element  over  the  operational  frequency  band.  This  noise  can  be  due 
to  thermal,  amplifier,  and/or  spatial  noise.  The  total  noise  power  in  the  combined  array 
output  is  given  by: 

w  =  X>,W-  (4) 

i=0 

This  assumes  that  noises  in  the  different  array  elements  are  independent. 

Finally,  the  user  must  determine  design  objectives  and  constraints.  In  the  current  version 
of  awd,  design  objectives  are  allowed  to  be  one  (or  a  weighted  sum)  of  the  following: 


•  max  max  max  |S(0,  <f>,  /) |,  maximum  value  that  the  magnitude  of  the  array 
/j  </</2  4>\  9\ 

space  factor  assumes  between  azimuth  angles  0X  and  02 ,  for  any  elevation  angle  be¬ 
tween  <f>i  and  fa,  and  any  (normalized)  frequency  between  fx  and  /2. 

'The  array  space  factor  is  sometimes  referred  to  as  the  array  directivity  pattern. 

2 From  now  on  it  is  assumed  that  all  frequencies  are  normalized  (/  =  ///nominal)  frequencies  unless 
otherwise  stated. 
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•  total  noise  power  in  the  combined  array  output  ( cf .  equation  (4)). 

Design  constraints  are  constructed  from  the  following  expressions: 

•  max  max  max  \S(0, 6,  /)  j ,  same  as  above. 

•  real,  imaginary,  or  magnitude  of  aS(0,  <f>,  /),  where  a  is  some  complex  weighting  factor, 
6  is  some  azimuthal  angle,  is  an  elevation  angle,  and  /  is  a  (normalized)  frequency. 

•  real,  imaginary,  or  magnitude  of  orW,-,  where  a  is  some  complex  weighting  factor,  and 
W{  is  the  complex  weight  factor  for  element  i. 

All  the  expressions  above  can  be  shown  to  be  convex  functionals  of  the  real  and  imaginary 
parts  of  the  weight  factors  W.-’s. 

Another  design  constraint  can  be  imposed  by  setting  S(0,<f>,f)  to  zero  at  given  incident 
angles,  6  and  (f>,  for  some  (normalized)  frequency,  /.  This  is  termed  placing  a  null  at  0 ,  (j), 
and  /.  It  can  be  shown  that  placing  a  null  at  angles,  9  and  (f>,  and  frequency,  /,  is  an  affine 
constraint  on  the  real  and  imaginary  parts  of  the  weight  factors  M^’s. 


4  THE  ANTENNA  WEIGHT  DESIGN  SPECIFICATION 
LANGUAGE  (AWDSL) 

The  Antenna  Weight  Design  Specification  Language  (AWDSL)  that  the  program  awd  com¬ 
piles  and  runs  consists  of  three  major  sections,  array  description,  objective  function,  and 
design  constraints .  Each  of  these  sections  consists,  as  will  be  detailed  below,  of  a  key  word, 
followed  by  a  left  brace  ‘{’,  followed  by  a  group  of  statements  (some  of  which  is  characteristic 
to  the  certain  section),  and  ended  by  a  right  brace 
A  typical  AWDSL  program  will  look  as  follows: 

array_description  { 

/*  description  of  array  geometry  and  element  noise  powers*/ 

> 

minimize  { 

/*  description  of  design  objective  */ 

} 

subject _to  { 

/*  description  of  design  constraints  */ 

> 

Each  of  these  sections  is  described  below  in  full  detail. 
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4.1  Statement  Syntax 

Syntax  is  similar  to  the  C  language:  statements  are  separated  by  semi-colons  (‘;’V  and 
groups  of  statements  are  enclosed  by  braces  ‘}’).  The  input  is  free  format:  any  number 
of  spaces,  newlines,  and  tabs  can  replace  a  space.  Awd  passes  the  source  file  filename 
(which  contains  the  AWDSL  program)  through  the  C  preprocessor  cpp ,  so  all  the  standard 
preprocessor  functions  are  available,  for  example  #def  ine  and  #include.  C-style  comments 
(enclosed  between  “/*”  and  “*/”)  are  permitted.  As  in  C,  comments  do  not  nest. 

A  statement  can  be  any  of  the  following. 

•  Multiple  statements:  as  in  C,  a  statement  can  be  replaced  by  any  number  of  state¬ 
ments  enclosed  in  braces  (‘{’,  ‘}’). 

•  Assignment:  real  variables  can  be  assigned  the  values  of  expressions.  See  below. 

•  For  loop:  looping  takes  the  following  forms: 

for  var  =  exprl  to  expr2  do  stmt 

for  var  -  exprl  to  expr2  step  expr3  do  stmt 

/*  the  word  ‘do*  is  synonymous  to  'sum,'  thus 
the  above  can  be  written  as:  */ 
for  var  =  exprl  to  expr2  sum  stmt 
for  var  =  exprl  to  expr2  step  expr3  sum  stmt 

The  variable  var  (a  lower  case  letter)  is  assigned  the  value  of  expression  exprl.  Expr2 
sets  the  final  value  of  var.  Expr3  sets  the  increment  of  var  each  time  around  the 
loop.  The  default  step  size  is  1.0.  The  expressions  expr2  and  expr3  are  only  evaluated 
when  the  for  loop  starts.  Thus,  variable  step  sizes  are  not  possible.  For  example,  the 
following  loops  execute  the  statement  stmt  10  times: 

for  i  =  0  to  9  do  stmt 

for  i  =  0  to  0.9  step  0.1  do  stmt 

•  Array-description  statement:  is  used  to  describe  the  array  geometry  and  character¬ 
istics.  Array-description  statements  can  only  appear  in  the  array  description  section. 
An  array-description  statement  can  be  one  of  the  following: 

n_elements  =  expr\ 
lower_freq  =  exprl; 
upper_freq  =  exprl; 
element  _loc_x(eipr^)  =  expr3; 
element  _Loc_y(expri?)  =  expr3; 
element  J.oc.z(expr2)  =  expr3; 
eleraent_pwr(erpr^)  =  cxpr3; 
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eleraent_w  (expr2)  =  (exprS,  exprJ^j', 
reference  =  (.expr5,expr6,expr  7); 
mag_2_lin  =  cxpr8 ; 

where  expr  is  the  number  of  antenna  elements  in  the  array,  exprl  is  any  non-negative 
valued  expression,  expr2  represents  a  given  element  index  and  must  evaluate  to  a  value 
between  0  and  n.elements  — 1,  and  exprS  and  exprj  are  any  real  valued  expressions. 
Expr5,  expr6,  and  expr7  give  the  azimuth,  elevation  and  frequency  information  of  the 
reference  signal  for  pre-steered  arrays.  Expr8  assumes  the  value  of  2  or  4  depending 
on  whether  it  is  wished  to  approximate  a  circle  by  a  square  or  an  octagon. 

•  Functional  statement:  can  be  either  a  constraint  statement ,  which  is  used  to  specify 
a  design  constraint  in  the  constraints  section  of  the  source  file,  or  an  objective  term , 
where  the  optimization  objective  is  the  sum  of  the  objective  terms  that  are  specified 
in  the  objective  section  of  the  source  file.  These  are  explained  below. 

4.2  Constants,  Variables  and  Expressions 

Numerical  constants  are  the  same  as  in  C.  The  mantissa  can  have  an  optional  fraction  part. 
An  exponent  can  be  specified  as  ‘e’  or  ‘E\  followed  by  an  optional  *+’  or  followed  by 
an  integer.  There  can  be  no  spaces,  tabs,  or  new-lines  in  a  constant.  “V.pi”  and  u,/,e”  are 
replaced  by  7r  (3.14...)  and  e  (2.718...)  respectively.  Also,  “infinity”  is  replaced  by  (l.OelO). 
In  addition,  since  angles  in  the  program  are  assumed  to  be  in  degrees,  two  constants  are 
made  available  to  the  user.  These  two  constants  are  “deg_2_rad”  and  “rad_2_deg,”  which 
are  replaced  by  (tt/  180)  and  (180/tt)  respectively. 

One  more  constant,  called  “correction-factor,”  can  only  be  used  in  the  objective  and 
constraints  sections,  but  not  in  the  array  description  section.  This  constant,  a s  will  be 
discussed  below,  has  to  do  with  the  fact  that  awd  approximates  magnitude  constraints,  which 
are  quadratic  in  nature,  by  a  number  linear  constraints.  The  constant  “correction-factor” 
is  internally  set  to  l/yxos(22.5)  if  the  magnitude  constraint  is  to  be  approximated  by  4  pairs 

of  linear  constraints,  and  to  l/^/cos(45)  if  the  magnitude  constraint  is  to  be  approximated 
by  2  pairs  of  linear  constraints. 

Variables  (which  can  take  on  only  scalar  floating-point  values)  are  restricted  to  names 
with  one  lower  case  letter.  The  variables  can  be  assigned  real  numbers  (doubles)  and  used 
in  scalar  expressions.  For  example, 

x  =  1  +  2  *  5.2; 
y  =  x  +  3; 

will  assign  11.4  to  x  and  14.4  to  y. 

The  binary  operators  *+’,  c*’,  ‘/\  (plus,  minus,  multiply,  divide  and  exponentiation) 

operate  between  two  scalar  expressions.  They  associate  left  to  right  and  have  the  usual 
precedence  (“°  highest,  followed  by  ‘/’  (equal),  followed  by  *+’,  (equal)).  Parentheses 
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can  be  used  in  the  usual  manner  to  force  order  of  evaluation.  For  example,  all  of  the  following 
statements  assign  8.0  to  the  variable  a. 

a=2*2+2*2; 
a  =  2  *  (1+3); 
a  =  0.5  *  -3; 
a  =  2  “  3; 
a  =  4  /  2  *  4; 

The  functions  sin,  cos  and  tan  return  the  corresponding  trigonometric  function  values 
of  degree  arguments.  The  functions  asin,  acos  and  atan  return  the  inverse  trigonometric 
function  values  (as  in  the  C  mathematical  library),  except  that  the  angles  are  in  degrees.  If 
atan  is  called  with  two  arguments  (comma  (V)  separated),  it  will  be  equivalent  to  the  C 
library  function  atan2  (with  the  result  converted  to  degrees,  resulting  in  an  answer  between 
—  180  and  180  degrees).  The  functions  exp,  log,  loglO,  dB,  inv_dB,  sqrt  and  abs  return 
the  exponential  and  logarithm  (base  e),  the  logarithm  (base  10),  the  decibel  (201og10)  and 
its  inverse,  positive  square  root  and  absolute  values  of  their  arguments  (respectively).  In 
addition,  the  functions  min  and  max  take  two  arguments  (comma  (‘,’)  separated)  and  return 
the  smaller  or  larger  argument  respectively.  All  of  these  functions  take  expressions  which 
evaluate  to  real  scalars  as  arguments.  Finally,  floor  finds  the  greatest  integer  less  than  the 
given  argument.  For  example,  all  of  the  following  statements  assign  a  number  close  to  2.0 
to  the  variable  b. 

b  =  sqrt (4) ; 
b  =  loglO(lOO) ; 

b  =  2  *  min(5  *  cos(le-lO),  exp(max(0*x,-l))) ; 
b  =  f loor(2 . 5) ; 
w  =  rad_2_deg  *  '/.pi; 

b  =  2  *  (cos(w/4)  *  cos(w/4)  +  sin(w/4)  *  sin(w/4)); 

Expressions  in  awd  can  be  used  in  place  of  a  constant  anywhere  in  the  source  file  filename. 
Thus,  for  example,  expressions  can  be  used  for  array  indices  (the  compiler  will  round  indices 
to  the  nearest  integer). 

4.3  The  Array  Description  Section 

This  section  must  appear  first  in  the  program.  As  the  name  suggests,  this  section  contains 
information  about  the  array  structure.  The  number  of  elements  (n.elements)  in  the  array  is 
the  first  thing  to  be  specified.  Optionally,  the  upper  (upper_freq)  and  lower  (lowerjfreq) 
bounds  on  the  operational  frequency  band  are  specified  next.  If  either  upper_freq  or 
lower_freq  is  not  specified,  it  will  be  defaulted  to  unity  (the  nominal  frequency).  Thus 
for  narrowband  elements,  the  user  can  ommit  the  specification  of  such  bounds  and  they  will 
be  defaulted  to  the  nominal  frequency.  The  array  geometry  is  then  given  by  specifying  the 
x,  y,  and  2  locations  (in  rectangular  coordinates)  of  the  array  elements  in  units  of  nominal 
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wavelength,  Anominai.  Finally,  the  noise  power  of  each  element,  p,  is  specified.  Optionally,  the 
initial  clement  weights  can  be  specified  in  this  section.  The  above  information  is  set  using 
array-description  statements  of  the  form  discussed  above. 

If  the  array  is  pre-steered  to  a  certain  direction  for  a  given  (scaled)  signal  frequency, 
such  (optional)  information  should  also  be  made  available.  If  no  pre-steering  information 
are  given,  kre{  will  be  assumed  to  be  a  zero  vector.  That  is,  an  infinite  reference  frequency 
will  be  assumed.  An  array  description  section  thus  has  the  form: 

array_descri.pt  ion{ 

/*  specify  the  number  of  antenna  elements  in  the  array  */ 
n_elements  =  number  of  elements ; 

/*  specify  the  pre-steering  information  */ 

reference  =  ( reference  azimuth, reference  elevation,  reference  frequency ); 

/*  how  to  approximate  magnitude  constraints*/ 
mag_2_lin  =  2  for  a  square  or  4  for  an  octagon; 

/*  specify  the  lower  and  upper  bounds  on  operational  frequency  band  */ 
lower_freq  =  lower  bound; 
upper-freq  =  upper  bound; 

/*  specify  the  x  locations  */ 
element_x_loc(0)  *  x-loc  of  element  0; 

element_x_loc(n_elements  -  1)  =  x-loc  of  element  n_elements-/; 

/*  specify  the  y  locations  */ 
element_y_loc(0)  =  y-loc  of  element  0; 

element-y_loc(n_elements  -  1)  =  y-loc  of  element  n.elements-/ ; 

/*  specify  the  z  locations  */ 
element-z_loc(0)  =  z-loc  of  element  0; 


element-Z_loc(n_elements  -  1)  =  z-loc  of  element  n_elements-/; 

/*  specify  the  element  noise  powers  */ 
element-pwr(O)  =  power  of  element  0; 


element_pwr(n_elements  -  1)  =  power  of  clement  n_elements-/; 
/*  specify  the  initial  weights  */ 
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element_w(0)  -  ( real  part  ,  imaginary  part ) ; 


} 


element_w(n_elements  -  1) 


(.real  part  ,  imaginary  part)  ; 


Note  that  the  array  elements  are  indexed  from  0  to  n_elements  — 1. 

The  mag_2_lin  (optional)  value  specifies  to  the  program  the  method  of  approximating  a 
magnitude  constraint.  A  value  of  2  means  that  the  magnitude  constraint  (which  resembles 
a  circle)  is  to  be  approximated  by  two  pairs  of  linear  constraints  (a  square),  while  a  value 
of  4,  the  default,  indicates  that  the  magnitude  constraints  should  be  approximated  by  four 
pairs  of  linear  constraints  (an  octagon). 

The  element-w  statements  above  allow  an  initial  estimate  of  the  solution  to  be  given.  If 
awd  is  being  re-run  on  a  source  file  which  has  had  a  small  number  of  constraints  changed,  it 
is  likely  that  the  old  optimal  weight  values  will  be  a  good  starting  point  for  the  optimization 
phase.  This  will  speed  up  the  optimization  phase.  If  element.w  is  not  specified,  a  zero 
starting  point  will  be  used. 

The  user  is  advised  to  store  the  above  coefficients  in  separate  files  that  can  be  conveniently 
included  in  the  file  filename  using  the  “#include”  facility  of  the  C  preprocessor.  For  example, 
if  the  coefficients  for  x-locations  are  in  the  file  “problem .  x,”  then  the  following  should  appear 
in  the  source  file: 


#include  "problem. x" 

It  should  be  noted  that  the  C  preprocessor  (cpp)  requires  that  the  symbol  ‘#’  is  the  first 
character  on  the  line. 

The  following  is  a  typical  array  description  section,  which  describes  a  circular  array  of 
ten  narrowband  elements: 

array_description  { 

/*  the  number  of  elements  is  the  first  thing  to  declare  */ 
n_elements  =  10; 

/*  specify  a  circular  array  of  radius  equal  to  half  wavelength  */ 
for  i  =  0  to  n_elements  -1  do  { 

element_x_loc(i)  =  0 . 5*cos (360/n_elements  *  i) ; 
element_y_loc(i)  =  0 . 5*sin(360/n_elements  *  i) ; 
element_z_loc(i)  =  0; 

} 

/*  the  other  information  is  included  in  separate  files  */ 

#include  problem. p 
•include  problem. w 
> 
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The  file  problem. p  contains  the  noise  power  of  the  array  elements,  and  the  file  problem. w 
contains  initial  starting  weights.  Note  that,  as  in  C,  to  loop  over  two  or  more  statements, 
the  statements  must  be  enclosed  by  braces. 

4.4  The  Objective  Section 

This  is  the  second  section  of  the  AWDSL  program  and  is  optional.  If  the  objective  section 
is  missing,  awd  will  understand  that  the  problem  is  a  feasibility  problem  and  simply  try  to 
satisfy  the  constraints  in  the  constraints  section  (discussed  below),  provided  the  latter  section 
exists.  The  objective  function  is  formed  by  summing  all  the  objective  terms  appearing  in 
the  objective  section. 

An  objective  term  is  of  the  form: 

expr  *  objective-functional ; 

The  optional  scalar  expression  expr  allows  relative  weighting  of  the  different  functionals.  The 
objective-functional  must  be  one  of  the  following: 

•  noise_pwr:  total  noise  power  in  the  comoined  array  output  (as  discussed  above). 

•  max_mag_S (0 -range)  (4>-range)  (/ -range):  the  maximum  value  that  the  magnitude  of 
the  space  factor  assumes  on  certain  points  within  the  given  azimuthal  ( 0 ),  elevation 

and  frequency  (/)  ranges.  If  one  of  the  ranges  is  omitted,  it  is  replaced  by  a 
single  default  value.  The  default  value  for  the  azimuth  is  0°,  the  value  for  the  elevation 
is  90°  (parallel  to  the  x-y  plane),  and  the  default  value  for  the  frequency  range  is  1 
(corresponding  to  the  nominal  frequency). 

Each  range  consists  of  a  lower  bound,  upper  bound,  and  an  integer  resembling  the 
number  of  points  at  which  the  given  range  is  to  be  sampled.  The  number  of  interme¬ 
diate  sample  points  is  optional.  There  are  no  default  values  for  the  lower  and  upper 
bounds  for  any  of  the  above  ranges,  but  the  value  for  the  number  of  sample  points 
is  optional  and  is  defaulted  to  100  points  for  the  azimuth  range,  50  for  the  elevation 
range,  and  10  for  the  frequency  range.  Such  values  are  arbitrary  and  thus  the  user 
should  avoid  using  such  defaults. 

The  following  is  a  sample  objective  section  which  explains  the  use  of  the  above  functionals: 

minimize  { 

/*  total  output  noise  */ 
noise_pwr ; 

/*  magnitude  of  space  factor  on  a  given  azimuth  angle  range  */ 
max_mag_S(l5 ,  345); 

/*  the  above  is  equivalent  to 
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max_mag_S(15 ,  345,100)()(); 

or 

max_mag_S(l5 ,  345 , 100) (90) (1) ; 

*/ 

/*  magnitude  of  space  factor  for  given  azimuthal  and  elevational  ranges  */ 
max_mag_S(-75,  -30,  40) (40,  140,  40); 

/*  the  above  is  equivalent  to 

max_mag_S(-75,  -30,  40) (40,  140,  40)(); 

or 

max_mag_S (-75 ,  -30,  40) (40,  140,  40) (1); 

*/ 

/*  magnitude  of  space  factor  for  given  azimuthal  and  frequency  ranges  */ 
max_mag_S(75 ,  105,  40)0(0.8,  1.2,  20); 

/*  the  above  is  equivalent  to 

max_mag_S(75 ,  105,  40) (90) (0.8,  1.2,  20); 

*/ 

> 

4.5  The  Constraints  Section 

This  section  is  used  to  specify  design  constraints.  This  is  the  third  and  final  part  of  the  pro¬ 
gram.  If  no  such  section  exists,  the  problem  will  be  an  unconstrained  optimization  problem. 

The  following  functionals  may  be  referred  to  in  the  constraints  section: 

•  re(sc/_/ctr*W(0),  im(sc/_/ctr*W(0),  mag(sc/_/ctr*W(z)):  the  real,  imaginary,  and 
magnitude  of  the  complex  weight  of  element  i. 

•  null(0,  <f> ,  /):  the  space  factor  has  a  null  at  azimuth  angle  0,  elevation  angle  <f>,  and 
(normalized)  frequency  /,  i.e.  S(0,  <j>,  f)  =  0.  If  /  is  omitted,  it  will  be  assumed  to 
have  unity  value  (the  nominal  frequency).  If  <f>  is  omitted,  it  will  be  defaulted  to  90°, 
that  is  a  plane  parallel  ot  the  x-y  one.  The  default  for  0  is  0°. 

•  re(scLfctr*S(0 ,  <j>,  /)),  im(scLfctr*S  (0 ,  <j>,  /)),  mag(sc/_/ctr*S(0,  <f> ,  /)):  the 
real,  imaginary,  and  magnitude  of  the  space  factor  with  respect  to  incident  azimuth 
angle  6 ,  elevation  angle  and  (normalized)  frequency  /.  If  /  is  omitted,  it  will 
be  assumed  to  have  unity  value  (the  nominal  frequency).  If  (f>  is  omitted,  it  will  be 
defaulted  to  90°,  that  is  a  plane  parallel  ot  the  x-y  one.  The  default  for  0  is  0°. 

•  max_mag_S  {0-range)  (<f>-rangc)  {/-range) :  same  as  the  functional  discussed  in  the  ob¬ 
jective  section. 
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In  the  above,  scLfctr  is  a  complex  (optional)  scale  factor  that  has  the  form  ( real  part, 
imaginary  part).  If  no  scale  factor  is  given,  it  is  assumed  to  be  (1.0,  0.0),  i.e  1. 

Functional  inequalities  are  allowed  to  take  the  forms 

g  <=  expr  ; 
expr  >=  g  ; 
g  >=  expr  ; 
expr  <=  g  ; 
g  ==  expr  ; 
expr  ==  g  ; 

expr l  <=  g  <=  expr2  ; 
exprl  >=  g  >=  expr2  ; 
abs  (  g )  <=  expr  ; 
expr  >-  abs  (  g  )  ; 

where  g  is  one  of  the  above  functionals  (other  than  those  involving  mag,  null,  and  max_mag_S), 
and  expr,  exprl  and  expr2  represent  any  scalar  expressions. 

The  last  two  forms 

abs  (  g )  <=  expr  ; 

and 


expr  >-  abs (p)  ; 
are  equivalent  to 

-expr  <=  g  <=  expr  ; 

The  special  functionals  involving  mag,  and  maxjnagJS,  can  only  be  used  in  the  form 
g  <=  expr  ; 

where  expr  is  a  non-negative  valued  expression. 

A  mag  functional  inequality  is  replaced  by  the  four  inequalities 

-  expr/ correct  ion.f  actor  <=  re  <=  expr/ corrections  actor ; 

-expr/correction_f actor  <=  im  <=  expr/correctionSactor; 

-  expr/ correct  ion_f  actor  <=  (re  +  im)/sqrt(2)  <=  expr/correction_f  actor ; 
-expr/correction-f actor  <=  (re  -  ira)/sqrt(2)  <=  expr/correction-factor; 

where  “correctionSactor”  is  a  constant  that  is  set  to  l/^/cos(22.5),  and  the  above  ap¬ 
proximation  yields  a  maximum  of  ±0.34  dB  error  if  the  mag_2_lin  parameter  is  set  to  4  (the 
default).3 

A  sample  constraints  section  is  given  below: 

3If  “correction-factor”  is  set  to  l/>/cos(45)  if  mag_2_lin  is  set  to  2,  and  the  program  uses  only  the 
first  two  inequalities  in  approximating  a  magnitude  constraint. 
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subject_to  { 

/*  bounds  on  W  coefficients  */ 

for  k=0  to  n_elements-l  do 
mag(W(k) )  <=  1.0; 

/*  a  null  at  90  degrees  incident  angle  and  nominal  frequency  */ 
null (90) ; 

/*  the  above  null  can  be  rewritten  as  */ 
re(S(90))  ==  0.0; 
im(S(90) )  ==  0.0; 

/*  or  it  can  be  rewritten  as  */ 
mag(S(90))  <=  0.0; 

/*  constrain  magnitude  of  space  factor  on  given  azimuth  and  frequency  ranges  */ 
max_mag_S(60 ,  300,  100) (90) (0.8,  1.2,  5)  <=  0.1; 

/*  require  that  the  array  space  factor  has  a  magnitude  higher 

them  one  in  the  zero  direction  for  a  range  of  frequencies  */ 
for  f  =  0.7  to  1.3  step  0.05  do 
re(S(0,90,f ))  >=  1.0; 

> 

5  OPTIMIZATION 

Once  the  design  problem  has  been  compiled  into  a  mathematical  program,  the  compilation 
phase  is  complete  and  the  optimization  phase  can  be  executed.  This  is  where  the  mathemati¬ 
cal  program  is  solved.  Using  LSSOL  [GHM*86],  a  convex  quadratic  program  solver  from  the 
Systems  Optimization  Laboratory  of  the  Department  of  Operations  Research  at  Stanford 
University,  the  objective  function  is  minimized  subject  to  the  affine  constraints  given. 

The  optimization  phase  usually  ends  in  one  of  the  following  results: 

•  found  strong  minimum:  the  unique  feasible  optimal  weight  values  (VU)  were  found. 

•  found  weak:  minimum:  non-unique  feasible  weight  values  (VU)  that  minimize  the  ob¬ 
jective  function  were  found. 

•  solution  appears  to  be  unbounded:  happens  only  if  the  weight  values  (H/)  in  the 
AWDSL  program  do  not  have  explicit  bounds.  It  is  generally  a  symptom  of  an  ill- 
formulated  design  problem. 

•  no  feasible  point  found:  there  are  no  weight  values  (W)  that  satisfy  the  con¬ 
straints.  Some  (or  all)  of  the  constraints  in  the  constraints  section  must  be  relaxed. 


IX- 1 4 


With  the  first  two  outcomes,  the  optimal  weight  values  (W),  are  written  into  the  file 
filename. v.  This  file  is  written  in  the  same  format  as  required  in  the  array  description 
section.  Thus,  W  values  can  be  conveniently  used  as  the  starting  point  for  the  optimization 
in  a  subsequent  run  of  awd  by  placing  the  line 

•include  "filename. w" 

in  the  array  description  section.  This  will  make  the  optimization  phase  of  a  subsequent  run 
of  awd  run  faster  (if  the  constraints  or  objectives  have  not  been  changed  significantly). 

With  the  remaining  outcomes  of  the  optimization  phase,  the  last  W  values  are  written 
into  the  file  filename. last  .v  instead.  This  allows  the  values  to  be  inspected  without  the 
result  of  any  successful  run  (in  filename,  v)  being  destroyed. 

Regardless  of  the  outcome  of  the  optimization  phase ,  the  array  description,  along  with 
the  produced  weights,  the  state  of  all  constraints  and  the  final  objective  function  value  are 
written  into  the  file  filename. out.  The  constraint  functionals  are  printed  in  that  file  in  the 
same  order  in  which  they  were  listed  in  the  source  file. 

As  far  as  the  file  filename. out  is  concerned,  the  constraint  functionals  can  be  divided 
into  two  categories;  those  for  which  only  one  line  in  the  file  filename. out  is  printed,  and 
those  with  one,  or  more  printout  lines  (depending  on  the  results  of  the  optimization). 

For  the  first  type  of  constraint  functionals,  the  line  is  formatted  into  six  (6)  columns 
separated  by  spaces.  A  description  of  these  columns  follows. 

Column  1:  one  of  the  following  functionals 

•  re(sc/_/cfr*W(0 ) , 

•  im  (  scLfctr  *W  (i)  ) , 

•  mag  (scLfctr*  W(i)), 

•  re(sc/_/c<r*S(0 ,  <f>,  /)), 

•  im(scLfctr*S(9 ,  /)), 

•  mag  ( scLfctr  *S(0 ,  <j>,  /)), 

•  null (5,  (ft,  /). 

If  <j)  and/or  /  are  not  printed,  it  is  implied  that  they  assume  the  default  values  of  90° 
(the  horizontal  plane)  and  unity  (the  nominal  frequency)  respectively. 

Column  2:  the  value  of  the  functional. 

Column  3:  the  value  of  the  lower  bound  on  the  functional. 

Column  4:  the  value  of  the  upper  bound  on  the  functional. 

Column  5:  the  value  of  the  Lagrange  multiplier  for  the  functional. 
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Column  6:  This  field  is  empty  if  the  functional  is  within  its  bounds,  or  the  upper  and  lower 
bounds  are  equal  (an  unviolated  equality  constraint).  Otherwise,  it  will  contain  one 
of: 


•  lb 

•  ub 

•  violates  lb 

•  violates  ub 


according  to  whether  the  functional  is  at  its  lower  or  upper  bound,  or  exceeds  its  lower 
or  upper  bound,  respectively. 

For  example,  the  following  fragment  from  the  constraints  section  of  the  source  file  filename 

re(S(0))  >=  1.0; 
mag(S(0) )  <=  1.5; 
null (270, 90, 0.9); 


may  produce  the  following  lines 

in  filename,  out 

re(S(0)) 

1 

1 

le+10 

1.04 

lb 

mag(S(0) ) 

1.082 

0 

1.5 

0 

null (270, 90, 0.9) 

0 

0 

0 

0 

The  second  type  of  constraint  functionals  includes  the  functional  max_mag_S(0-ran<7e)  (<j>- 
range)  (/ -range).  For  this  functional,  the  first  line  is  formatted  exactly  as  in  the  case  above. 
If  the  bounds  are  met  or  violated,  more  lines  will  be  printed.  The  first  extra  line  would  be 
cither, 


Above  const,  mt  attains  bounds  at: 


or, 

Above  constraint  violates  bounds  at: 

as  appropriate.  The  lines  that  follow  would  be  five-column-printouts  with: 

Column  1:  the  azimuthal  angle  at  which  the  functional  attained  or  violated  its  constraints. 

Column  2:  the  elevational  angle  at  which  the  functional  attained  or  violated  its  constraints. 

Column  3:  the  (normalized)  frequency  at  which  the  functional  attained  or  violated  its 
constraints. 

Column  4:  the  value  of  the  functional  at  that  point. 

Column  5:  the  Lagrange  multiplier  of  the  functional  at  that  point. 
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For  example,  the  following  fragment  from  the  constraints  section  of  the  source  file  filename 
max_mag_S (90 ,270,50) (90) (1.1)  <=  0.1; 
may  produce  the  following  lines  in  filename .  out 

max_mag_S (90 , 270 , 50)  0.1035  0  0.1  1.87  ub 

Above  constraint  attains  bounds  at: 

90  90  1.1  0.1035  1.24 

180  90  1.1  0.0999  0.63 

It  is  only  fair  to  warn  the  user  that,  due  to  the  nature  of  optimization  algorithms  in 
general,  it  may  be  the  case  that  the  algorithm  did  not  detect  all  of  the  constraints  at  which 
the  bounds  are  attained.  Thus,  in  the  above  example,  it  may  be  the  case  that  the  bounds 
are  also  attained  at  other  angles  not  shown  above.  However,  all  violations  are  detected. 


6  FILE  SUMMARY 

Input  files 
filename 
filename  .lsoptn 

Output  files 
filename  .v 
filename,  last  .w 
filename,  out 
filename,  lslist 

7  HINTS 

7.1  LSSOL  Options 

There  are  many  optional  parameters  that  the  user  can  specify  in  the  file  filename,  lsoptn, 
[GHM*86].  In  particular,  crash  tolerance  and  feasibility  tolerance  have  significant 
effects  on  the  number  of  iterations  that  LSSOL  takes  to  find  the  optimal  solution.  Thus,  it 
may  be  unwise  to  let  the  optimization  phase  of  the  program  run  for  long  periods  of  time, 
since  there  may  be  an  option  that  can  be  set  in  the  lsoptn  file  to  make  the  program  run 
faster  on  your  specific  problem.  See  the  LSSOL  manual  for  all  the  allowable  options  and 
their  effects  on  the  behavior  of  the  code. 

Options  in  the  filename,  lsoptn  file  can  be  set  with  statements  of  the  form: 

#  this  is  a  sample  lsoptn  file 
crash  tolerance  =  1.0e-4 


the  source  file  name  given  as  argument  to  awd. 
extra  LSSOL  options  (optional). 

W  values  after  successful  optimization. 

W  values  after  unsuccessful  optimization, 
listing  file  for  variable  and  constraint  values. 
LSSOL  listing  if  -1  option  given. 
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feasibility  tolerance  =  1.0e-6 
feasibility  phase  iteration  limit  =  100 
optimality  phase  iteration  limit  =  100 
print  level  =  5 

Lines  starting  with  the  character  ‘#’  are  treated  as  comments  and  are  ignored. 


7.2  Problem  Size 

The  user  should  note  that  awd  may  generate  a  large  number  of  linear  constraints.  Thus, 
it  is  advisable  to  be  economical  with  the  size  of  problems  handled  using  the  code.  This 
limitation  is  due  to  LSSOL  (see  [GHM*86]). 


7.3  Speed 

Most  of  the  number  crunching  is  done  in  the  optimization  code  LSSOL.  LSSOL  uses  a 
standard  set  of  “BLAS”  (Basic  Linear  Algebra  Subroutines).  Those  can  be  freely  replaced 
by  any  appropriate  BLAS  fine  tuned  to  the  user’s  machine.  In  particular,  by  profiling  the 
code,  we  found  out  that  most  of  the  time  is  spent  in  two  specific  BLAS,  “DAXPY”  and 
“DDOT.”  Running  awd  on  the  “narrowband”  example  shown  below,  it  was  found  that 
99.4%  of  the  time  was  spent  in  solving  the  mathematical  program,  39%  of  the  total  time  was 
spent  in  “DAXPY,”  and  20%  of  the  time  was  spent  in  “DDOT.” 

7.4  Assuring  Sensible  Answers 

Although  awd  will  generate  correct  results  for  any  problem  of  reasonable  size,  the  designed 
weights  need  not  be  sensible  from  an  engineering  point  of  view.  Awd  is  just  an  optimization 
code,  and  correct  mathematical  programming  results  do  not  necessarily  imply  acceptable 
results.  The  results  will  be  acceptable  only  if  the  problem  is  correctly  formulated.  The  user 
should  keep  in  mind  that  a  mathematical  program  solver  like  LSSOL  will  minimize  the  given 
objective  function  subject  to  the  given  constraints.  Thus,  any  feature  of  concern  must  be 
included  either  as  a  constraint  or  as  an  objective.  Otherwise,  the  code  will  not  “care”  about 
such  a  feature.  Ignoring  some  features  or  merely  specifying  a  feasibility  problem  may  very 
well  yield  an  array  space  factor  that  is  not  “acceptable.” 


7.5  Awd  Error  Handling 

All  syntax  errors  are  caught,  but  the  resulting  error  messages  are  sometimes  hard  to  connect 
to  the  actual  error  in  the  source  code.  No  effort  was  made  to  accurately  report  syntax  errors 
or  recover  from  them.  If  awd  produces  a  long  stream  of  syntax  errors,  look  only  at  the  first 
one  that  it  prints.  Fix  this  and  then  re-run  awd. 
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8  MATLAB 


PRO-MATLAB  is  not  a  requirement  for  using  awd,  but  it  may  simplify  specifying  the 
array  structure  and  simulating  the  produced  results.  Thus,  some  supporting  C-programs 
and  PRO-MATLAB  scripts  are  available  to  the  user.  In  the  following,  we  will  detail  the 
description  of  each  of  these  scripts. 

8.1  Transforming  from  PRO-MATLAB  to  Awd 

MATLAB  is  a  convenient  tool  for  setting  the  array  description.  The  following  support 
utilities  are  available  for  converting  MATLAB  format  into  ASCII  format  acceptable  to  awd. 

•  m2awd :  is  a  C-program  that  takes  a  character  string,  say  “filename”  as  its  argument, 
searches  for  variables  x,  y,  z,  p ,  and  re  in  a  file  called  “filename. mat,”  and  writes 
whichever  of  these  variables  it  finds  to  corresponding  ASCII  files  called  “ filename. x,” 
“ filename. y,”  “ filename. z,”  “filename. p,”  and  “ filename. w”  respectively.  These  files 
then  contain  the  x,  y,  and  z  locations  of  the  antenna  array  elements,  the  noise  power 
of  the  elements  p,  and  the  weight  factors  w  written  in  the  same  format  that  awd 
accepts  in  its  array  description  section.  These  ASCII  files  can  be  conveniently  included 
in  the  source  code  of  the  AWDSL  program  using  the  “tinclude”  facility  of  the  C 
preprocessor. 

•  maL2.awd:  is  a  PRO-MATLAB  script  that  can  be  run  while  using  MATLAB  and  it 
allows  the  transformation  from  MATLAB  to  awd  to  be  handled  internally.  This  script 
looks  for  a  MATLAB  ASCII  variable  called  “f  ile_name”  whose  value  is,  say,  ‘ filename'. 
The  script  then  saves  the  x,  y,  z,  p,  and  w  existing  variables  to  “filename. mat,”  and 
runs  m2awd  to  produce  the  required  ASCII  files. 

8.2  Transforming  from  Awd  to  PRO-MATLAB 

The  following  utilities  for  converting  from  awd  output  format  to  MATLAB  format  are  also 
available. 

•  awd2m:  is  a  C-program  that  takes  a  character  string,  say  “filename,”  as  its  argument, 

searches  for  ASCII  files  called  “filename,  x,”  “filename. y,”  “filename,  z,”  “filename. p,” 
and  “ filename. w.”  These  files  must  contain  the  x,  y,  and  z  locations  of  the  antenna 
array  elements,  the  noise  power  of  the  elements  p,  and  the  weight  factors  w  written  in 
the  same  format  that  awd  accepts  in  its  array  description  section.  It  then  stores  the 
values  corresponding  to  whichever  of  these  files  it  finds  as  variables  x,  y,  z,  p,  and  w, 
respectively,  in  a  file  called  “ filename.  This  file  can  be  then  conveniently  loaded 

to  MATLAB. 

•  awd-2-mat:  is  a  PRO-MATLAB  script  that  can  be  run  while  using  MATLAB.  This 
script  allows  the  transformation  from  MATLAB  to  awd  to  be  handled  internally. 
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Figure  2:  Sample  layouttool  window 


This  script  looks  for  a  MATLAB  ASCII  variable  called  “f  ile_name”  whose  value  is, 
say,  ‘filename'.  The  script  then  runs  the  program  awd2m  to  create  a  file  called  “file¬ 
name.  mat”  that  contains  the  x,  y,  z,  p ,  and  w  variables.  Awd.-2.mat  then  loads  the 
above  “ filename. mat”  file  to  MATLAB. 


8.3  Graphical  Layout  of  Array  Elements 

The  vector  handling  capabilities  built  into  Matlab  provide  a  convenient  means  of  specifying 
regularly  shaped  antenna  arrays  (e.g.,  linear  arrays,  circular  arrays,  or  cylindrical  arrays). 
However,  experimentation  with  more  irregular  array  shapes  can  be  cumbersome.  The  fol¬ 
lowing  support  utilities  are  available  to  design  antenna  arrays  interactively: 

•  layouttool :  is  a  C-program  that  takes  a  character  string,  say  “filename,"  as  its  argu¬ 
ment,  and  searches  for  an  ASCII  file  by  that  name.  This  file  must  contain  the  number 
of  array  elements,  the  x,  y ,  and  2  locations,  the  noise  power  p,  and  the  weight  factors 
tv  of  the  antenna  array  elements,  written  in  the  same  format  that  awd  accepts  in  its 
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array  description  section.  If  no  input  file  argument  is  given,  layouttool  will  start  up 
with  an  empty  array,  layouttool  starts  up  a  Suntools  window,  providing  the  following 
functionality: 

—  placing  new  elements:  click  the  left  mouse  button  on  a  free  space  inside  the 
drawing  area 

—  moving  elements:  click  at  an  element  with  the  left  mouse  button,  move  the  mouse 
to  the  desired  location,  and  release  the  mouse  button 

—  deleting  elements:  click  the  middle  mouse  button  at  the  element  to  be  deleted 

—  editing  element  parameters  numerically,  click  the  right  mouse  button  at  the  el¬ 
ement  to  view  a  pop-up  window  containing  the  element  parameters,  fill  in  the 
appropriate  fields,  and  acknowledge  by  clicking  the  ok  button 

—  loading  an  antenna  array  specification:  fill  in  the  dir  and  file  fields  and  click  the 
load  button 

—  saving  an  antenna  array  specification:  fill  in  the  dir  and  file  fields  and  click  the 
save  button 

—  clear  all  elements:  click  the  clear  button 

—  redraw  screen  (in  case  the  display  got  corrupted):  click  the  redraw  button 

—  increase  resolution:  click  the  zoom  in  button 

—  decrease  resolution:  click  the  zoom  out  button 

—  quit  layouttool:  click  the  quit  button  (make  sure  to  save  the  array  specification 
before  quitting) 

—  scroll  to  the  right:  click  at  the  scale  to  the  left  of  the  drawing  area 

—  scroll  to  the  left:  click  at  the  scale  to  the  right  of  the  drawing  area 

—  scroll  down:  click  at  the  scale  at  the  top  of  the  drawing  area 

—  scroll  up:  click  at  the  scale  at  the  bottom  of  the  drawing  area 

•  layout:  is  a  PRO-MATLAB  script  that  can  be  run  while  using  MATLAB.  This  script 
saves  the  Matlab  variables  n ,  x,  y ,  z,  p,  and  w  (number  of  elements,  z,  z/,  2  positions, 
noise  power,  and  element  weights,  respectively)  into  the  file  “array.mat”,  converts  it 
into  an  awd  file  “array.awd”,  and  starts  up  the  layouttool  with  that  file  as  an  in¬ 
put.  Upon  quitting  the  layouttool ,  the  layout  script  converts  “array.awd”  back  into 
“array.mat”,  and  loads  that  file  into  MATLAB. 

8.4  Simulating  Awd  Results  with  PRO-MATLAB 

The  following  functions  are  available: 
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•  gets[w,x,y,z,0  ,(f>,f  ,kTetJ:  is  a  PRO-MATLAB  function  that  can  be  run  while  using 
MATLAB.  This  function  calculates  the  antenna  array  space  factor  S  (as  described 
above)  at  azimuth  angle  6  (in  degrees),  elevation  angle  <f>  (in  degrees)  and  (normalized) 
frequency  /  with  weight  factors  w  and  geometry  described  by  x,  y,  and  2.  Pre-steering 
information  (if  any)  is  given  in  kTe f.  If  this  last  input  is  omitted  it  will  be  defaulted  to 
zero.  The  result  will  be  a  complex  number. 

Note  that  0,  <j>,  and  /  can  be  vectors.  In  this  case,  the  output  of  the  function  will  be 
a  (complex)  matrix  of  length(/)*length((^)  rows  and  length($)  columns.  The  first 
length(/)  rows  will  contain  function  values  obtained  using  the  first  element  of  <j>,  and 
so  on. 

•  db-plot[what,w,x,y,z,0 ,kTe(]:  is  a  PRO-MATLAB  function  that  can  be  run  while 
using  MATLAB.  This  function  assumes  that  6  and  <f>  are  given  in  degrees.  The  function 
first  calls  geLs[w,x,y,z,6 ,<f>,f  ,kTe(]  to  calculate  the  antenna  array  space  factor  S. 

If  “what”  is  set  to  1,  /  should  be  a  scalar.  If  “what”  is  set  to  2,  0  should  be  a  scalar. 
Otherwise,  <f>  should  be  a  scalar.  The  function  evaluates  the  absolute  value  of  the  S 
matrix  in  decibels  and  plots  it  as  a  three-dimensional  object  (using  PRO-MATLAB’s 
mesh  function).  If  “what”  is  set  to  1,  such  object  will  be  a  plot  of  the  magnitude  of  the 
array  in  decibels  versus  0  and  <j>.  If  “what”  is  set  to  2,  such  object  will  be  a  plot  of  the 
magnitude  of  the  array  in  decibels  versus  /  and  <j>.  Otherwise,  the  three-dimensional 
object  will  be  a  plot  of  the  magnitude  of  the  array  in  decibels  versus  /  and  6. 

If  one  (but  not  both)  of  the  vectors  that  should  enter  the  3-D  plot  is  actually  a  scalar, 
a  two  dimensional  plot  of  the  magnitude  of  S  in  decibels  as  a  function  of  the  remaining 
vector  entry  will  result. 

The  function  will  return  the  output  of  get_s. 

•  polar-plot[w,x,y,z,0,<f>,f  ,kre{]:  is  a  PRO-MATLAB  function  that  can  be  run  while  using 
MATLAB.  This  function  requires  that  (f>  and  /  be  scalars.  The  function  calls  gets  to 
calculate  the  value  of  the  space  factor  S.  The  function  then  plots  the  magnitude  of 
the  space  factor  in  polar  coordinates. 

•  polar.db_plotfw,x,y,z,0 ,<f>,f ,krc[J:  is  a  PRO-MATLAB  function  that  can  be  run  while 
using  MATLAB.  This  function  is  the  same  as  polar.plot  except  that  the  plots  will  be  in 
decibels.  It  should  be  noted  here  that  the  value  of  the  magnitude  of  5,  will  be  scaled 
to  overcome  certain  difficulties  arising  from  the  use  of  MATLAB. 

•  plot-array (x,y,z,  what):  is  a  PRO-MATLAB  function  that  can  be  run  while  using 
MATLAB.  If  “what”  is  set  to  1,  the  function  plots  the  projection  of  the  array  geometry 
on  the  x-y  plane.  If  “what”  is  set  to  2,  the  the  projection  of  the  array  geometry  on 
the  x-z  plane  will  be  plotted.  Finally,  if  “what”  is  set  to  3,  the  function  plots  the 
projection  of  the  array  geometry  on  the  y-z  plane. 
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9  EXAMPLES 


9.1  Linear  Array  of  Narrowband  Antenna  Elements 

In  this  section,  awd  is  used  to  solve  an  antenna  design  example  for  a  linear  array  of  20 
narrowband  elements  equally  spaced  along  the  y- axis  between  0  and  9.5A.  The  array  is 
assumed  to  lie  in  a  plane  parallel  to  the  x  —  y  plane  at  z  =  2A.  It  is  assumed  that  the 
element  noise  power  is  unity  for  all  elements. 

The  design  objective  is  to  minimize  the  combined  array  noise  power.  The  array  space 
factor  is  to  have  at  least  a  unit  gain  at  —15  degrees,  with  no  sidelobe  gain  of  more  than  —30 
dB  outside  a  20  degree  sector  centered  at  the  target  angle.  It  is  also  desired  to  have  a  —50 
dB  null  in  the  40  degree  sector  from  —75  to  —35  degrees. 

A  source  file  for  such  a  program  is  as  follows.  Note  that  the  (normalized)  frequency  index 
is  omitted  since  array  elements  are  assumed  to  be  narrowband. 

/* 

Linear  array  of  20  elements  spaced  at  0.5  wavelength  along  the  y-axis. 

*/ 


/*  basic 

definitions:  */ 

#def ine 

TARGET 

-15 

#def ine 

NULL 

-55 

#def ine 

BEAM_WIDTH 

20 

#def ine 

NULL_WIDTH 

40 

#def ine 

NULL_GAIN 

-50 

#def ine 

SIDE.GAIN 

-30 

array_description{ 

n_elements  =20;  /*  the  number  of  antenna  elements  in  the  array  */ 

for  i  =  0  to  n_elements-l  do  { 
element_loc_x(i)  =  0.0; 
element_loc_y(i)  =  i*0.5; 
element_loc_z(i)  =  2; 
element_pwr(i)  =  1.0; 

} 

> 


minimize! 

noise_pwr ; 


> 
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subject_to 

{ 

/*  target  angle  TARGET  as  defined  above  */ 
re(S(TARGET))  >=  1; 

/*  a  null  centered  at  NULL  degrees  and  of  width  NULL.WIDTH  */ 

raax_mag_S (NULL-NULL.WIDTH/2 , NULL+NULL_WIDTH/2 , 100)  <=  inv_dB (NULL.GAIN) 

/*  no  sidelobe  has  gain  hi<-  r  than  SIDE_GAIN  dB  */ 

max_mag_S(-90,  TARGET-BEAM_WIDTH/2)  <=  inv_dB(SIDE_GAIN) ; 
max_mag_S (TARGET+BEAM_WIDTH/2 ,  90)  <=  inv  dB(SIDE  GAIN); 

} 


The  .  out  file  for  such  a  program  looks  as  follows. 


awd  version  3.0 
copyright  Prometheus  Inc.  1988 


element 

x-loc 

y-loc 

z-loc 

noise_pwr 

complex  weight 

0 

0 

0 

2 

1 

0.01317 

0.008598 

1 

0 

0.5 

2 

1 

0.01626 

0.01105 

2 

0 

1 

2 

1 

-0.003668 

0.02777 

3 

0 

1.5 

2 

1 

-0.02633 

0 . 02438 

4 

0 

2 

2 

1 

-0.04747 

-0.004941 

5 

0 

2.5 

2 

1 

-0.03337 

-0.04614 

6 

0 

3 

2 

1 

0.01016 

-0.06534 

7 

0 

3.5 

2 

1 

0.06185 

-0.04119 

8 

0 

4 

2 

1 

0.07566 

0.01723 

9 

0 

4.5 

2 

1 

0.04208 

0.06901 

10 

0 

5 

2 

1 

-0.02257 

0.07739 

11 

0 

5.5 

2 

1 

-0.06809 

0.03611 

12 

0 

6 

2 

1 

-0.06931 

-0.02297 

13 

0 

'6.5 

2 

1 

-0.02624 

ID 

O 

O 

14 

0 

7 

2 

1 

0.02151 

-0.05211 

15 

0 

7.5 

2 

1 

0.04692 

-0.01678 

16 

0 

8 

2 

1 

0.03172 

0.01725 

17 

0 

8.5 

2 

1 

0.008442 

0.02679 

18 

0 

9 

2 

1 

-0.01457 

0.01529 

19 

constraint 

0 

9.5 

2 

value 

IX-24 

1 

lo_bnd  up. 

-0.01105  0.01084 

_bnd  lambda  sta* 

re(S(-15) ) 

1 

1 

le+10 

0.1296 

max_mag_S(-75,-35 , 100) (90) (1) 

0.00329 

0 

0.003162 

-0.5443 

Above 

constraint 

attains 

bounds  at : 

-69.4 

90 

1 

0.003217 

-0.02565 

-68.6 

90 

1 

0.003241 

-0.009777 

-59.8 

90 

1 

0.003071 

-0.03465 

-59.4 

90 

1 

0.003122 

-0.02704 

-51.4 

90 

1 

0.003173 

-0.07514 

-50.6 

90 

1 

0.003138 

-0.01113 

-44.6 

90 

1 

0.003107 

-0.09142 

-43.8 

90 

1 

0.003176 

-0.0215 

-38.6 

90 

1 

0.003175 

-0.1272 

-38.2 

90 

1 

0.003228 

-0.01943 

-35 

90 

1 

0.00329 

-0..014 

max_mag_S(-90, -25 , 100) (90) (1) 

0.03272 

0 

0.03162 

-0.03585 

Above 

constraint 

attains 

bounds  at : 

-27.6 

90 

1 

0.03046 

-0.01309 

-26.95 

90 

1 

0.03272 

-0.005512 

-26.3 

90 

1 

0.03042 

-0.01725 

max_mag_S(-5,90, 100) (90) (1) 

0.0329 

0 

0.03162 

-0.1475 

Above 

constraint 

attains 

bounds  at : 

-5 

90 

1 

0.03054 

-0.006633 

-4.05 

90 

1 

0.03224 

-0.0285 

-0.25 

90 

1 

0.03176 

-0.01261 

0.7 

90 

1 

0.03113 

-0.01891 

5.45 

90 

1 

0.0329 

-0.02123 

11.15 

90 

1 

0.03171 

-0.01002 

12.1 

90 

1 

0.0304 

-0.002987 

16.85 

90 

1 

0.03081 

-0.00645 

17.8 

90 

1 

0.0325 

-0.00292 

23.5 

90 

1 

0.03077 

-0.006671 

30.15 

90 

1 

0.03079 

-0.004516 

31.1 

90 

1 

0.03134 

-0.0007145 

37.75 

90 

1 

0.03052 

-0.004529 

46.3 

90 

1 

0.03211 

-0.003496 

47.25 

90 

1 

0.03065 

-0.001102 

55.8 

90 

1 

0.03109 

-0.00217 

57.7 

90 

1 

0.03041 

-0.003682 

71.95 

90 

1 

0.03047 

-0.005647 

72.9 

90 

1 

0.03049 

-0.004698 

found  strong  minimum  in  751  iterations 
objective  function  value  =  0.061177 
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total  noise  power  =  0.0611725 

The  next  three  graphs  contain  plots  of  the  array  geometry,  the  magnitude  of  the  array 
space  factor  in  polar  coordinates  and  the  magnitude  of  the  array  space  factor  in  dB  versus 
the  incident  angle  The  resulting  main  lobe  3  dB  beam  width  is  7°. 


Plot  of  Normalized  Element  Locations  for  Antenna  Array 


x-axis 


Figure  3:  PRO-MATLAB  plot  of  array  element  locations  for  example  9.1. 


9.2  Linear  Array  of  Wideband  Antenna  Elements 

In  this  section,  awd  is  used  to  solve  an  antenna  design  example  for  a  linear  array  of  20 
wideband  elements  equally  spaced  along  the  y-axis  between  0  and  9.5A.  The  elements  lie  in 
a  plane  parallel  to  the  x  —  y  plane  at  z  =  0.5A.  Taking  the  nominal  frequency  as  the  center 
frequency,  the  (normalized)  frequency  band  of  operation  is  assumed  to  be  [0.9,  1.1].  It  is 
assumed  that  the  element  noise  power  is  unity  for  all  elements. 

The  design  objective  is  to  minimize  the  combined  array  noise  power.  Signals  of  concern 
are  those  that  lie  in  the  plane  of  the  array.  The  array  space  factor  is  to  have  at  least  unit 
gain  at  all  band  frequencies  in  the  —15  degrees  direction,  with  no  sidelobe  gain  (at  any 
frequency)  of  more  than  —20  dB  outside  a  20  degree  sector  centered  at  the  target  angle.  It 
is  also  desired  to  have  a  —40  dB  null  in  the  40  degree  sector  from  —75  to  —35  degrees  for  all 
frequencies. 

A  source  file  for  such  a  program  is  as  follows. 

/* 
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Mag(S)  in  dB 


Magnitude  of  Array  Space  Factor  vs.  Angle 


Figure  4:  PRO-MATLAB  polar  plot  of  magnitude  of  array  space  factor  for  example  9.1 


Magnitude  of  Array  Space  Factor  vs.  Angle 


Angle  of  Incidence  in  Degrees 


Figure  5:  PRO-MATLAB  plot  of  magnitude  of  array  space  factor  in  dB  for  example  9.1 


Linear  array  of  20  elements  spaced  at  0.5  wavelength  along  the  y-axis. 
♦/ 


/*  basic  definitions:  */ 


#def ine 

TARGET 

-15 

/*  target  direction 

♦/ 

#def ine 

NULL 

-55 

/*  null  direction 

*/ 

#def ine 

BEAM.W 

20 

/♦  beam  width 

*/ 

#def ine 

NULL.W 

40 

/*  null  sector  width 

*/ 

#def ine 

NULL_G 

-40 

/♦  null  gain 

*/ 

#def ine 

SIDE_G 

-20 

/*  side  lobe  gain 

♦/ 

#def ine 

L.FREQ 

0.9 

/♦  lower  frequency 

*/ 

#def ine 

U_FREQ 

1.1 

/♦  upper  frequency 

*/ 

array_descript ion{ 


n_elements  =20;  /*  the  number  of  antenna  elements  in  the  array  ♦/ 

lower_freq  =  L_FREQ;  /*  lower  bound  on  operational  frequency  band  */ 
upper_freq  =  U_FREQ;  /*  upper  bound  on  operational  frequency  band  */ 

for  i  =  0  to  n_elements-l  do  { 
element_loc_x(i)  =  0.0; 
element_loc_y(i)  =  i*0.5; 
element_loc_z(i)  =  0.5; 
element _pwr(i)  =  1.0; 

> 

> 

minimize{ 

noise_pwr; 

} 

subject_to 

{ 

/♦at  least  unit  gain  for  all  frequencies  at  angle  TARGET  ♦/ 
for  f  =  lower_freq  to  upper.freq  step  0.02  do 
re (S (TARGET, ,f))  >=  1; 

/♦a  null  centered  at  NULL  degrees  and  of  width  NULL_W  for  all  frequencies*/ 
max_mag_S(NULL-NULL_W/2 ,NULL+NULL_W/2 , 20) () (lower_freq,upper_f req, 10) 

<=  inv_dB(NULL_G) ; 
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/*  no  sidelobe  has  gain  higher  than  SIDE_G  dB  for  any  frequency  */ 

max_mag_S(-90 ,  TARGET-BEAM_W/2 ,  40) () (lower_f req,upper_f req, 10) 

<=  inv_dB(SIDE_G) ; 

max_mag_S(TARGET+BEAM_W/2 ,  90,  40) () (lover_f req,upper_f req, 10) 

<=  inv_dB(SIDE_G) ; 

> 


The  resulting  .out  file  for  the  above  source  program  looks  as  follows. 


avd  version  3.0 
copyright  Prometheus  Inc.  1988 


element 

x-loc 

y-loc 

z-loc 

noise_pwr 

compl 

ex  weight 

0 

0 

0 

0.5 

1 

0.02534 

0 . 04824 

1 

0 

0.5 

0.5 

1 

0.04649 

0.0596 

2 

0 

1 

0.5 

1 

-0.01651 

0.05755 

3 

0 

1.5 

0.5 

1 

-0.06919 

0.04503 

4 

0 

2 

0.5 

« 

X 

-0.08411 

-0.02889 

5 

0 

2.5 

0.5 

1 

-0 . 04948 

-0.08397 

6 

0 

3 

0.5 

1 

0.03232 

-0.09897 

7 

0 

3.5 

0.5 

1 

0.08869 

-0.05132 

8 

0 

4 

0.5 

1 

0.1017 

0.02769 

9 

0 

4.5 

0.5 

1 

0.05035 

0.08539 

10 

0 

5 

0.5 

1 

-0.0223 

0.09507 

11 

0 

5.5 

0.5 

1 

-0.0767 

0.04564 

12 

0 

6 

0.5 

1 

-0.07906 

-0.01906 

13 

0 

6.5 

0.5 

1 

-0.03457 

-0.06201 

14 

0 

7 

0.5 

1 

0.01689 

-0.05426 

15 

0 

7.5 

0.5 

1 

0.04123 

-0.01935 

16 

0 

8 

0.5 

1 

0.03069 

0.01361 

17 

0 

8.5 

0.5 

1 

0.008726 

0.0176 

18 

0 

9 

0.5 

1 

-0.01201 

0.01062 

19 

0 

9.5 

0.5 

1 

-0.007842 

0.00767 

constraint 

value 

lo_bnd 

up_bnd 

lambda 

status 

re(S(-15,90,0.9)) 

1 

1 

le+10 

0.1168 

lb 

re(S(-15 , 90 ,0 . 92) ) 

1.121 

1 

le+10 

0 

re (S(- 15, 90, 0.94)) 

1.22 

1 

le+10 

0 

re(S(-15,90,0.96)) 

1.293 

1 

le+10 

0 

re(S(-15,90,0.98)) 

1.338 

1 

le+10 

0 

re(S(-15) ) 

1.352 

1 

le+10 

0 

re(S(-15 ,90 , 1 . 02) ) 

1.336 

1 

le+10 

0 

re(S(-15 ,90 , 1 . 04) ) 

1.289 

1 

le+10 

0 

re(S(-15 ,90 , 1 . 06) ) 

1.216 

1 

le+10 

0 

re(S(-15 ,90, 1 . 08) ) 

1.118 

1 

le+10 

0 

re(S(-15 ,90, 1.1)) 

1 

1 

le+10 

0.1533 

max.mag.S (-75 , -35,20) (90) (0.9,1 

•1,10) 

0.0104 

0 

0.01 

-0.7577 

Above 

constraint  attains 

bounds  at : 

-73 

90 

0.9 

0.01037 

-0.01158 

-37 

90 

0.9 

0.01034 

-0.07331 

-35 

90 

0.9 

0.01006 

-0.06621 

-35 

90 

0.92 

0.009767 

-0.03517 

-67 

90 

0.94 

0.01038 

-0.06528 

-53 

90 

0.96 

0.009714 

-0.05736 

-39 

90 

0.96 

0.0104 

-0.07326 

-45 

90 

0.98 

0.009965 

-0.04545 

-49 

90 

1.02 

0.009718 

-0.03688 

-65 

90 

1.04 

0.009613 

-0.07113 

-73 

90 

1.06 

0.009742 

-0.04784 

-35 

90 

1.06 

0.0104 

-0.06974 

-71 

90 

1.08 

0.009749 

-0.008774 

-75 

90 

1.1 

0.0104 

-0.02528 

-39 

90 

1.1 

0.009966 

-0.07045 

max_mag_S (-90 , -25 ,40) (90) (0.9,1 

.1,10)  0.104 

0 

0.1 

-0.267 

Above 

constraint  attains 

bounds  at : 

-25 

90 

0.9 

0.104 

-0.1223 

-28.25 

90 

0.92 

0.09642 

-0.1448 

max_mag_S(-5 ,90,40) (90) (0 .9 , 1 . 1 

,10) 

0.103 

0 

0.1 

-0.07177 

Above 

constraint  attains 

bounds  at : 

-5 

90 

0.9 

0.103 

-0.01779 

6.875 

90 

0.9 

0.09613 

-0.002469 

18.75 

90 

0.94 

0.09612 

-0.001434 

-0.25 

90 

0.96 

0.09653 

-0.007831 

11.63 

90 

0.98 

0.09614 

-0.00393 

-5 

90 

1.1 

0.09781 

-0.02107 

-0.25 

90 

1.1 

0.09653 

-0.01062 

4.5 

90 

1.1 

0.09654 

-0.006615 

found  strong  minimum  in  1319  iterations 
objective  function  value  =  0.115121 
total  noise  power  =  0.115121 


lb 

ub 


ub 


ub 


Figure  5  shows  the  magnitude  of  the  array  space  factor  (in  dB)  as  a  function  of  the 
(normalized)  frequency  and  the  incident  angle  (in  degrees). 
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Figure  6:  3-dimensional  PRO-MATLAB  plot  of  magnitude  of  array  space  factor  in  dB  for 
example  9.2. 

9.3  Cylindrical  Array  of  Narrowband  Elements 

In  this  example  we  show  how  awd  can  be  used  to  design  weights  for  a  3-D  array  with 
narrowband  elements.  We  assume  that  we  have  a  cylindrical  array  of  80  elements.  Cylinder 
height  is  4A  and  its  principal  axis  is  the  x-axis.  The  radius  of  the  cylinder  is  A.  The  sensors 
are  grouped  in  circles  of  10  elements  each.  Circles  lie  in  the  y-z  plane  and  are  spaced  0.5A 
a.  rt  along  the  x-axis.  Elements  are  spaced  uniformly  on  each  circle.  The  following  two 
figures  (Fig.  6  and  Fig.  7)  illustrate  the  geometry  of  this  array. 

Target  signals  arrive  with  an  azimuthal  AOA  of  15°  and  an  elevation  of  100°  (t.e.  10° 
below  the  x-y  plane).  It  is  desired  to  constrain  all  sidelobes  outside  an  azimuthal  band  of 
width  50°  and  an  elevation  band  of  width  15°  centered  around  the  target  directions.  Sidelobes 
gains  are  to  be  below  —  12dB. 

A  possible  awd  source  file  for  such  a  problem  is  as  follows. 

/* 

A  cylindrical  array  of  80  elements. 

*/ 

#def ine  N_X_ELT  8 

#def ine  N_RADIAL_ELT  10 

#def ine  HEIGHT  4 
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Projection  of  Antenna  Array  in  the  X-Y  plane 


x-axis 


Figure  7:  Output  of  ploLarray:  PRO-MATLAB  plot  of  array  element  location  projected  on 
the  x-y  plane  for  example  9.3. 


#def ine 

SIDE_GAIN 

-12 

#def ine 

B_WIDTH_A 

50 

#def ine 

B_WIDTH_E 

15 

#def ine 

TARGET.A 

15 

#def ine 

TARGET.E 

100 

#def ine 

N.AZIM 

140 

#def ine 

N.ELEV 

25 

array _description{ 


n_eleraents  =80;  /*  the  number  of  antenna  elements  in  the  array  */ 
mag_2_l in=2 ; 

for  i=0  to  N_X_ELT  -  1  do  { 

x  =  i*HEIGHT/N_X_ELT ; 
k  =  i*N_RADIAL_ELT ; 
for  j=  0  to  N_RADIAL_ELT-1  do  { 
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Figure  8:  Output  of  ploLarray:  PRO-MATLAB  plot  of  array  element  location  projected  on 
the  y-z  plane  for  example  9.3. 

t  =  j*360/N_RADIAL_ELT; 
element _pwr(k)  =  1.0; 
element_loc_x(k)  =  x; 
eleraent_loc_y(k)  =  cos(t); 
element_loc_z(k)  =  sin(t); 
k  =  k+ 1 ; 

> 

} 

•include  "ell.w" 

> 

minimize{ 

noise.pwr ; 

> 

subject_to 

{ 

/*  at  least  unity  gain  for  target  direction  */ 
re(S(TARGET_A ,  TARGET_E) )  ==  1; 
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/*  constrain  every  sidelobe  outside  a  sector  of  band  widths 

B_WIDTH_A  for  azimuth  and  B_WIDTH_E  for  elevation  */ 
max_mag_S(TARGET_A+B_WIDTH_A/2,TARGET_A+360-B_WIDTH_A/2,N_AZIM) 
(0 ,TARGET_E-B_WIDTH_E/2 ,N_ELEV)  <=  inv.dB(SIDE.GAIN) ; 
max_mag_S (TARGET. A+B_WIDTH_A/2 , TARGET. A+360-B_WIDTH_ A/2 ,N_AZIM) 
(TARGET_E+B_WIDTH_E/2 , 180.N.ELEV)  <=  inv.dB(SIDE.GAIN) ; 

/*  help  the  optimization  code  by  setting  some  (reasonable)  nulls  */ 
for  a=105  to  285  step  20  do  { 

for  t  =  0  to  60  step  30  do 
null(a,t) ; 

for  t  =  140  to  180  step  20  do 
null(a,t) ; 

} 

> 


Plots  of  the  awd  designed  array  space  factor  follows  in  Fig.  8  and  Fig.  9. 


Plot  of  Mag(S)  vs.  Azimuthal  and  Elevationai  AOI 


Figure  9:  3-dimensional  PRO-MATLAB  plot  of  magnitude  of  array  space  factor  for  example 
9.3 


IX-34 


Magnitude  of  Array  Space  Factor  vs.  Azimuthal  Angle 


Figure  10:  PRO-MATLAB  polar  plot  of  magnitude  of  array  space  factor  vs.  azimuth  for 
target  elevation  angle  for  example  9.3 
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Professional  Consultants  To  The  Scientific  Community 


Mr.  Benjamin  Cray 
Code  2133 
NUSC 

New  London,  CT  06320 


February  1,  1989 


Dear  Ben, 

Enclosed  are  the  results  we  discussed  over  the  phone  yesterday.  Figure  1  indicates  that  I  was  able 
to  reproduce  your  results.  Based  on  the  back-lobe  presence,  some  of  the  AWD  results  shown  minimized 
max  |S(0,/)|  over  the  region  [—167.5°,  —12.5°]  U  [12.5°,  167.5°].  It  should  be  obvious  from  the  AWD  plots 
when  this  was  done.  The  mesh  plot  in  Figure  2  was  a  complex  weight  design  with  /o  =  170Hz.  Figures  3a 
and  3b  show  the  results  of  optimizing  over  4  frequency  intervals  as  indicated.  Interestingly,  the  design  at 
the  highest  frequency  interval  was  significantly  poorer  than  the  others.  1  don’t  have  a  real  good  intuitive 
explanation  for  this  other  than  to  note  that  the  pre-steering  was  done  at  240Hz  for  all  4  designs.  I  am 
currently  working  on  designs  where  the  steering  is  done  at  the  mid-point  of  each  interval.  The  results  will 
follow  as  they  become  available. 

Note  well  that  the  plots  over  frequency  ranges  are  plots  of 

S(0 ,  /)  =  £  f 

■=i 

where  m  =  34  and  do  =  0°  for  the  cases  shown,  and  /o  is  indicated  in  the  figures.  This  implicitly  assumes 
the  analytic  or  complex  signal  is  available  (e.g.  via  Hilbert  transform).  If,  as  I  am  guessing,  the  current 
processor  configuration  simply  delays  the  real  received  signal  by  7o/4  where  To  =  l//o  to  generate  the 
quadrature  component,  the  plots  contained  herein  will  bear  little  (or  no)  resemblance  to  those  you  may 
have  generated  with  your  simulation  program  that  correctly  models  the  processor  structure.  A  simple  delay 
generates  the  quadrature  component  for  only  a  single  frequency  component  of  the  received  signal.  The 
overall  consequence  of  such  a  processing  strategy  is  to  create  a  frequency  dependent  array  manifold  from  a 
frequency  independent  one.  Since  AWD  is  not  currently  capable  of  handling  frequency  dependent  response 
functions  yet,  ... 

Hope  this  helps.  Please  let  me  know  how  it  is  received  at  NUSC.  Also,  I’d  appreciate  finding  out  who 
to  contact  about  the  upcoming  experiment  at  Seneca  Lake. 


Sincerly, 


Richard  Roy, 
Senior  Scientist 
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INC. 


Mr.  Benjamin  Cray 
Code  2133 
NUSC 

New  London,  CT  06320 


February  2,  1989 


re:  AWD  and  the  horseshoe  array 


Dear  Ben, 

...  and  here  are  the  results  for  centered,  design  over  the  four  bands.  Note  the  improvement!  Again  the 
results  assume  the  analytic  signal  is  available  over  the  band. 


More  as  it  becomes  available, 


Richard  Roy, 
Senior  Scientist 


end:  more  AWD  results  for  34  element  horseshoe  array 
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Appendix 

Miscellaneous  Information 


PART  I  -  THE  SCHEDULE 


SECTION  B  -  SUPPLIES  OR  SERVICES  AND  PRICES/COSTS 


0001  RESEARCH 

The  contractor  shall  furnish  the  level  of  effort  specified  In  Section  F, 
together  with  all  related  services,  facilities,  supplies  and  materials  needed 
to  conduct  the  research  and  prepare  the  reports  described  below.  The  research 
shall  be  conducted  and  reports  delivered  In  accordance  with  Section  F. 


0001 AA 

1.  Develop  the  Direct  Adaptive  Antenna  System,  Incorporating  analytic  null 
steering  methods  drived  during  Phase  I  Into  existing  adaptive  algorithms. 

2.  Produce  notch  filters  with  more  that  one  notch,  similar  to  Ideal 
single-notch  filters  developed  In  Phase  I. 

3.  Construct  and  code  wavelets  (harmonics)  for  high  frequency  high  resolution 
expansion  of  signal  arising  In  time-varying  environments. 

4.  Develop  a  convex  programming  approach  to  the  problem  of  choosing  antenna 
array  weights. 

5.  Consider  the  development  of  an  alternative  to  the  Cox  efficiency 
constraint  for  reducing  the  effect  of  shading  coefficient  errors. 

6.  Design  and  build  a  prototype  compiler  which  will  translate  specifications 
stated  In  a  natural  language  for  antenna  engineers  (PSL)  Into  executable  code 
In  the  form  of  a  convex  program  or  an  appropriate  linear  program. 

7.  Develop  a  Problem  Scenario  Generator  to  produce  Inputs  to  the  executable 
code  generated  by  the  compiler.  This  will  Include  a  Signal  Generator  and 
Operators  capable  of  combining  various  signals. 

8.  Develop  an  Antenna  Data  Analyzer  to  generate  plots  Illustrating  the  output 
of  the  executable  code. 

9.  Develop  a  user  Interface  to  enable  the  above  components  to  be  employed  in 
an  Iterative  and  Interactive  manner. 

10.  Expand  upon  the  new  peak  factor  results  obtained  In  Phase  I  by 
considering  the  general  problem  of  determining,  for  a  given  amplitude 
spectrum,  how  to  choose  the  phases  of  Fourier  coefficients  In  order  to 
minimize  the  peak  factor  In  the  corresponding  Inverse  Fourier  transform. 
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